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Tandem locomotives, teamed to pull a train : 
up a grade, are often called a double header. 
It would take a “35 header’ to equal the 
power of the RMI rocket engine shown here. 


- RMI propulsion systems are now being per- RMI since 1941 . . . reliability in application. 

a fected for use in ballistic and guided missiles Another reason why RMI’s program of 
as well as in the manned aircraft of the advanced research and engineering contin- 
future. These advanced propulsion systems ues to be a good investment for taxpayers, 
incorporate the principle which has guided government and industry! 
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Recent Advances in Real Gas Effects in Hypersonic Flow 


WAYLAND C. GRIFFITH 


Lockheed Missile Systems Division, Palo Alto, Calif. 


The author studied engineering 
science at Harvard University 
where he received an A.B. in 1946 
and a Ph.D. in 1949. His thesis 
was on an experimental measure- 
ment of vibrational relaxation 
times in gases. In 1950 he joined 
the shock wave laboratory of the 
Physics Department at Princeton 
University where he continued his 
research on physical problems in 
gas dynamics. In 1951 he was 
appointed assistant professor of 
physics; later, associate professor. 
The subjects on which he has worked extensively include tran- 
sonic flow, thermal boundary layers, shock structure, chemical 
kinetics, diffraction of blast waves, and supersonic flow. As one 
of the editors of the Journal of Fluid Mechanics, he has an inti- 
mate and up-to-date knowledge of the problem on which he writes. 
Dr. Griffith is currently manager of the Flight Sciences Division 
in the Missile Systems Division of Lockheed Aircraft Corp. 


F THE many problems confronting the aerodynamicist 
engaged in the development of satellites and ballistic 
missiles, two arise through departures in the behavior of air 
from the ideal gas relations. At high altitudes the atmosphere 
is so rarefied that the equations of continuum mechanics no 
longer hold, while the temperature produced by hypersonic 
flight through the lower atmosphere is sufficient to dissociate 
and ionize air-molecules. This-article deseribes recent prog- 
ress in our understanding of these effects. The bibliography 
is not intended to be at all complete, but is selected to serve 
as a guide to literature on the subject. 

The nearly exponential decrease of atmospheric density with 
altitude results in a rapidly increasing molecular mean free 
path so that, at about 200,000 ft, the average distance be- 
tween molecular collisions becomes comparable with the di- 
ameter of a missile. Under such conditions new aerodynamic 
phenomena arise. “Slip-flow’’ is used to designate this regime 
since the average velocity of the molecules near a surface 
has a finite value and the air appears to be slipping over the 
body. At higher altitudes free molecule flow occurs in which 
molecules striking the missile are not deflected back to it by 
further intermolecular encounters. 

The air densities, and hence the altitudes, at which these 
various flow regimes exist are somewhat difficult to specify. 
Knudsen’s number K, the ratio of the molecular mean free 
path to some pertinent body dimension, is usually used as a 
criterion. From the kinetic theory of gases it is easy to show 
that K is approximately equal to the ratio of Mach number 
to Reynolds number. Classically, the limits of slip flow 
have been chosen as 0.1 < K < 1.0 (see, for example, Ester- 
mann (5)!) but some care must be used in selecting an appro- 
priate “typical’? body dimension in any given problem. 
Tsien (24) discusses the problem in some detail, and more re- 
cently the group working on low density flows at Berkeley has 
reported experimental results on the onset of slip flow over 
flat plates and around spheres (see, for example, Laurmann 
(11)). In an excellent review of recent developments in hy- 


Received Jan. 24, 1958. 
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personic flow in this journal Lees (13) points out that the 
local air density in the boundary layer near a stagnation point 
is so much higher than the ambient density that any criterion 
for the onset of slip flow based only on free stream conditions 
is of doubtful validity. The slip flow regime is postponed to 
much lower densities than has previously been supposed. A 
theoretical problem of some interest arises in trying to fit the 
boundary conditions at a surface into either a kinetic theory 
of the flow or into the Navier-Stokes equations. Several 
articles on this problem have appeared, for example, Wang 
Chang and Uhlenbeck (25) and Yang and Lees (26). Much 
of the work has dealt with concentric rotating cylinders with 
comparatively small differences in speed and temperature. 
Extrapolation to missile conditions is clearly risky until the 
theory is more firmly established than at present. 

The gas dynamic problem in free molecule flow becomes 
one of finding out how various atoms and ions reflect from 
solid surfaces. In its details the process must be a very 
complex one involving the binding force between particle and 
surface atoms as well as the geometry of scattering from a 
microscopically rough surface. In practice, the custom has 
been to assume that some fraction of the particles reflect 
specularly and that the rest are momentarily absorbed and 
then emitted isotropically. The mean energy of the emitted 
molecules is specified by a parameter called the ‘“‘accommoda- 
tion coefficient’? which measures the extent to which these 
molecules reach thermal equilibrium with the wall. Much 
new data is becoming available on various gas and solid com- 
binations (Merlic (17)) but progress on fundamentals seems 
to be temporarily at a standstill. This is unfortunate if for 
no other reason than that a similar concept of partial specular 
and diffuse reflection is being used quite widely in studying 
the intermediate regime mentioned above. 

Since a ballistic missile spends most of its flight time in the 
outer atmosphere one might at first suppose that even the 
small amount of air encountered there would have some ef- 
fect on the missile’s trajectory. This is not the case, how- 
ever, and the missile actually does ‘re-enter’ the lower at- 
mosphere before aerodynamic and thermal problems become 
significant. Satellites, on the other hand, because of their 
prolonged travel through the outer atmosphere, are gradually 
affected by the rarefied air and eventually fall back to earth. 
Thus the problems of free molecule and slip flow are asso- 
ciated primarily with satellites, while high temperature gas 
effects are predominantly the concern of ballistic missile 
designers. 

Aerodynamic and heating problems arising in hypersonic 
flight through the atmosphere are complicated by a number of 
molecular and atomic processes not encountered previously. 
Among these are vibrational and electronic excitation, disso- 
ciation, chemical reaction and ionization. New mechanisms 
for heat transfer are provided by radiation and by atomic re- 
combination on the surface of a missile. Other problems 
arise from the generation of an ionized wake. 

The equilibrium properties of any gas can be computed 
from statistical mechanics once all of the quantum mechanical 
energy levels are known. Fortunately nearly all of the neces- 
sary data have been accumulated in three decades of physical 
research. Tables giving thermodynamic properties of many 
pure gases have been prepared by the National Bureau of 
Standards (19) and of air by Feldman (6), Logan (15), and 
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Hochstim and Arave (8). Normal and oblique shock rela- 
tions are also available over a wide range of conditions. 

One notable exception to the fund of physical data men- 
tioned above is the dissociation energy of nitrogen. Spectro- 
scopic evidence indicated a value of either 7.38 or 9.76 ev. 
An unambiguous choice was made possible by a purely gas 
dynamic measurement. The velocity of a plane shock wave 
depends upon the state of the gas behind it so that the degree 
of dissociation, which in turn depends upon the dissociation 
energy, can be determined from a velocity measurement. 
Toennies and Greene (22) measured the velocity of incident 
and reflected shocks in a shock tube and found that their 
data led to the higher value of 9.76 ev. 

In discussing the properties of a gas, it is customary to refer 
to the ratio pV /RT as the “compressibility factor.’’ Two in- 
dependent effects cause departures from the ideal gas value 
of 1: intermolecular forces leading to the van der Waals 
coefficients, and a change in the number of particles by disso- 
ciation and ionization. The first effect is quite negligible at 
the densities encountered in flight, although it may make a 
detectable contribution in hypersonic wind tunnel work. 
The second effect is important in flight and in shock tube work, 
while it is absent in hypersonic tunnel experiments because of 
the low operating temperature. 

Because the excitation of the internal degrees of freedom 
of oxygen and nitrogen, as well as chemical reactions between 
them, proceeds at rates which are comparable to the rates of 
change in the flow about a hypersonic missile, knowledge of 
these rates is required for a complete solution of the gas 
dynamic problem. 

Vibrational excitation proceeds according to a linear re- 
laxation equation; that is, the rate of change of vibrational 
energy is proportional to the amount that vibrational energy 
is out of equilibrium with the mean translational tempera- 
ture. Vibrational relaxation times have been measured in 
the shock tube by Blackman (1) for pure oxygen and nitrogen, 
and the effectiveness of oxygen-nitrogen collisions in exciting 
oxygen has been determined. At 3000 K and 1 atm the 
pure gas relaxation times are 2 usec and 30 usec, respectively. 
A collision between nitrogen and oxygen molecules is about 
40 per cent as effective in exciting vibration in oxygen as an 
oxygen-oxygen collision. The effectiveness of oxygen in 
exciting nitrogen has not been measured but may be expected 
to have a lower value as well. Vibrational relaxation requires 
many collisions, not so much because the quantum level 
spacing is large compared to kT’, but because the probability 
of excitation in a collision which is sufficiently energetic is 
extremely small. Since the total energy per molecule which 
can be involved in vibration is.only kT at best, vibrational 
relaxation is not particularly important to hypersonic flow in 
itself. As a stage in the mechanism of dissociation however, 
it plays an important role. 

Lighthill (14) points out that dissociation is complete in 
any gas at a temperature far below D/k (where D is the dis- 
sociation energy) because of the very large volume of phase 
space available to two atoms as compared with one molecule. 
Thus oxygen and nitrogen are completely dissociated at 6000 K 
and 10,000 K, respectively, even though the values of D/k 
are 59,400 K and 113,300 K. Dissociation proceeds rela- 
tively slowly because only a few molecules in the tail of the 
Maxwellian distribution function have enough energy to break 
the molecular bond. 

The rate of dissociation is given by the fraction of collisions 
which are sufficiently energetic, (D/kT)*e~?*", multiplied 
by an efficiency factor which is less than unity. The param- 
eter s must also be determined by experiment. Measure- 
ments by Matthews (16), Eckermann (4), Camac (3), and 
Treanor (23) indicate that the efficiency factor is between 
3 and 3; for oxygen and nitrogen and that s has a value be- 
tween 2 and 3. Unfortunately none of the data thus far 
available is accurate enough to establish these values with 
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precision. The rates are highly dependent on temperature, 
having both a power of 7’ and the Boltzmann factor in the 
expression. Thus no part of the process is correctly describ- 
able by a dissociation relaxation time, contrary to statements 
by some authors. It is preferable to state either the instan- 
taneous rate itself or the time to reach, say, 90 per cent of the 
equilibrium value. 

Because the binding energies of oxygen and nitrogen differ 
by nearly a factor of two, it is clear that both the dissociation 
rates and equilibrium degree of dissociation are quite different 
at any given temperature. In air this means that the order of 
equilibration is: oxygen vibration, oxygen dissociation or 
nitrogen vibration, nitrogen dissociation. The cross-over of 
rates just mentioned may account for an apparent over- 
dissociation of oxygen behind a plane shock wave in air that 
was observed by Resler (21). More experimental work will 
be needed before a clear picture of these fundamental proc- 
esses is available. 

Although the very large energies absorbed in dissociation 
and the change in number of particles may be expected to 
have a strong influence on hypersonic flow patterns, as shown 
by Lighthill (14) and Eckermann (4), boundary layer heat 
transfer is not so strongly effected. Lees (12) investigated the 
role of atomic recombination in heat transfer through a hyper- 
sonic boundary layer and found that the rate depends mainly 
on the difference between the total stream enthalpy and the 
surface enthalpy. The rate for a boundary layer in complete 
thermodynamie equilibrium differs by only about 30 per cent 
from that for the case where atoms recombine only on the 
surface. The diffusion of atoms is limited by collisions with 
molecules to such an extent that little increase in heating 
occurs even though the energy of recombination may be 
delivered right at the surface. 

Compounds of nitrogen and oxygen, principally NO, NO» 
and NO*, appear in equilibrium in trace amounts at elevated 
temperatures. These are probably not important in them- 
selves but do serve as catalysts for the very important pro- 
duction of electrons and ions. The transient existence of 
similar compounds might also be of significance but at present 
too little is known to offer much concrete evidence. The reac- 
tion with the lowest energy barrier for electron production is 
N+0O0—NOt++e. So little is known about the cross sec- 
tion, however, that one cannot be sure that this is actually 
an important source of electrons. Bond (2) and Petsehek 
(20) have studied the ionization of argon behind strong shock 
waves and conclude that the rate of approach to equilibrium 
ionization is controlled by the electron population. If the 
same is true for air, then knowledge of the first source of elec- 
trons is extremely important. 

The energy radiated from an electronically excited and 
ionized gas can be computed from the particle wave functions 
and the Kinstein radiation coefficients. Such computations 
are often very laborious with the result that few have been 
carried out. Failing a purely theoretical attack it is possible 
to make some intelligent guesses on the strength of a whole 
electronic band from a knowledge of other band systems. 
Hindmarsh (7) discusses this problem in some detail. Mever- 
ott (18) makes use of the recent experimental data of Keck 
et al. (9) and Kivel et al. (10) on the emissivity of hot air to 
to show that the radiant heat transfer to a ballistic missile is 
only about 10 per cent of the convective heat transfer. Thus 
the radiation problem seems to be of primary importance 
only to recovery of a satellite vehicle, for example a translunar 
rocket. For an excellent summary of the radiant heat trans- 
fer problem the reader is referred to the paper by Meyerott 
(18). 
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Some Effects of Oxidizer Concentration and Particle Size 
on Resonance Burning of Composite Solid Propellants 


LEON GREEN Jr.’ 


Aerojet-General Corp., Azusa, Calif. 


A comparison of four experimental, composite propel- 
lant formulations with nominally constant burning rates, 
tested at 170 F in internal-burning tubular charges par- 
tially supported the widely held belief that the tendency 
toward irregular reaction increases as the energy content 
of the propellant is raised. However, the highest energy 
formulation showed a smooth reaction, apparently as a re- 
sult of a burning rate lower than anticipated on the basis 
of Crawford bomb measurements. The results suggest 
that the rate of energy release (the product of burning rate, 
propellant density and heat of explosion) affords a better 
measure of the relative stability behavior of various formu- 
lations than does either burning rate or energy content 
alone. A subsequent brief comparison at 170 F of three 
other composite propellants with constant energy content, 
but with varying burning rates (obtained by varying the 
oxidizer particle-size distribution), indicated that the tend- 
ency toward irregular reaction increased as the burning 
rate increased or as the particle size became finer. Owing 
to propellant processing difficulties, efforts to study the in- 
fluence of burning rate alone while controlling both oxi- 
dizer concentration and particle-size distribution were not 
successful. It is recommended that the apparent in- 
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fluence of a coarse oxidizer particle-size distribution in 
suppressing resonance burning be investigated further. 


Introduction 


HE tendency of solid propellant charges to exhibit res- 

onance burning is obviously a function of propellant 
composition, since some propellants never show this behav- 
ior, while others almost always do, at least in tests of in- 
ternal-burning charges conducted at extreme conditioning 
temperature. On the basis of development experience, a 
qualitative proposition has been advanced; namely, resonant 
burning is encountered only in motors employing high energy 
propellants, i.e., propellants with a high heat of explosion or a 
high specific impulse. Examples of such propellants are the 
double-base compositions using nitrocellulose of high nitration 
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Fig. 1 Particle-size distribution of ammonium perchlorate 
employed in propellants A, B, C and D 


(such as JPN ballistite) and composite propellants employing 
ammonium perchlorate as oxidizer. The lower energy pro- 
pellants that have not reportedly manifested unstable com- 
bustion include the slow-burning, double-base formulations, 
such as Russian cordite, and the German and Japanese pro- 
pellants described by Wimpress (1),* as well as the composite 
propellants using potassium perchlorate, ammonium picrate, 
sodium nitrate, potassium nitrate and ammonium nitrate. 
A general description of these various types of propellant is 
given in (2). Wimpress states that the trend toward resonance 
increases with increasing burning rate, but this observation 
does not apply to the fast-burning potassium perchlorate 
formulations. Although the correlation with energy content 
has been observed in homogeneous double-base propellants 
as well as in composite ammonium perchlorate and organic 
fuel systems, the burning rates of both propellants systems also 
generally increase with energy content, and it was not possible 
to state whether the observed trend in reaction stability was 
attributable to either energy content or burning rate alone. 
In the composite systems, moreover, the burning rate in- 
creases with decreasing oxidizer particle size. The objective 
of the present investigation was to distinguish the effects of 
these variables. 


Experimental Propellants 


It was desired to study the effect of (i) varying energy con- 
tent (determined by oxidizer concentration) in propellants of 
approximately equal burning rates but with particle-size dis- 
tribution not controlled, (ii) varying burning rate in propel- 
lants of equal energy content and constant particle-size dis- 
tribution and (iii) varying particle-size distribution in propel- 
lants of equal energy content but with burning rate not con- 
trolled. This choice of experiments was determined by the 
fact that the effectiveness of burning-rate accelerators (added 


’ Numbers in parentheses indicate References at end of paper. 
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Propellant D 


Propellant Cc 


Reference Length = 0.1mm.——a k— 


Fig. 2 Photomicrographs of samples of propellants A, B, C and 
D (taken by transmitted light) 


in amount of a few weight per cent or less) in compensating 
for the effect of oxidizer concentration and particle size is 
limited; i.e., it was not considered possible to control both 
burning rate and particle size while varying the oxidizer con- 
centration significantly, or to control both burning rate and 
concentration while varying the particle size. 

The first objective was achieved with a moderate amount of 
laboratory development. Four experimental propellants, 
identified as A, B, C and D, were formulated using different 
concentrations and oxidizer grinds. Particle-size distributions 
of the ammonium perchlorate grinds employed and _photo- 
micrographs of propellant samples are presented in Figs. 1 
and 2, respectively. The comparison of interior ballistic 
properties presented in Table 1 shows that the strand burning 
rates at standard conditions (60 F and 1000 psia) in a Craw- 
ford bomb (3) were approximately equal. Subsequently, how- 
ever, a test temperature of 170 F was selected (as discussed 
later), and at this temperature a greater divergence of burn- 
ing rates was manifested, with propellant D exhibiting a 
markedly lower temperature sensitivity than did the other 
formulations. This divergence was accepted, since the de- 
velopment of another set of propellants would have required 
an effort exceeding the authorized scope of the program. 

The second objective was not achieved. Three experimental 
propellants (designated E, F and G) were formulated using a 
conventional (medium) oxidizer grind, but having the burn- 
ing-rate accelerators present in zero, medium and relatively 
high concentrations. Unfortunately, it was found that these 
propellants would cure satisfactorily only when the burning- 
rate accelerator (which also apparently was required to ac- 
celerate the curing process) was present in the relatively high 
concentration. After considerable effort to obtain satisfactory 
curing properties proved unsuccessful, the development of 
propellants E, F and G was discontinued. 

The third objective mentioned above was accomplished by 
the formulation of propellants H, I and J with constant 
(fairly high) concentration of burning-rate accelerators, but 
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Table 1 Comparative ballistic properties of experimental 
composite propellants with varying energy content but 
approximately constant burning rate under standard 


conditions 
| Propellant A B C D 
Total oxidizer con- 
tent, per cent 70 75 80 85 


| 

| Density, lbm/in.* 

| Strand burning rate 
| rate at 1000 psi and 
| 

| 


60 F, in./sec. 0.35 0.37 0.34 0.34 
Heat of explosion,! 
cal/g 698 864 1010 1050 


Heat of explosion/- 

unit vol., Btu/in.? 72.9 92.5 112 
Nominal rate of heat 

release at 60 F and 

1000 psi, Btu/in.? 

sec 25.5 34.2 38.1 41.1 
Flame temperature, 

°R 3570 4320 5030 5330 
Specifie impulse,? Ibf- 

sec/lbm 211 226 239 251 
Impulse density, 

sec /in.* 12.8 13.5 14.7 16.0 
Average specific heat 

ratio, y 1.22 1.24 1.22 
| Molecular weight of 
gases in combustion 
chamber 21.9 23.6 25.5 27.2 

Characteristic gas ve- 

locity, c, fps 4350 4600 4800 4860 


1 Enthalpy change between propellant at 25 C and reaction 
products (in relative concentration corresponding to equi- 
librium composition of flame temperature) at 25 C, with water 
and hydrogen chloride in liquid state. 

2 Assuming mobile equilibrium of products during expansion 
from 1000 psi. 


with the three different oxidizer particle-size distributions 
shown in Fig. 3. Photomicrographs of propellant samples are 
shown in Fig. 4, and the interior-ballistic properties of these 
propellants are given in Table 2. It may be noted that even 
the coarsest-grained propellant showed a higher burning 
rate than any of those described in Table 1. 


Test Procedure 


Internal-burning, tubular charges of each of the experi- 
mental propellants were statically tested at 170 F and at 
various nominal pressure levels, as determined by the ratio, K, 
of burning surface area to nozzle throat area. The test tem- 
perature was selected on the basis of preliminary firings which 
indicated that irregular reaction of propellants of the general 
class involved was most likely to be exhibited at high tem- 
peratures. The internal charge diameter was 3.0 in., the web 
thickness 0.5 in. and the length 16.5 in. This short grain 
length was selected in order to minimize any effects due to 
high aft-end gas velocity (4). 

In the tests of propellants A, B, C and D, the conventional 
pressure measurements were supplemented by high frequency 
response instrumentation, in order to check the possible 
modes of acoustical oscillation associated with any rough 
burning. The detailed results of these measurements are 
available elsewhere (5) and will not be discussed here. In 
summary, it will suffice to say that oscillations of frequencies 
corresponding to both longitudinal and transverse modes were 
detected, but the amplitudes generally did not exceed a few 
psi, even in the tests exhibiting irregular low frequency pres- 
sure vs. time histories. This behavior is consistent with the 
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Fig. 3  Particle-size distribution of ammonium perchlorate 
employed in propellants H, I and J 


Propellant H 


Propellant TIT Propellant J 


Fig. 4 Photomicrographs of samples of propellants H, I and J 
(taken by reflected light) 


Table 2 Comparative ballistic properties of experimental 
composite propellants with varying burning rate but with 
constant energy content 


Propellant H I J 
Total oxidizer content, not in- 
eluding burning-rate accelera- 
tors, per cent 76 76 
Density, lbm/in.* 
Strand burning rate at 1000 


psia and 60 F, in./sec 0.70 0.48 1.24 
Heat of explosion, cal/gm 932 932 932 
Heat of explosion/unit volume, 

Btu/in.’ 104 101 105 


Nominal rate of heat release at 
60 F and 1000 psia, Btu/in.? 


sec 73.0 48.3 130 
Flame temperature, °R 4830 4830 4830 
Specific impulse, lbf sec /Ibm 233 233 233 
Average specific heat ratio 1.20 1.20 1.20 
Molecular weight of gases in 

combustion chamber 26.2 26.2 26.2 
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observation that the amplitude of acoustical oscillation in the 
cavity (“sonant” burning, the existence of which is a neces- 
sary but not sufficient condition for the occurrence of gross 
irregularities in the mean pressure-time history, or ‘“‘resonant”’ 
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Fig. 5 Pressure-vs.-time curves obtained in the static tests of 
propellant A 
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Fig. 6 Pressure-vs.-time curves obtained in the static tests of 
propellant B 
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Fig. 7 Pressure-vs.-time curves obtained in the static tests of 


propellant C 


burning) does not necessarily determine the magnitude of the 
gross irregularities (6, 7). In the tests of propellants H, I and 
J, no high frequency measurements were made. 


Experimental Results 


The pressure-time curves obtained from the tests of propel- 
lants A, B, C and D are presented in Figs. 5, 6, 7 and 8, re- 
spectively, and the interior ballistic data are summarized in 
Table 3. These results generally confirm the belief stated 
earlier that irregular reaction tends to be more pronounced in 
high energy propellants than in low energy ones. The firing 
curves obtained with propellant A (which show a sensitivity 
of pressure to K ratio so high that the measured K values 
were not sufficiently accurate) are uniformly smooth, while 
those obtained with propellants B and C show an increasing 
degree of irregularity manifesting resonant burning. The ex- 
tended “‘tail-off’’ observed in the tests of propellants A and |3 
is attributable to a nonuniform distribution of oxidizer along 
the length of the grain, because of the settling of the particles 
when cast in low concentration in a fuel of low viscosity; the 
high concentration end of the grain burns out before the low 
concentration end does. This effect becomes less pronounce:| 
as the nominal concentration level is increased, but some tail- 
off is always observed with grains of constant web thickness 
because of erosive-burning effects. The presence of tail-oif 
in many of the tests required a compensation of this effect in 
reduction of the firing curve data. The mean pressures and 
mean burning rates reported herein were computed by extra- 
polating the principal portion of the pressure-time decline to 
the zero-pressure ordinate. The point of intersection was used 
to define the effective duration of reaction at the nominal 
mean pressure level, thus disregarding the effects of tail-off. 

In the case of propellant C, pronounced “secondary pres- 
sure peaks” (as distinguished from the initial “ignition peak’’) 
occur early in the tests corresponding to curves 2 and 4. The 
secondary pressure peaks observed with propellants B and C 
must be distinguished from the smooth, broad “hump” which 
is observed at the midpoint of the pressure-time histories 
obtained with propellant A. This burning anomaly is a 
property of the particular class of composite propellants em- 
ployed. These often show a nonuniform distribution of burn- 
ing rate across the grain web, with the maximum value 
occurring at the center of the web. This effect is most pro- 
nounced in the formulations employing relatively low oxidizer 
concentrations and becomes negligible in formulations with 
high oxidizer concentration. In composite propellants, this 
effect may be attributable to a nonuniform distribution of 
oxidizer, but a similar burning-rate nonuniformity has also 
been observed in extruded, double-base propellants (1). 
Surprisingly, the highest energy charges (propellant D) 
showed uniformly smooth curves. This result is attributed to 
two factors: First, the coarse particle size employed in this 


CURVE NO RATIO (Avg) | 
| 
4 326 | 
308 
5 200 | 
6 201 
1400 
y ‘200 | 
2 
w 
800 
3 400 
° 2 a 6 8 1o 12 14 16 18 20 
TIME (sec) 
Fig. 8 Pressure-vs.-time curves obtained in the static tests of 
propellant D 
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re 3 Summary of data obtained from the static testing 
of internal-burning, tubular charges of experimental 
propellants A, B, C and D 


Nominal 
heat 
Nozzle Avg. release 
throat Avg. burning rate, 
Curve Avg. area, pressure, rate, Btu/ 
No. K in.? psia in./sec —_in.?/sec 


PROPELLANT A (Fig. 5) 
1 355 0.482 1154 0.412 30.1 
2 353 0.483 1152 0.411 30.0 
3 354 0.482 1030 0.382 27.9 


PRoPELLANT B (Fig. 6) 


1 378 0.450 1900 0.662 61.2 
2 328 0.522 804 0.324 30.0 
3 236 0.728 663 0.299 27.7 
PROPELLANT C (Fig. 7) 
1 280 0.617 1939 0.643 72.0 
2 281 0.615 1439 0.534 59.8 
3 242 0.703 966 0.425 47.6 
4 246 0.703 1054 0.494 55.3 
5 241 0.694 953 0.460 51.5 
6 247 0.694 969 0.460 51.5 
7 225 0.770 898 0.470 52.6 


PropELLANT D (Fig. 8) 


1 355 0.484 1105 0.331 40.0 

326 0.523 996 0.319 38.5 
3 326 0.524 992 0.314 30.9 
4 308 0.557 824 0.268 34.5 
5 200 0.850 425 0.223 26.9 
6 201 0.848 424 0.220 26.6 


formulation apparently tends to promote stable combustion 
(as has been suggested by development experience) and thus 
may counteract the effect of the increased energy content. 
Second, the burning rates of the propellant D charges were 
found to be considerably lower than the values indicated by 
strand burning measurements using the Crawford bomb, an 
effect which is not readily explainable. Accordingly, the rates 
of energy release (here taken as the product of heat of ex- 
plosion, the propellant mass density and the linear burning 
rate) experimentally observed with the propellant D charges 
are actually lower than the corresponding rates observed with 
propellant C, as is shown in Table 3. These results suggest 
that the rate of energy release may possibly afford a measure 
of the relative reaction stability. 

The pressure-vs.-time curves obtained in the static tests of 
propellants H, I and J are grouped in Figs. 9, 10 and 11, re- 
spectively. The results shown herein followed the anticipated 
trend; the finest grained, fastest burning propellant, J, 
showed the most irregular behavior. In general, however, the 
irregularities were observed in the tests at low pressure levels, 
with the charges fired at low AK ratios. The inconsistent 
relation between average pressure and K ratio for curves 4 
and 5 of Fig. 11 is attributed to grain fracture during reso- 
nance. Itisnot known whether the nozzle separation indicated 
in Fig. 11 was attributable to a high initial resonant-burning 
peak or to the presence of undetected, internal fissures within 
the propellant. 

The charges of propellant H, with the intermediate dis- 
tribution of oxidizer particle size, showed a smoother reaction, 
with only mild irregularities being observed. The extended 
regression (tail-off) observed in curve 5 of Fig. 9 indicates 
that some fracturing or fissuring of the grain occurred during 
the run. The tests of the slower burning propellant I showed 
the least tendency toward resonance, but the results were 
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Fig. 9 Pressure-vs.-time curves obtained in the static tests of 

propellant H 
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Fig. 10 Pressure-vs.-time curves obtained in the static tests of 
propellant I 
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Fig. 11 Pressure-vs.-time curves obtained in the static tests of 
propellant J 
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partially obscured by the difficulty in fabricating satisfactory 
charges of this propellant, owing to the poor castability which 
is characteristic of a very coarsely ground oxidizer fraction. 

Of the five tests performed with this propellant, three yielded 
pressure-time histories which, although smooth, showed initial 
pressures excessive enough to blow the gaskets sealing the two 
plugged fore-end bosses for the unused high frequency re- 
sponse pressure pickups, and thus exhibited an extremely re- 
gressive character. For comparison, the results of one such 
test are included as curve 1 of Fig. 10. Also, curve 2 of this 
figure represents a test in which the pressure was just high 
enough to cause a leak at one fore-end pickup location, thus 
producing a moderately regressive pressure-time characteris- 
tic. Only one test of propellant I (curve 3 of Fig. 10) was com- 
pletely satisfactory. Nevertheless, it is emphasized that the 
erratic pressure-vs.-time curves obtained with propellant I 
are not characteristic of resonant burning, but are attributed 
to poor propellant consolidation. The low pressure-time 
integral shown by curve 1, for instance, manifests a significant 
loss of unburned propellant through the nozzle. 


Conclusion 


The study described above consituted an effort to investi- 
gate systematically certain trends in relative reaction stability 
with propellant composition which had been suggested by the 
experience accumulated in various rocket development pro- 
grams. The study was hampered by a limited scope and by 
difficulties in processing experimental propellants of widely 
varying compositions which did not conform to conventional 
manufacturing practice. Nevertheless, the results permit 
the conclusion that the tendency of a propellant to exhibit 
resonance burning may not be judged on the basis of either 
energy content or burning rate alone, but that a nominal rate 
of energy release involving the product of the two quantities 
affords a better measure for estimating relative reaction sta- 
bility. The apparent effectiveness of a coarsely ground oxi- 


dizer fraction in suppressing resonance was not investigated 
thoroughly enough to permit any judgment as to whether this 
beneficial effect is due to the coarse particle size per se, or to 
its influence in lowering the burning rate. Further efforts to 
distinguish particle-size and burning-rate effects are recom- 
mended. 
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A Generalized Optimization Procedure for 
N-Staged Missiles 


LEON WEISBORD' 


Wright Aeronautical Division, Curtiss-Wright Corp., Wood-Ridge, N. J. 


This report is an extension of the Goldsmith report (Ref- 
erence 2) and outlines a simple method for determining 
the optimum weight distribution for a staged rocket with 
different structural factors and propellant specific impulses 
in each stage. It is applicable for a rocket of any number 
of stages, and is not limited to either a horizontal or verti- 
cal trajectory. The report can be used in a first design 
estimate to determine the weight distribution among the 
various stages and the weight breakdown within each 
stage. A sample problem is included to illustrate the 
method of application. 


Nomenclature 
C = propellant exhaust velocity 
t 
G = V+ sin oat 
= structural factor 
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M = payload 

M, = gross weight 

Mm = miscellaneous weight (structures, engines, controls, etc.) 
M, = burnout weight 

M, = propellant weight 

= thrust 

tp = burning time 

= velocity 

0 = angle between velocity vector and horizontal 
Subscripts 


1, 2, 3 = number of stage 
N = last stage 
N — 1 = next to last stage, ete. 


Introduction 


T IS necessary for the rocket designer to know the best 
method of apportioning the weight of a staged rocket. It 
is possible to distribute the weight in an infinite variety of 
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ways to achieve the same final performance; however, there 
is only one minimum weight distribution. Knowledge of this 
theoretical minimum weight distribution is obviously an im- 
portant bit of information in the preliminary design of the 
staged rocket. 

The problem has been treated by Malina and Summerfield 
(1).2. Their analysis held for a rocket of an, arbitrary number 
of stages which used the same impulse and structural factor in 
each stage. Under these conditions the following relation 
holds for minimum take-off weight 

M Mey Moy-1 Mz, 


However, the problem of optimum staging when structural 
factors and propellant specific impulses differed remained un- 
solved. Goldsmith (2) realized this lack and devised a method 
for the optimization of two-stage rockets under conditions of 
different impulses and structural factors. This pointed up 
the need for a general method of optimization for rockets of 
any number of stages. 

Such a method is presented in this report. It is valuable as 
a first design weight breakdown and is quick and easy to 
apply 


Analysis 


The rocket velocity is given by Equation [2] which is the 
relation for drag-free flight with constant exhaust velocities 


t M 
G=V+ sin = C, In Mo, + C2In Me, 
Moy 
Cy In hd [2] 


It can be seen that G in Equation [2] is equal to the velocity 
that would be attained if the rocket were fired horizontally. 
If the rocket were fired vertically, tw g sin @dt = gt,. 
An evaluation of the term tg g sin 6dt is difficult since it 
depends on the particular trajectory and the engines used. 
In using this optimization procedure, i” g sin @dt must be 
approximated, but it will be between 0 and gt, depending upon 
the trajectory followed. 

Each stage of the rocket consists of propellant, payload 
and miscellaneous additional hardware consisting of tanks, 
structure, guidance, engines, etc. The miscellaneous weight 
will be lumped together and included in the calculation by 
defining a structural factor K as 


miscellaneous weight = K(M, — payload)....... [3] 


This analysis will first derive expressions for a three-stage 
rocket, and follow with the general case. 

The empty weight M is defined as the payload plus the 
miscellaneous weight. 

Expressions will now be written for the burnout weight at 
the end of each stage. 


Mo, = M+ KxM,, — M) = K:M,, + (1 — Ks) ) 
Me, = KiM,, — 
Mo, = KiM,, + (1 — Ki)M,. } 


Equation [2] can be rewritten as 


M,, My 
C1 


Cy In 


Mo, 
K3M,, + (1 — K3)M 


? Numbers in parentheses indicate References at end of paper. 
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or 
C, In My, — Ci In [KiM,, + (1 — K1)M,,] + C2 In Ma, 
— C2 In [K2M,, + (1 — K2)M,,] + In M,, — 
C; In [Ks + (1 — K3)M] —G =0......... [5a] 


It is desired to minimize M,, for a given velocity. There- 
fore, in Equation [5 or 5a], M,, and M,, become the variables 
and the following must be satisfied for minimum take-off 
weight 


My, _ 
oM,, 


It is extremely difficult to write [5] in the form M,, = 
f(M,,, M,,) so that the partial derivatives can be found ex- 
plicitly. However, it will be recognized that [5a] is of the 
form 


where G is a constant for any particular mission. 

This function can be differentiated implicitly quite easily. 
The necessary partials (Equation [6]) can be obtained from 
the following relations holding for Equation [7] 


af af 
Oz Or oy 
ox of oy of 
oz Oz 
In our problem, Equation [5a] can be expressed as 
[9] 
and we will find 
af 
mM, a 2s 
aM, aM, 
The above partials will be set equal to zero or 
M ra) oM of 
— = 0 = —.. 
aM, 
The differentiation gives 
of Ki) 
Man KiM,, + (1 — Ki)Mn | 
M,, KeM,, + (1 — Kz)M,, | 
[12] 
Of Cf{1 — Ke) 4 
oM,, K2M,, + (1 — | 
C; 


| 
M,, K:M,,+(1 — K3)M} 


Setting Equations [12] equal to zero as indicated by Equa- 
tions [11] for minimization, obtain after algebraic manipula- 
tion 


— — Ke)(C2 — C1) + — 
Ma 
and 
* (1 — Ks)(1 — Ks)(C3 — C2) + — Ks) 
93 


igated 
er this poe 
rts to 
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1 
0 6] 
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te.) 
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Rewriting Equation [5] as 


1 


a. + C2 In + : 
M,, My, 
Ki+(1- K2+(1 


G =C,In 


C3 In 


M 


93 


we can now determine the necessary take-off weight and in- 
dividual stage weights. 

For any given problem it is necessary to establish C, C2, Cs, 
Ky, Ke, Kzand G. The steps in the solution are: 


(a) Solve [12b] for 
My, 
M,, 


Mo, 
My 

(c) A short trial and error procedure follows. Select a 
numerical value for M,,/M,, and obtain numerical answers 
from Equations [12a and 12b]. Using the values found for 
M,./M,,, M,,/M,, and M/M,,, solve [5b] for G. If the value 
of G thus found is not the value desired, repeat the process 
using a different value of M,,/M,,. 

The method outlined is, in actual practice, quite fast. The 
trial and error process takes five to ten minutes to reach a 
solution. 

(d) Finally, obtain a maximum value of 

M Mya M 


My, Mo, Mo, 


(b) Solve [12a] for 


This value of //M,, is the optimum theoretical value for 
the conditions prescribed. 
It will be immediately noted that when 


C,=C.=C; K, = Ke = K; 


the Malina-Summerfield conditions are encountered. Under 
these conditions Equations [12a and 12b] break down and 
become simply 


which is recognized as Equation [1], the Malina-Summerfield 
optimization criterion. 

The analysis can now be easily extended to cover any num- 
ber of stages.* If N represents the number of stages, N being 
the last stage to fire, Equations [12a and 12b] can be extended 
to give the N — 1 relations necessary for an N stage problem 
solution. That is 


M — Ky-1) 


= Ky-1)(1 — Kw)(Cw — Cy-1) + CyKw-1(1 — Ky) 


Moy Cy-2Ky-1(1 — Ky-2) 


Mon-1 


T 


STAGE Missile | MISSILE 2 MISSILE 3 
No. PROP c PROP c PROP c 
! LOx-uP 8,950 Fo-NHs 6,950 O2-H2 u,780 
2 LOX-JP 10,000 11,200 O2-H2 12,380 
| 3 LOX-JP 10,000 Fa-NH3 11,200 O2-H2 12,380 


MISSILE | 
50 


r 


40| | 

30| 

MISSILE 3 | 
| 

07 08 09 410 12 


K, 


Fig. 1 Optimized three-stage missiles for the case K. = K, + 
0.01, K; = K, + 0.02, G = 30,000 fps 


group of widely differing structural factors and exhaust veloci- 
ties will be selected. Let us use 


C, = 9,000 K, = 0.2 G = 30,000 
C2 = 15,000 Ky = 0.25 
C3 = 25,000 K; = 0.40 


Substitute the appropriate numerical values into Equa- 
tions [12a and 12b]. Then, perform the trial and error process 
mentioned before by first selecting a value of M,,/M,, and 
then solving Equation [5b] for G using the values determined 
from Equations [12a and 12b]. 

This procedure yields M,,/M,, = 0.678, M,,/M,, = 0.260 
and M/M,, = 0.238. Therefore, M,,/M = 23.8, or for every 
pound of payload, 23.8 lb of take-off weight are needed. In 
addition, the size and weight of propellant for each stage can 
be determined easily by using the structural factors as indi- 
cated in [4]. Table 1 gives the complete weight breakdown 
for a payload of 1000 lb and approximate engine sizes and 
burning times. The engine thrusts were arbitrarily made 
equal to twice the gross weight. 

Actually, it was never shown whether the procedure leads to 

a minimum or a maximum gross 
weight. It is easiest to decide this 
question by checking any other com- 
bination of stage weights. Since the 
.. [13] procedure outlined describes either 
a maximum or a minimum, if any 


Moy-1 


etc. The Nth equation necessary for solution is, of course, 
Equation [2], written in the form of Equation [5b]. 
Sample Calculation and Check for Minimization 


In order to illustrate the use of the method, a sample prob- 
lem will be performed. This sample problem is not intended 
to illustrate any particular missile application. Rather, a 
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(1 — Ky-2)(1 — Ky-1)(Cw-1 — Cn-2) + Cy-1Ky-2(1 — Ky-1) 


Mon-2 


other combination of weights gives 
higher ratio of gross weight to pay- 
load than 23.8, then the procedure is 
one of minimization of the gross take-off weight as desired. 
Selecting M,,/M,, = 0.600 and M,,/M,, = 0.400, obtain 
by solution of Equation [5b] 17/M,, = 0.1585 Then M,,/M 


5 Vertregt (3) analyzed the problem through the use of La- 
rangian multipliers. His analysis presents an implicit solution 
or the general case and a formal solution for the case of equal 

exhaust velocities. 
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= 26.3, which shows that this procedure is indeed one of 
minimization. 


Table 1 Weight breakdown for 1000-lb payload 


Mise. Prop Gross Burnout Thrust, Burning 
Stage weight weight weight weight 2 mg time 
1 1540 6160 23,800 17,640 47,600 37 
2 3065 9195 16,100 6,905 32,200 133 
3 1136 1704 3,840 2,136 7,680 172 


Fig. 1 shows optimized sizes for several propellant combina- 
tions vs. structural factors for a three-stage rocket when the 
structural factors are such that K. = K, + 0.01 and K; = Ke 
+ 0.01. It will be noted that the missile gross weight is not 
as sensitive at low structural factors as at higher values, and 
the higher the specific impulse the less the effect of structural 
factor for a value of G = 30,000. 

Another conclusion may be drawn by an examination of 
Equation [2]. It will be seen that the burning times can be 
distributed among the stages in any fashion and not change 
the final velocity. In the sample problem the thrust was ar- 
bitrarily made equal to twice the stage gross weight which re- 
sulted in unequal burning times. The reverse procedure 
could have been taken. That is, the total burning time could 


have been divided equally between the various stages. This 
would then set the engine thrusts by the familiar relation 


Of course, the entire procedure is based upon an accurate 
guess for the value of S” g sin 6dt. The value of the integral 
can usually be evaluated with some accuracy by a study of the 
planned trajectory. As pointed out previously, it will lie be- 
tween 0 and gt,. In familiar trajectories the value can be 
guessed with fairly good accuracy. 

In addition, the value of G can be increased slightly, if de- 
sired, to take into account an estimated drag loss. Thus, it 
can be seen that the inclusion of the integral, while trouble- 
some, still permits a degree of flexibility in the use of the opti- 
mization procedure. 
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Escape From a Circular Orbit Using Tangential Thrust 


D. J. BENNEY' 


Canterbury University College, Christchurch, New Zealand 


In this article we consider the problem of escape of a 
rocket from the gravitational field of a planet. Initially 
the space ship is taken to be in a circular orbit about the 
planet. The motors are then switched on and the thrust 
is taken as tangential to the direction of motion at each 
instant. Clearly the motion will remain in the plane of the 
circular orbit. It has been shown that a tangential thrust 
should be a close approximation to the optimal escape 
program.? In calculating the mass ratio it has been 
found convenient to use the r-s equation of the orbit, 
where r is the radius vector and s the length of path. This 
is most useful because the escape condition is very simply 
expressed in terms of s. We conclude by comparing the 
results we obtain for tangential thrust with those for 
circumferential thrust, which were found by Tsien.* We 
will find that a tangential thrust program leads to a 
slightly smaller value for the mass ratio. In particular, 
two cases are considered in some detail, namely, those of a 
large thrust and of a small or micro thrust. 
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Nomenclature 


constant rocket thrust per unit mass 

distance from the center of attraction to the rocket 
distance traversed by the rocket 

polar angle, measured from the initial line 

angle between the radius vector and the tangent to the orbit 
angle between the initial line and the tangent to the orbit 
velocity of the rocket 

time measured from take-off 


gravitational force per unit mass 


radius of curvature 
radius of the initial circular orbit 


BIE ne BG oe 


HE initial conditions att = Oarer = p=a,s=6=0, 
o=y=7/2. We shall require the following results of ele- 
mentary calculus 
dé ds 


sin @ = r— 


dr 
cos @ = — 
ds 


Equations of Motion 


During its flight the space ship is under the action of two 
forces, namely, the gravitational attraction and the tangential 
thrust, as indicated in Fig. 1. 
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Fig. 1 Diagram showing the variables used in the problem 


Resolution along the tangent and normal leads to the equa- 
tions 


For convenience we introduce the corresponding dimension- 
less quantities by replacing r by r/a, s by s/a, and writing 
f = Ag where g = yu /a*, so that g is the acceleration due to 
gravity at the point of projection. Distances are now in 
units of the initial radius a, velocity is in units of initial circu- 
lar velocity v,, time is in units of a/v; and ) represents the 
ratio of thrust to local weight. The above equations then 
become 


dv 1 dr 
ds r? ds 
1 dé 
p r ds 
The initial conditions now att = Oarer =v = 1,s =0=0, 


@= = 2/2. 
Equation [1] may be integrated immediately and using the 
initial conditions yields 


At the escape position the orbit will be parabolic, and so if we 
use the subscript E to denote escape values we have 


ve? = [4] 
TE 
and therefore 
1 
= Dr ee via vac ke [5] 


We now solve Equation [2] for v? in terms of r and the 
derivatives of r with respect to s. We have 


rds dy do as 


Therefore 


Equating values of v? from Equations [3 and 6], and letting a 
prime denote differentiation with respect to s, we obtain the 
differential equation of the orbit 
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r r(1 —r’? — rr”) 


[(2As — 1)r + 1} [1 — r’2] = [(2As — 1)r + 2)]rr”.... [7] 


This differential equation has to be solved subject to the 
initial conditions that when s = 0,r = l andr’ = 0. 


The General Solution 
We attempt to find a series solution for r in powers of the 
parameter \ by assuming the expansion 


r=1+) 


where each r; is a function of s, and where the boundary condi- 
tions require that when s = 0, 7; = r;’ = 0. Substituting this 
series into Equation [7] and equating the coefficients of \ and 
d? we obtain 
mn” +r = 2s 
re” + re = — 8) 


The solutions of these differential equations are 


= 2(s — sin s) 


1 1 
ro = 4 (« s? cos s — 58 sins + 2 cos s — 2) 

It is readily seen that each r; is the solution of a differential 
equation of the type 


where 7;(0) = 7;/(0) = 0. These equations can be solved in 
turn by elementary means; but clearly the solution obtained 
will be of doubtful convergence at escape where s = 1/2X. 

Therefore, we consider two special cases according to 
whether A is large or small. 


Approximation for Large \ 


Let x = 2As. Then initially where s = 0, 2 = 0; and at 
escape where s = 1/2A, x = 1. On making z the independent 
variable in Equation [7] we obtain 


2 2 
dx dx? 


We assume a series solution for r in inverse powers of \?, 
namely 


with the initial conditions r; = dr;/dx = 0, when z = 0. Sub- 
stituting this series into Equation [8] and equating the co- 
efficients of the constant and of \~? we obtain 


x ax? (9) 
dr, \? dr, 


On integrating these equations we obtain the solutions 
= — Que“ — 1] 
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2u 
where u = log(1 + 2). 


A similar procedure could be used to calculate 73, 74, etc. 
To calculate the mass ratio we use Equation [3], i.e. 


2 


Putting 7’ = Xt this becomes 


Therefore 


2), + Mi+a2 
+ 3r? + | [11] 
Let 
a ae 1 
Tr =a + [12] 


On substituting for r; and r2 a somewhat lengthy calculation 
yields 


a = V2 —1 = 0.414212 


48 V2 - log 2 = 0.001615 
104 + | + 
7 25 9 
1 24 4. 104 17 
4 ‘ 5 9 2 
5| 8 48 16 
ie. a, = —0.000064 


Therefore 


0.414212 + + 5, 113] 


If Mp, is the initial mass of the rocket and fuel, then 


M qa 
log ve Tr 
E 


c 


where c is the exhaust velocity. We compare this result with 
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that for circumferential thrust obtained by Tsien,* where for 
large 


0.002349 0.000048 1 
Tg = 0.414212 + = + 0 


This shows that a somewhat smaller mass ratio is required for 
escape when a tangential thrust is employed. 


Approximation for Small \ 


Perhaps of more practical interest is the problem of escape 
from a circular orbit by means of a microthrust. We now ob- 
tain a rather crude approximation for this case, i.e., when X is 
small. Intuitively r’ and r” will both be small for the first part 
of the orbit. We simplify our equation by assuming that r” 
= 0 throughout the motion. The approximate solution of 
Equation [7] is 


This satisfies the boundary conditions to order \. 
Equations [5 and 7] we see that at escape, 1 — r’? = 2rr”. 
The approximation given by [14] then yields 


1 
(20)1/4Q 1/2 
As ds? = dr? + r°d@?, we have 


8E 
(1 — 2As)ds 
0 


So the approximate number of revolutions to escape nz is 
given by 


] 
SE [1 — (20)! 


Therefore 


1 

Further as 

2 

(*) =2\s+--1 
dt r 
we find that 
Te = Mg = 1 — [15] 


Thus, the mass ratio for escape when X is small is given by 
log (% Vga (1 — (20)'8\1/4) 
Me c 


Tsien found that 7, 1 — for circumferential thrust; 
so that again we obtain a somewhat smaller mass ratio for 
escape. 
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The Optimization of Nozzle Area Ratio for Rockets 
Operating in a Vacuum 


M. GOLDSMITH! 
The Rand Corporation, Santa Monica, Calif. 


The correct choice of nozzle area ratio should be con- 
sidered in rocket-powered missiles even when intended to 
operate at high altitudes, where ambient pressure may be 
neglected. Under such conditions, specific impulse of the 
propellant increases monotonically with increasing area 
ratio (exit area/throat area), but the increased weight of 
larger nozzles degrades performance. Hence an optimum 
area ratio can be sought which will provide maximum per- 
formance for the stage in question. Gains in stage velocity 
increment due to correct choice of area ratio are usually 
only a few per cent, which may be of importance in some 
applications. The present investigation results in a linear- 
ized analytic approach to the problem of calculating the 
the optimum area ratio. Design charts are presented. 


Nomenclature 
Ae = nozzle exit area 
Ag = nozzle area at start of exit cone 
At = nozzle throat area 
Cr = nozzle thrust coefficient 
c* = propellant characteristic velocity 
Ts» = specific impulse 
gq = sea-level acceleration of gravity 
k = Ajtp/sin a 
m, = weight of inert parts (total burnout weight) 
Mn = weight of exit cone 
My, = weight of propellant 
ms = inert (burnout) weight of vehicle (without nozzle exit 
cone) 
ms’ = Ms — keg 
t = nozzle wall thickness 
a = nozzle exit half-angle 
y = ratio of specific heats 
AV = stage velocity increment 
= Ag/As 
€0 Ao/At 
= + mp) 
vo = Mp/(ms’ + mp) 
p = nozzle wall weight density 
Subscript 


opt = optimum 


Introduction 


N STUDIES of missile systems in which one or more rocket- 
powered stages operate at very high altitude, the question 
of correct. choice- for nozzle area ratio always arises. Where 
the ambient pressure is sufficiently small that it can be neg- 
lected, elementary rocket theory shows that specific impulse 
of the propellant increases monotonically with increasing 
area ratio (exit area/throat area). However, the larger 
nozzles impose a weight penalty on the missile, which will 
degrade over-all performance. Therefore, an optimum area 
ratio can be sought which will provide maximum performance 
for the stage in question. 


Received July 15, 1957. 
1 Engineer, Aeronautics Department. Mem. ARS. 


Analysis 


The velocity increment for the stage is a useful criterion of 
performance. For a rocket with constant effective exhaust 
velocity (constant specific impulse) during burning, the follow- 
ing well-known equation may be written for the stage velocity 
increment 


where J,, is specific impulse, g is sea-level accleration of 
gravity, and p is the ratio of propellant weight to gross weight 
of the stage. The effects of gravity and drag are ignored in 
Equation [1]. 

For increasing values of ¢, the nozzle area ratio, the value of 
/,, will increase; likewise, the value of v will be decreased, 
owing to the increased weight of the nozzle. Examination of 
Equation [1] then reveals that some optimum value of « 
should be chosen to obtain a maximum velocity increment. 
The purpose of this note is to develop a systematic method of 
determining this optimum area ratio. 

To find the maximum velocity increment, using Equation 
[1], we differentiate with respect to e, all other variables being 


considered as constant, and set the differential equal to zero. 
_2AV AV dy 
de ov de 


Therefore, for best performance 


dI dp 
In (1 v) ( rs ) (2) [2] 


From rocket theory, it is shown that 
= — Ce 


where c* is a characteristic constant of a propellant system, 
and C, is the over-all thrust coefficient. Therefore, Equation 
[2] may be rewritten as 
din Cr dv/de 
de (1 —v)In(1 —») 


Unfortunately, it is not possible to obtain an analytic ex- 
pression for In Cy in terms of the area ratio. Hence, the 
function d In Cp/de must be computed numerically. 

Now, an analytic representation for y may be found. Con- 
sider a rocket nozzle with a conical exit section of half-angle 
a, and constant wall thickness ¢. The parameter v is defined 
as 


Mp 


Mm, + my 


where m, is the weight of propellant carried, and m, is the 
weight of inert parts (structure, payload, guidance, etc.). 
Let us break m; into two parts, the weight m, of the nozzle 
exit cone aft of section Ao, where the conical exit begins, and 
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m,, the weight of all the rest of the inert parts. 
It is easily shown that m, is accurately expressed by 


= (A, — Adip—— 
sin @ 


where ¢ is the thickness of the nozzle wall, p is the density of 
the material, and a@ is the exit cone half-angle. Then write 


where A, is the throat area. Therefore 


Ad 
ms + my, + aod (€ — €) 
sin @ 


For convenience, the following quantities are defined 


Ajtp Mp 
ms’ = m, — key 
sin @ m,’ + Mp 


Introducing the expression for v given in Equation [5] into 
Equation [3], the following relation is obtained 


d In Cer = 
de 


+ my + ke)? 


| 1 + (k/m,’e | 1 + (k/m,')e |} 
{1 + [k/(m,’ + m,) Je 1 + [k/(m,’ + m,) Je 


Now, if the quantity ek/m,’ < 1, ie., if the weight of the 
nozzle exit cone is much less than the remainder of the inert 
weight, the following approximate expression, better suited 
for computation, can be obtained from Equation [6] 


din Cr Yok 


=~ 


de {In (1 — vo) 


2 
| (— + m, 


If the function d In C;/de could be represented with suffi- 
cient accuracy by a polynomial] of only order two, then an ex- 
plicit expression could be obtained for the optimum value of 
area ratio. Unfortunately, it is necessary to use a higher 
order polynomial to represent this function with sufficient 
accuracy, and the solution of Equation [7] is best found by 
graphical means. 

The functions of Cy vs. € have the ratio of specific heats, y, of 
the combustion gas as a parameter. Fortunately, many rocket 
propellants have an effective y ranging from 1.20 to 1.25, and 
using a value of y = 1.23 will give a sufficiently accurate 
representation of Cy for most examples. Fig. 1 shows 
a curve of d In Cp/de vs. € prepared for this value of y. 

Curves showing the right-hand side of Equation [7] are best 
parameterized by the value of vo, which is close to v, the ratio 
of propellant weight to gross weight, and k/m,’, which relates 
nozzle exit cone weight to the weight of the rest of the inert 
parts of the rocket. If Equation [7] is rewritten. the utility 
of these parameters can be seen 


de In [1/(1 — m%)] \m,’ ms! 
[20 ~») + (>. 
In — J § 


Hquation [1] can also be rewritten in terms of these parameters 


t+ — ke]. 


m,' In (1 — 


= c* 
AV = c*Cy In {1 
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Fig. 1 Thrust coefficient Cy and d In Cp/de vs. area ratio « for 
y = 1.23 and expansion into vacuum 
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Fig. 2. Optimum area ratio «,,, as a function of » and k/m’s; 


In Fig. 2, the values of optimum area ratio €o:, as deter- 
mined by graphical solution of Equation [7a], are given in 
terms of v and k/m,’. Thus, the designer can enter this 
figure directly with his own design numbers to obtain the 
proper value of é . Substituting this value of €,: into 
Equation [la] yields the velocity increment obtainable for 
optimum nozzle design. 


Examples 


Two numerical examples will be given to demonstrate the 
utility of this method, and the results obtained from the 
approximate treatment will be compared with the results of 
exact calculations. 


CasE 1 
mp = 8000 Ib 
m, = 2500 |b 
At = 30 in.? 
p = 0.275 lb/in.* 
t = 0.30 in. 
a = 15 deg 
€o = 2 
Therefore 
30 X 0.275 X 0.30 
= — = 9.56 lb 
0.259 
m,’ = 2500 — 2 X 9.56 = 2481 lb 
8000 
= = (0.763 
” 3000+ 2481 
171 


[4] 
0.016} t t 16 
Ss 
A lp 
Mn = (€ — 
sin @ 0.008 + t ’ 14 
2 
_ 
| 
: 
[1a] 


= 0.00382 


€opt = 24.6 from Fig. 2 


Next, to afford a check with an exact calculation, values of 
AV/c*, the performance parameter, have been computed for a 
range of ¢, using Equation [la]. The results are shown in 


Fig. 3(a). 
CasE 2 
Mp = 26 1b 
m, = 9.0 lb 
At = 1.5 in.? 
p = 0.275 lb/in.*. 
= 0.1 in. 
a = 15deg 
€o = 2 
1.5 0.1 0.275 
0.259 
m,’ = 9.00 — 0.159 K 2 = 8.68 lb 
= 0.0184 
Ms 
26.0 
= ——_ = 071 
26.0 + 8.7 
€opt = 8.2 from Fig. 2 


Examination of Figs. 3(a) and 3(b) reveals, as is to be ex- 
pected, that the choice of area ratio has little effect on missile 
performance when the nozzle weight is a very small portion 
of the burnout weight. The importance of correct choice of 
area ratio increases as the nozzle weight becomes more 
significant. Nozzle weights, in conventional practice, seldom 
become large enough in relation to the rest of the inert weight 
to invalidate the approximation made in Equation [7]. The 
gains in velocity increment realized by correct choice of area 
ratio are usually only a few per cent, which may be important 
in some applications. In other cases, nozzle design may still 
be dictated by structural, cooling, etc., requirements. 
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An Analysis of Fuel-Oxidant Mixing in Screaming 
Combustors 


WILLIAM R. MICKELSEN! 


Lewis Flight Propulsion Laboratory, NACA, Cleveland, Ohio 


Vapor fuel-oxidant mixing is analyzed for standing 
transverse acoustic fields simulating those found in 
screeching or screaming combustors. The additional mix- 
ing due to the acoustic field is shown to be a function of 
sound intensity, oscillation frequency and stream ve- 
locity. The effects of these parameters are shown graphi- 
cally for a realistic range of combustor conditions. The 
fuel-oxidant ratio at various combustor stations is shown 
to have a cyclic fluctuation which is in phase with the 
pressure fluctuations. Possible mechanisms contribut- 
ing te screech and scream are discussed. 
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Nomenclature 


speed of sound, fps 

concentration of oxidant 

activation energy 

frequency, cps 

fuel-air, or fuel-oxidant, ratio by weight 
gravitational conversion factor 

Bessel function, first kind, order n 
Lagrangian scale of turbulence, ft 
sound pressure fluctuation, psi 

ambient mean static pressure, psia 


root mean square of sound pressure fluctuation, psi 


= root mean square of sound pressure fluctuation at 


pressure antinode, psi 
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mn, psi 
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USION 


- N+ - + 
LONGITUDINAL RADIAL TANGENTIAL 
Fig. 1 Modes of oscillation in closed-end cylinders 
R = radius of cylindrical duct, ft 
= gas constant 
R = Lagrangian turbulence correlation coefficient 
r = radial coordinate in cylindrical coordinate system, ft 
= flame temperature 
l = residence time in fluid stream, sec 
l = stream velocity, fps 


v = resultant particle velocity, fps 
= resultant particle velocity at duct axis, fps 


Ve 

Ur = radial particle velocity, fps 

 - = tangential particle velocity, fps 

Vor = root mean square of particle velocity, fps 
Vi = intensity of turbulence, fps 

we = fuel flow rate from point source, lb/sec 


X, Y, Z = position of particular point in rectangular coordi- 
nates, ft 


r = axial coordinate in cylindrical coordinate system, ft 

y = lateral displacement of mixing wake center, ft 

Ynax = maximum lateral displacement of mixing wake center, 
ft 

7 = specific heat ratio 

6 = angle in cylindrical coordinates 

pi = density fluctuation at point source, lb/cu ft 

po = mean ambient density, lb/cu ft 

T = time scale of acoustic oscillations, sec 

® = turbulent spreading coefficient, sq ft 

Introduction 


HE effect of acoustic fields on vapor mixing may be of 

importance to the combustion process in screeching or 
screaming combustors. If the acoustic oscillations change 
the fuel-air or fuel-oxidant ratio, the combustion process will 
be directly affected. To illustrate the effect of acoustic oscil- 
lations, an analysis (1)? has been made of vapor mixing in 
cylindrical ducts under conditions found in screeching ramjets 
and afterburners and in screaming rockets. 

Acoustic oscillations in cylindrical chambers have three 
general modes: longitudinal, radial and tangential. These 
modes are illustrated in Fig. 1, where pressure contours are 
shown at an instant in time. The pressure maximums alter- 
nate with sign in time, and higher modes can exist for each 
type. Oscillations in ramjet combustors, afterburners, and 
rockets may be any combination of these types, depending 
on a number of acoustic, aerodynamic and combustion factors 
(2-5). 

The pressure fluctuations in each of the modes of oscillation 
give rise to acoustic particle velocities which can affect the 
fuel-air or fuel-oxidant mixing process in a manner somewhat 
similar to that of turbulent velocity fluctuations. It has 
been shown (6) that the action of a transverse acoustic par- 
ticle velocity is to produce a sinusoidal-shaped deformation of 
the concentration pattern downstream from a source of heat 
or mass, as illustrated in Fig. 2. The analysis showed that 
if the deformation is sufficiently large, the concentration of 
heat or mass in the mixing wake will be changed. The pres- 
sure fluctuations encountered in screeching (3) and screaming 


? Numbers in parentheses indicate References at end of paper. 
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SOURCE 


STREAM 


TRANSVERSE VELOCITY 


Fig. 2. Deformation of mixing wake by standing transverse 
acoustic wave 


(5) combustors are sufficiently large to imply particle veloci- 
ties of the same order of magnitude as the stream velocity. 
Mixing wake deformations caused by acoustic fields of these 
strengths are great enough to cause a considerable change in 
the fuel-oxidant ratio throughout the mixing region. 


Properties of the Acoustic Field 


Experiment has indicated that screeching and screaming 
combustors have acoustic fields approximated by the closed 
end cylinder case. If the ends are acoustically hard, and the 
gas is acoustically homogeneous, the transverse modes of 
oscillation are preserved throughout the length of the cylinder. 
In realistic combustors, the acoustic field may depart from 
the ideal, uniform mode. Data from a simulated afterburner 
(3) show an attenuation and phase shift of the tangential mode 
upstream of the flame zone. A theoretical analysis (7) of a 
similar combustion system has predicted an attenuation up- 
stream, but no phase shift. The attenuation in such a case 
is due to the change in speed of sound between the burned 
and unburned gases, and to a weak closed-end condition at 
the combustor inlet. The present analysis will be confined to 
the ideal, uniform, standing, tangential mode and to pressure 
amplitudes in the acoustic range. 

The general solution of the wave equation (8, p. 398) pro- 
vides an expression for the rms pressure fluctuation at any 
point in the first pure standing tangential mode 


52 
V cos R 


where 7/ p? is the rms pressure fluctuation at the cylindrical 
coordinates (x, r, 0), VW. p3 is the rms pressure fluctuation at 
the pressure antinodes, J;( ) is the Bessel function of the first 
kind and first order, and R is the chamber radius. The 
pressure parameter given by Equation [1] is shown in Fig. 3 
as a function of the fraction of duct radius. 

The particle velocity is related to the pressure fluctuation 


by (8, p. 295) 


where v is the particle velocity, g, the gravitational conversion 
factor, p the fluid density, and the oscillation frequency & is 
given by (8, p. 398) 

where a is the speed of sound. Noting that g/(ap) = a/(ypo), 
where po is the mean pressure of the fluid, the particle velocity 
components can be derived (1) from the general solution of the 
wave equation 


V 13 R 
a Vp sin 0 r R 


YPo Tr r 
of 1.841 — } — J2{ 1.841 — } }.... [5] 
av cos R R 
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where vs and »v, are the tangential and radial components, 
respectively, and Jo( ) and J2( ) are Bessel functions of the 
first kind and zero and second orders, respectively. The par- 
ticle velocity components given by Equations [4, 5] are shown 
in Fig. 3 as a function of the fraction of duct radius. As illus- 
trated in Fig. 3, the resultant transverse particle velocity re- 
duces to a simple expression at the chamber axis 


7 Po 
PRESSURE 
CONTOURS 


° 
71 


cos 6 
2h PR 
cos 
0 2 4 6 8 1.0 
FRACTION OF DUCT RADIUS, r/R 


PARTICLE VELOCITY PARAMETERS 


PRESSURE FLUCTUATION AND 


Fig. 3 Pressure fluctuation and particle velocity profiles in 
first pure standing tangential mode 
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SOURCE CYCLE, st 


DIMENSIONLESS TRANSVERSE 
POSITION, y/[W2 


ie) 5 1.0 5 2.0 
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Fig. 4 Position of mixing wake center at various times through 
one cycle 
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Ne | 
| FREQUENCY, 
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Fig. 5 Position of fuel-oxidant ratio profiles. (a) Fr = 0.50; 
(b) Fr = 0.75 


Mixing Wake Displacement 


A previous analysis (6) provides equations describing the 
motion of the mixing wake caused by a plane transverse sound 
wave. The analysis assumes that the rms particle velocity 
V/ y2 is constant in the region encompassing the deformed 
wake. The particle velocity profiles in Fig. 3 suggest that 
this condition is approximately met if the source is located on 
the chamber axis, and if the wake is not displaced too far 
from the chamber axis. For vapor injection on the chamber 
axis, then, the displacement of the mixing wake center is (‘)) 


y = V2 Ve ain ( sin — ) 


where y is the displacement of the wake center from the 
chamber axis at the axial station x, 7 is the time scale of the 
particle velocity fluctuation, and U is the stream velocity. 
The oscillatory motion of the mixing wake is illustrated in 


Fig. 4, where the dimensionless displacement y/ 2 V 
(1) | is plotted against distance downstream from the inje«tor 
for various times through one cycle of the oscillation. As 
shown in Fig. 4, the transverse motion of the mixing wake is 
cyclic, and produces definite nodes at whole number multiples 
of the dimensionless axial distance $r/UT. he maximum dis- 
placement of the wake center occurs at the wake antino:les, 
and is defined by ymax = V 

An estimate of fidelity to the condition of uniform particle 
velocity in the mixing region may be made by expressing the 
wake displacement in terms of pressure amplitude. This ex- 
pression is obtained by substituting Equations [3, 6] into the 
definition of ymax 


ymax _ 1.32 V pp 

R Y Po 
Comparison of Equation [8] with Fig. 3 shows that up to 
moderately high acoustic amplitudes 


0< Vpi/p <02 


the assumption of uniform particle velocity is fairly good for 
the case of vapor injection on the chamber axis. 


Mixing From Source on Chamber Axis 


The equation for concentration downstream from a point 
source emitting fuel vapor at the weight flow rate w, into a 
gas stream of homogeneous turbulence (9) may be written 
in terms of fuel-oxidant ratio 


pof/A 1 ¥24-29)/4e 
2 ZZ 9 

wy (4mr)*/? Jo wi 
where po is the mean density of the gas, f/A is the fuel-oxidant 
(or fuel-air) ratio at the rectangular coordinates (X, Y, Z), ¢ 
is the residence time, and w is the turbulent spreading coeffi- 
cient. The spreading coefficient w is a function of the time 


t, the turbulence intensity V vb and the Lagrangian correla- 
tion coefficient®. Assuming that 


= exp 03.) 


where £ is the Lagrangian scale of turbulence, the equation 
for the turbulent spreading coefficient [9] becomes 


(10) 


The effect of the transverse particle velocity may be in- 
cluded in the mixing equation by accounting for the trans- 
verse displacement of the wake center (6) 
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where y is given by Equation [7]. 

To investigate the effect of oscillation amplitude and fre- 
quency, fuel-oxidant profiles were calculated from Equation 
{11] by numerical integration with a high speed digital com- 
puter using Simpson’s rule. The profiles were calculated for 
a fixed axial station X = 1.25 ft and for two instants in time 
gr = 0.50 and 0.75. The turbulent stream properties were 


fixed at U = 500 fps, Wv2/U = 0.05 and & = 0.1 ft. The 
profiles were calculated for a range of particle velocities and 
oscillation frequencies. The profile planes were positioned 
normal to the wake center as shown in Fig. 5. The fuel- 
oxidant profiles are shown in Fig. 6. As the particle velocity 
is increased, the peak fuel-oxidant ratio decreases in all 
cases shown in Fig. 6. In addition to this general decrease 
in fuel-oxidant ratio, a cyclic change in fuel-oxidant ratio 
occurs since the profiles in Fig. 6(a) will be repeated at $7 = 
0, 1.0, 1.5,.. . and those in Fig. 6(b) at $7 = 0.25, 1.25, 1.75, 

As shown previously (6) the fuel-oxidant ratio will be high- 
est in portions of the wake having zero slope (dy/dx) and 
lowest in portions having maximum slope. This indicates 
that the fuel-oxidant ratios in Fig. 6(a) for ¢7 = 0.5 will form 
the upper limit, and those of Fig. 6(b) will form the lower 
limit of the cyclic variation. 

[t has been shown (1) that a time and position relation exists 
between the fluctuations in fuel-oxidant ratio and the pres- 
sure. When the wake is at its maximum lateral displacement, 
the fuel-oxidant ratio at the antinodal stations is at a maxi- 
mum, and the pressure is also at a maximum. This relation 
is expressed by an equation which defines the pressure change 
experienced by the wake center during its transverse motion 


YPop 
== —2.95 sin | —— } sin | — —— cos (2757) 
Pr U l 


The pressure fluctuation given by Equation [12] is com- 
pared in Fig. 7 with values of fuel-oxidant ratio obtained from 
Fig. 6. The comparison is made for two sets of axial stations 
= 0.50, 1.50, ... and 0.75, 1.75, .... At the antinodal 
stations ¢X/U = 0.50, 1.50, . . . , the pressure fluctuations 
are exactly in phase with the fuel-oxidant ratio fluctuations. 
At stations away from the antinodal, such as #X/U = 0.75, 
1.75, ..., the pressure and fuel-oxidant ratio fluctuations are 
no longer in phase, but the magnitude of the pressure fluc- 
tuation is considerably reduced. 

Fuel-oxidant ratio profiles were also calculated from Equa- 
tion [11] for a range of stream velocities and particle velocities. 
The profiles were calculated for a fixed axial station of XY = 
1.25 ft and for two instants in time $7 = 0.50 and 0.75. The 


turbulence properties were fixed at Vv3/ U = 0.05 and £ = 
0.1 ft, and the oscillation frequency was fixed at ¢ = 1000 cps. 
The profile planes were positioned in the wake as shown in 
Fig. 8, and the fuel-oxidant profiles are shown in Fig. 9. The 
percentage reduction in peak fuel-oxidant ratio with increas- 
ing particle velocity is greatest for the lowest stream velocity 
shown in Fig. 9. Comparison of Figs. 6 and 9 shows that 
the effect of stream velocity is similar to that of oscillation 
frequency, but more pronounced. A cyclic variation in 
fuel-oxidant ratio with time is also evident in Fig. 9. 


Injection Away From the Chamber Axis 


In the case of vapor injection at points away from the duct 
axis, pressure fluctuations at the injector may cause cyclic 
changes in the fuel-oxidant ratio. The mixing equation may 
be corrected for this effect by noting that the fuel-oxidant 
ratio of a fluid particle leaving the vicinity of the injector is 
proportional to po/(po + p;) where po is the ambient time- 
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Fig. 6 Effect of oscillation frequency on mixing in first pure 
standing tangential mode. Stream velocity, 500 fps; turbulence 
intensity, 5 per cent; turbulence scale, 0.1 ft; axial station, 
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Fig. 7 Time phase relation between pressure and fuel-oxidant 

fluctuations on wake center during transverse motion. (a) At 

antinodal stations |X /U = 0.5, 1.5,...; (b) at stations = 
0.75, 1.75, ... 
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Fig. 8 Position of fuel-oxidant ratio profiles. 
(b) Fr = 0.75 
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Fig. 9 Effect of stream velocity on mixing in first pure standing 
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Fig. 10 Mixing in first pure standing tangential mode with 
injection at points not on pressure node 
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Fig. 11 Typical heat release curve 


mean density, and p; is the density fluctuation at the injector, 
This density ratio may be inserted into the mixing equation 


af/A_ fe 


1 1 


= HY — 9)? 4-29) 13 
where p, is the pressure fluctuation at the injector 
= 
— 0 QnFr)..... 
Pi 2.43 V pe ) os cos (24Fr) 14] 


With injection away from the chamber axis, the solid body 
motion expressed by Equation [7] no longer obtains. The 
pressure gradient encountered by the wake exerts an alter- 
nately compressive and expansive motion in the transverse 
direction. Along the line @ = 0, the exact expression for 
particle motion is 


d(y/R) 


0.66 PR 
Jo (1.81 wi x) 


R 


sin(2m$r)dr..... 15] 


where y; and y2 are the particle locations at the times 7; 
and 72, respectively. Due to the complexity of Equations 
[14, 15], fuel-oxidant profiles were not calculated for the case 
of injection away from the chamber axis. By inspection of 
Equation [14], a qualitative picture of the effect of pressure 
fluctuation at the injector can be formed as shown in Fig. 
10 where the fuel-oxidant wake is shown for successive in- 
tervals of time, along with the pressure fluctuation curve. 
The width of the wake represents the relative magnitude of 
the fuel-oxidant ratio. The fuel-oxidant ratio at the in- 
jector is high when the pressure fluctuation is negative and 
is low when the pressure fluctuation is positive. The figure 
shows that when the richest portion of the wake has reached 
its maximum excursion toward the chamber wall, the pressure 
fluctuation has reached its maximum positive value. In this 
way, the pressure fluctuations at the injector tend to accen- 
tuate the fuel-oxidant ratio fluctuations caused by wake 
deformation. 


Possible Mechanisms Contributing to Combustor 
Oscillations 


The role of pressure fluctuations in driving screech has been 
discussed elsewhere (3), where it is suggested that a periodic 
change in pressure is expected to cause an even larger periodic 
change in heat release rate. Consideration of the chemical 
reaction rate equation indicates that fuel-oxidant fluctu- 
ations could also cause fluctuations in heat release rate. For 
example, the second order reaction expression may be written 


reaction rate ~ C4 f/A e~ 


where C, is the oxidant concentration, F is the activation 
energy, R, is the gas constant, and T is the flame temper- 
ature. The flame temperature is a strong function of the 
fuel-oxidant ratio so that the heat release has the form shown 
in Fig. 11. 

Unpublished data taken in a typical combustor indicated 
that screeching occurred at over-all fuel-air equivalence ratios 
between 0.25 and 0.45 based on unburned inlet air. In the 
absence of screeching, the peaks in the fuel-air ratio profile 
at the flame front may be expected to be lower than stoichio- 
metric for such lean over-all fuel-air ratios, as shown in Fig. 
12(a). As pointed out previously, transverse acoustic oscil- 
lations cause a general decrease in peak fuel-air ratio, and also 
cause a periodic fluctuation in fuel-air ratio as illustrated in 
Fig. 12(b). In the sketch, the solid line represents the upper 
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limit of the fuel-air ratio fluctuation, and the dotted line the 
lower limit. Since the mean value of the peak fuel-air ratio 
is lower than stoichiometric, then a fluctuation toward higher 
fuel-air ratio will result in an increased heat release rate as 
shown by Fig. 11. At antinodal distances ¢X/U downstream 
from the spray bar, the increase in heat release due to the 
fuel-air ratio rise is directly in phase with the increase in heat 
release due to pressure rise, so that screech might be driven 
to even higher amplitudes. If the mean peak fuel-air ratio 
were higher than stoichiometric, then a rise in fuel-air 
ratio would tend to decrease the heat release rate, and thus 
exert a suppressing effect on the screech driving mechanism. 
To obtain a mean peak fuel-air ratio higherthan stoichiometric, 
the number of injection points might be reduced, with a 
higher fuel flow rate through each, resulting in a profile such 
as shown in Fig. 12(c). An alternate method could be to 
move the spray bar closer to the flame front as shown in 
Vig. 12(d). 

The effect of off-center injection on the screaming mecha- 
nism is not immediately apparent. It was shown earlier that 
the pressure fluctuation at the off-center injector tends to 
accentuate the fuel-oxidant fluctuations caused by the de- 
formation of the wake by particle velocity. Since the par- 
ticle velocity encountered by the wake decreases as the point 
of injection is moved away from the chamber axis, the fuel- 
oxidant ratio fluctuations caused by wake deformation will 
also decrease. Conversely, as the injector is moved away 
from the chamber axis, the contribution to fuel-oxidant fluc- 
tuation by the pressure fluctuation at the injector should 
increase. In addition to these two effects on fuel-oxidant 
ratio, the mixing wake will encounter more severe pressure 
fluctuations during its transverse motion downstream from 
an off-center injector. Resolution of these three effects 
awaits further analysis or experiment. 

One more possible contribution to the screech mechanism 
in the first pure standing tangential mode deserves mention. 
In the case of ramjet combustors and afterburners where 
injectors are positioned upstream of flameholders, the axial 
distance between the injector and flameholder stations may 
be of importance to the screech mechanism. As shown by 
Fig. 4 and Equation [7], the mixing wake may be periodically 
displaced a considerable lateral distance at antinodal axial 
distances from the point of injection. If the injector to flame 
front distance is at or near a wake antinodal distance, then 
the fuel-air mixture might pass alternately back and forth 
across a flameholder and cause significant periodic changes 
in heat release in the two opposing flame fronts. The cure 
or control of this type of screech mechanism may lie in the 
injector to flame front distance. By changing this distance 
by a fraction of the internodal distance, the time and position 
phase relation between the fuel-oxidant ratio fluctuations and 
pressure fluctuations might be altered to attain a suppressing 
effect on the screech driving mechanism. 


Conclusions 


In summary, equations have been derived for the calcula- 
tion of fuel-oxidant ratio profiles downstream from injectors 
simulating point sources of vapor in cylindrical combustors 
with the first pure standing tangential mode of acoustic os- 
cillation. From the results of this analysis, the following con- 
clusions are drawn. 

1 Acoustic oscillations of the magnitude present in screech- 
ing or screaming combustors can cause considerable transverse 
motion and distortion of the fuel-oxidant mixing wake. 

2 In general, the peak fuel-oxidant ratio at all axial sta- 
tions is lowered by the transverse acoustic oscillations. 

3 At all axial stations, the peak fuel-oxidant ratio has a 
cyclic variation with time. The phase relation between the 
peak fuel-oxidant ratio and pressure fluctuations is a function 
of the position of the station, the stream velocity and the 
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Fig. 12 Fuel-air ratio profiles in combustor 


oscillation frequency. At the antinodal axial stations, the 
pressure and fuel-oxidant fluctuations are exactly in phase. 


Acknowledgment 


The assistance of Lois T. Dellner of the Lewis laboratory, 
who performed the programming and calculations with the 
high speed digital computer, is gratefully acknowledged. 


References 


1 Mickelsen, William R., “Effect of Standing Transverse 
Acoustic Oscillations on Fuel-Oxidant Mixing in Cylindrical 
Combustion Chambers,’? NACA TN 3983, 1957. 

2 Truman, John C., and Newton, Robert T., “Why Do High- 
Thrust. Engines Screech?” Aviation Age, vol. 23, no. 5, May 
1955, pp. 136-143. 

3 Blackshear, Perry L., Rayle, Warren D., and Tower, 
Leonard K., “Study of Screeching Combustion in a 6-Inch Simu- 
lated Afterburner,’? NACA TN 3567, 1955. 

4 Berman, Kurt, and Logan, Stanley E., ‘Combustion 
Studies with a Rocket Motor Having a Full-Length Observation 
Window,” JeT Proputston, vol. 22, no. 2, March-April, 1952, 
pp. 78-85. 

5 Smith, R. P., and Sprenger, D. F., “Combustion Insta- 
bility in Solid-Propellant Rockets,’ Fourth Symposium (Inter- 
national) on Combustion, Williams and Wilkins, 1953, pp. 893- 
906. 

6 Mickelsen, William R., and Baldwin, Lionel V., ‘“Aerody- 
namic Mixing Downstream from Line Source of Heat in High- 
Intensity Sound Field,” NACA TN 3760, Aug. 1956. Also, 
“Aerodynamic Mixing in a High-Intensity Standing-Wave 
Sound Field,’ Journal of the Acoustical Society of America, vol. 
29, 1957, p. 46. 

7 Moore, Franklin K., and Maslen, Stephen H., ‘Transverse 
Oscillations in a Cylindrical Combustion Chamber,’’ NACA TN 
3152, 1954. 

8 Morse, Philip M., “Vibration and Sound,” 2nd edit., 
McGraw-Hill, New York, 1948. 

9 Mickelsen, William R., ‘An Experimental Comparison of 
the Lagrangian and Eulerian Correlation Coefficients in Homo- 
geneous Isotropic Turbulence,’”” NACA TN 3570, 1955. 


177 


| 
U— 
(a) 
| 
| 
R 
| 
(c) 
| 
T2 
: 


Heat Transfer in Laminar Pipe Flow With Uniform Coolant 
Injection 


S. W. YUAN? and A. B. FINKELSTEIN?® 


Polytechnic Institute of Brooklyn, Brooklyn, New York 


The effect of coolant injection at the wall on the tem- 
perature distribution of the laminar flow of a fluid in a 
porous wall pipe has been investigated. The temperature 
distributions in the pipe for various fluid injection rates 
are obtained by solving the energy equation in which the 
velocity distribution terms used are taken from Reference 
(1). The heat transfer results, expressed as a local Nusselt 
number, are presented in terms of the Péclet number, the 
dimensions of the pipe and the coolant flow rate. The 
relation between the wall temperature and the rate of 
coolant injection is also calculated. The results also 
show that the heat transfer coefficient at the wall decreases 
with an increase in the rate of coolant injection, whereas it 
was previously found that the friction coefficient at the 
wall increases with an increase of fluid injection. Thus 
the usual analogy between heat transfer and momentum 
transfer does not hold in porous wall cooling of pipe flow. 


Nomenclature 


axial distance in the main flow direction 

radial distance normal to the main flow direction 
azimuthal angle in cylindrical coordinate system 
diameter of the circular pipe 

radius of the circular pipe 

length of the pipe 

r/R 

x/R 

fluid velocity in the x direction 

fluid velocity in the r direction 

maximum velocity in the Poiseuille flow (vo = 0) 
average velocity in the Poiseuille flow (U = 2um) 
coolant injection velocity at the wall 

fluid temperature (absolute) 

heat flow per unit area 

pressure 

dynamic viscosity 

kinematic viscosity 

mass density of the fluid 

specific heat at constant pressure 

thermal conductivity 

heat transfer coefficient 

dissipation function 

axial Reynolds number (u,D/v) 

injection Reynolds number 

Prandtl number (ycp/k) 

Péclet number (PrRe) 

local Nusselt number (hD/k) 


Subscripts 

0 coolant fluid 

1 refers to plane z = 0 
w wall of pipe 

M mean value 
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Introduction 


HE transpiration cooling problem, i.e., to cool the com- 

bustion chamber liners by means of injecting coolant 
through a porous wall, has been solved to a great extent 
in the case of a laminar flow over a porous plate (2).4 It 
has been assumed that the results obtained from the porous 
plate solution would yield some indication of the actual 
case in the combustion chambers. On the other hand, how- 
ever, the problem of the nonisothermal flow through a por- 
ous-wall pipe with injection has not yet been investigated. 
The purpose of this work was to obtain basic characteristics 
of this type of flow which in turn would provide guidance 
for the investigation of transpiration cooling in a turbulent 
pipe flow, and to verify whether the above assumption made 
for the plate investigation is valid. 

The problem of heat transfer in a steady laminar flow in a 
circular pipe without injection has been studied by Graetz 
(3) and Nusselt (4) since the latter part of the last century. 
In these studies it was assumed that the wall temperature is 
constant and changes discontinuously at x = 0 (Fig. 1), that 
the physical properties of the fluid are independent of tem- 
perature, and that the Poiseuille velocity distribution is main- 
tained throughout the motion. But up to the present, because 
of the complexity of the problem, there has not been much 
further improvement in the solution beyond that of an em- 
pirical nature. 

In the present investigation an exact solution of the energy 
equation in cylindrical coordinates was obtained in which the 
velocity component terms used were taken from (1). The 
solution gives the temperature distribution as functions of co- 
ordinate axes, Péclet number, the diameter of the pipe and 
the rate of coolant injection. 

The assumptions made in the present study were: 

1 The fluid is incompressible and the free convection, ra- 
diation and dissipation are neglected. 

2 The molecular heat conduction in the axial direction can 
be neglected in comparison with that in the radial direction. 

3 The coolant flowing through the porous wall is uniform 
throughout. 

4 The wall temperature is constant and changes dis- 
continuously at z = 0. 


Fundamental Equations 


If a curvilinear coordinate system is introduced with the 
origin at the center of the cross section of a circular pipe 
where x is taken in the direction of the flow, r in the radial 
direction and ¢ the azimuthal angle, then with axial sym- 
metry of flow the Navier-Stokes equations become 


1 Ou 


4 Numbers in parentheses indicate References at end of paper. 
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Fig. 1 Temperature distribution along a porous circular pipe 


The continuity equation 


Oru)  Arv) 
or or 


The energy equation 
y eq 


oT oT 
dp 10 ar 

The velocity distributions of the fluid flow in a porous-wall 
circular pipe through which uniform injection is applied were 
obtained by the solution of Equations [1, 2, 3] with appro- 
priate boundary conditions. These solutions were given in (1) 
and will be used in the present investigation. 

For an incompressible fluid it can be shown that the terms 
in Equation [4] due to the pressure gradient and the dissi- 
pation can be neglected, and, furthermore, it is assumed that 
0?7’/Ox? may be neglected in comparison with 


Hence Equation [4] can be simplified in the nondimensional 


form 
Udn Udé PrReLnodn\ on 
The boundary conditions are 


From (1), the nondimensional axial and radial velocity com- 
ponents are given respectively for 0<A<1 (neglecting terms 
of order \? with respect to d) as follows 


u 1 
and 
v r 1 
where 


= (1 — 9°) + (—2 + — + + ... 


With these values for the velocity components, Equation 


[5] becomes 
1 r 00 
F = (A718) 


h 1 «1 fla/ 
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Method of Solution 


The energy equation written in the form of Equation [10], 
again neglecting terms of order \? after multiplying both sides 
by 1—(A/18), yields immediately a solution in the form of an 
infinite series 


—Cj2/4Prr 
6= A; (1 C;)....- [11] 


j=l 


where the M; (n,C;) are the solutions of 


aM 1 | dM 
C2 [« — + (74? — + | M = 0... [12] 


and the C; are the eigenvalues of Equation [12] which cor- 
respond to the boundary condition (8), = 1 = 0. The series 
solution of Equation [12] corresponding to the eigenvalue 
C;, which is free from singularities at 7 = 0 is 
C;* 1 [ C;* 
M; =1—— 
i 4 1 al 4 
CF? (pa 36 A+ | ni +...[13] 


Multiplying both sides of Equation [12] by the appropriate 
factor, we can write it as the Sturm-Liouville equation 


d 
ay + C2p(n)M; =0.......... [14] 


where 


m(n) = nePn® — 


p(n) = mn — + — + | 


Hence M; (n, C;:) and M; (n, C;) for i # j are orthogonal 
functions with respect to the weight function p (9); i.e. 


1 
f, = 0 


The coefficients of the series expansion [11], the A;, are 
determined from the boundary condition [7] applied to 
Equation [11], which is 


[15] 


j=1 


Then multiplying both sides of Equation [16] by p(y): 
M; (n, C;) and integrating from 0 to 1 


1 
p(n) Mj(n, 


p(n) M Cy )dn 
From the differential Equation [14], it can be seen that 
1 1 dM; 


and 
1 v 
fy Can = 


2C; dy \OC 


Thus 
2 


The results of the calculations for the C; and A, are given in 
Table 1 for different values of A. 
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Fig. 2 Local heat transfer coefficient along the length of a pipe 


Application of the Solution 


The heat transfer coefficient for the flow in a pipe is us- 
ually calculated with the difference between the mean tem- 
perature of the fluid and the wall temperture. The mean 
temperature over the cross section considered, weighted with 
respect to the axial velocity, frequently referred to as the 
“mixed-mean’”’ fluid temperature, may be defined by 


The heat transfer coefficient h is then defined by 


where the heat flow per unit area at the wall is 


oT k 06 
q= =— — Ty) (2) 


In nondimensional form, the heat transfer coefficient, or 
local Nusselt number Nu is 


™ \ PrRe D 
1 


2 JB 1 + 


The temperature distribution given in Equation [11] may 
be used to determine the amount of coolant required to 
cool the wall to a predesignated temperature. From the 
condition of heat balance at the wall, one obtains 


T 
0 or /r=R 


Using Equations [23, 11] in conjunction with [25], one obtains 
the temperature difference ratio as follows 


AT — To)l....... 


A\My'(1) 

T, — T. 
4Prr 


A2M,'(1) 
C, {1 (1 +4 i) 


~ 


where 


3.3 
Nu 
3:2 
3.1 
3.0 
.75 1.0 


Fig. 3 Local heat transfer coefficient for various rates of cool- 
ant injection; (1/PrRe) (1/D) = 0.075 
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Fig. 4 Temperature difference ratio for various coolant param- 
eters; (1/PrRe) (1/D) = 0.075 


Table 1. Eigenvalues C; and coefficients A; for various 
injection Reynolds numbers A | 
C, Ai As | 
0 2.705 6.66 1.477 —0.810 
0.25 2.765 6.67 1.507 —0.861 
0.5 2.825 6.67 1.540 —0.920 
0.9 2.922 6.67 1.595 —1.020 
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Discussion 


In Fig. 2 the local Nusselt number along the tube length was 
shown for two cases of coolant Reynolds number A (see (5) 
for \ = 0). The Nusselt number, in both cases, starts with 
an infinitely large value at the beginning of the heated tube 
and decreases in flow direction to a constant value asymp- 
totically when (1/PrRe)(1/D) approaches infinity; however, 
this constant value is approximately reached when (1/PrRe)- 
(//D) = 0.06. 

It can be seen in Fig. 3 that the Nusselt number decreases 
almost linearly with an increase of the coolant Reynolds 
number. This is due to the fact that the range of coolant 
Reynolds numbers considered in the present investigation 
is rather small (A S 1). As the coolant Reynolds number 
increases further (A > 1) the Nusselt number then decreases 
more gradually. This phenomenon was also obtained in the 
case of nonisothermal flow over a plate with coolant injec- 
tion. 

The ratio of temperature difference, (7; — Tw)/(Tw — To), 
was plotted against the coolant flow ratio v) l/um D in Fig. 4. 
Yor a predesignated wall temperature, the amount of coolant 
required per unit time can be readily determined provided 
the entrance temperature 7, and the coolant temperature 7 
are known. 

The results obtained from the studies of nonisothermal 
laminar flow over a plate with coolant injection (2) has shown 
that the friction coefficient at the wall decreases with the 
increase of coolant injection. It has also indicated that the 


heat transfer coefficient decreases with the increase of coolant 
injection. 

Thus, there is a definite direct relationship between the 
friction coefficient and the heat transfer coefficient. On the 
other hand the results obtained in (1) indicate that in a 
pipe flow the effect of fluid injection at the wall is to accel- 
erate the main stream velocity. Hence the velocity gradient 
at the wall, which determines the wall friction, increases. 

This phenomenon thus indicates that the analog between 
the heat transfer and momentum transfer does not exist in 
porous wall cooling of pipe flow. In view of the results ob- 
tained in the present investigation it is doubtful that the re- 
sults obtained from the porous plate solution would yield the 
same trend of the actual case in the combustion chambers. 
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Influence of Pressure on the Combustion of Liquid 
Spheres 


GEORGE A. AGOSTON? BERNARD J. WOOD* and HENRY WISE‘ 


Stanford Research Institute, Menlo Park, Calif. 


An experimental study was made of the influence of pres- 
sure on the steady-state mass burning rates of liquid fuels 
coating the surface of porous Alundum spheres of various 
diameters. Observations were made on the combustion of 
a series of liquid organic hydroxy compounds and several 
heavy hydrocarbons at atmospheric and subatmospheric 
pressures. These data in addition to flame distance 
measurements were successfully correlated revealing 
quantitatively the influence of free convection. At sub- 
atmospheric pressures the change in burning rate is so 
modest that for practical purposes the burning rate of 
hexadecane drops (for example) less than 1000 » in diam 
may be considered independent of pressure under free 
convection conditions. 
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Nomenclature 
Li L 

Ce. = specific heat (constant pressure) of air at temperature 7',, 
eal/(g)(° C) 

Cri = specific heat (constant pressure) of liquid at mean tem- 
perature between initial liquid temperature and 7's, 
cal/(g)(° C) 

= drop diameter, cm 

D,. = flame diameter based on flame fadius at forward stagna-~ 


tion point, cm 
Deo = flame diameter (spherical flame) with conductive heat 
transfer in absence of convection, em 


g = gravitation constant, cm/sec? 

Gr = Grashof number = = (Te — Th = 2D3 

T» 
i = stoichiometric oxygen-fuel ratio g/g 
L = heat of vaporization (sensible plus latent heat), cal/g; 

L=L;+ — Te) 

Ls = heat of vaporization at 7's, cal/g 
m = mass burning rate with natural convection, g/sec 
mo = theoretical mass burning rate with conductive heat trans- 


fer in the absence of convection, g/sec 
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Table 1 Steady-state mass burning rate as a function of pres- 
sure and diameter of porous sphere; oxidizer = air at 20 C 


Sphere [Pressure | Burning Rate m/D, 
Fuel Diameter, P, m x 10°, 
cm mm Hg g/sec (g/sec cm) x 10° 
Hexadecane 0.440 47 2.91 6.61 
64 2.91 6.61 
141 2.75 6.25 
250 2.72 6.18 
450 2.76 6.28 
550 3.06 6.96 
650 3.10 7.05 
760 3.41 7.75 
0.940 18 5.81 6.18 
; 28 5.43 5.78 
58 5.43 5.78 
101 5.43 5.78 
187 6.28 6.68 
273 6.82 7.25 
538 7.72 8.21 
760 9.0 9.56 
1.208 17 6.86 5.69 
26 6.86 5.69 
30 6.52 5.40 
62 6.67 5.53 
62 6.87 5.69 
103 7.56 6.26 
220 8.80 7.29 
411 10.75 8.9 
604 11.9 9.85 
604 11.9 9.85 
760 12.7 10.5 
Xylene 1.208 100 16.1 13.3 
181 15.1 12.5 
400 16.1 13.3 
431 20.4 16.9 
760 23.3 19.5 
760 23.5 19.3 
Cente 
Pr = Prandtl number = x 
q = heat of combustion (fuel vapor to vapor products), cal/g 
T, = temperature of ambient air, ° K 
Ts + Ta 
T> = mean temperature, °K; 7, = 


T, = adiabatic flame temperature, ° K 


Ta = log mean temperature, °K; 23 log 
e = log mean temperature, ° K; 7, = 23 log (T./Ta) 
T; = temperature at liquid surface, ° K 
Y, = weight fraction of oxygen in ambient gas, g/g 
Aq = thermal conductivity of air at temperature 7g, cal/(sec) 
(em)(° C) 
Ae = thermal conductivity of air at temperature 7, cal/(sec) 
(em)(° C) 
wa = viscosity of air at temperature 7g, g/(sec)(cm) 
Me = viscosity of air at temperature 7,, g(/sec)(cm) 
pa = density of air at temperature Ja, g/cc 
QQ = diameter-free portion of Grashof number, cm~* 
Introduction 


OMBUSTION in gas turbines at high altitudes is of great 
concern in aircraft operation. One factor which has been 
considered to be of possible significance in combustor per- 
formance is the influence of pressure on the burning of fuel 
sprays. Studies of the combustion of individual fuel drops 
at pressures above atmospheric have indicated that the mass 
burning rate increases with pressure raised to a power be- 
tween 0.2 and 0.5 (1, 2).5 This pressure dependency suggests 
that natural convection plays a measurable role in the 
heterogeneous combustion process (3, 4). An object of the 
present investigation was to provide means for evaluating the 
mass burning rates of individual drops of various diameters at 
pressures less than 1 atm. 


5 Numbers in parentheses indicate References at end of paper. 
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Sphere | Pressure | Burning mate m/D, 
Fuel — (g/eec cm) x 103 

2, 2'-Oxydiethanol 0.312 760 2.46 7.89 

0.440 760 3.58 8.14 

0.638 730 5.38 8.43 

0.940 15 5.82 6.20 

27 6.05 6.44 

27 5.82 6.20 

40 5.55 5.91 

40 6.55 6.96 

40 6.33 6.74 

40 5.55 5.91 

50 5.70 6.06 

53 5.26 5.60 

53 5.55 5.91 

71 5.55 5.91 

71 5.82 6.20 

81 6.05 6.44 

101 6.38 6.78 

101 6.34 6.75 

162 6.77 7.20 

176 6.73 7.17 

302 7.34 7.80 

485 8.40 8.94 

667 9.13 9.71 

667 9.13 9.71 

722 9.40 10.00 

722 9.07 9.65 

724 9.13 9.71 

726 9.40 10.00 

1.208 13 6.85 5.67 

15 6.95 5.83 

15 7.05 5.75 

17 7.61 6.13 


Sphere |Pressure | Burning mate m/D 

Fuel Diameter m x 10 ‘ 
Hg g/eec (g/sec cm) x 10 

2, 2'- iethanol 1, 208 17 7.40 6.30 

_ 27 7.05 6.06 

27 7.34 5.83 

36 7.34 6.08 

64 8.00 6.62 

85 8.40 6.95 

105 8.62 7.14 

105 8.62 7.14 

155 9.15 7.58 

265 10.40 8.6 

423 11.75 9.71 

596 12.50 10. 34 

722 12.50 10.34 

1.590 760 20.4 12.83 

1.906 760 26.6 13.96 

Apparatus 


The use of wetted porous spheres in steady-state hetero- 
geneous combustion studies has previously been described (4, 
5). In this investigation this method was adapted to studies 
at subatmospheric pressures. A diagrammatic representation 
of the equipment is given in Fig. 1. Fuel is fed to the center 
of the porous Alundum® sphere and uniformly wets the sur- 
face. A flame is established by ignition of the fuel with an 
electric spark. The rate of fuel delivery to the porous sphere 
is then manually adjusted to the steady-state burning rate. 

The porous sphere is situated inside a glass cylinder, within 
which the pressure may be regulated by adjusting the intake 
and exhaust valves. An air intake of the order of 6000 em*/ 
min (S.T.P.) is required to maintain a continually fresh supply 
of oxidizer to the region of combustion; with this air flow, a 
high capacity Kinney Vacuum Pump permits the attainment 
of absolute pressures of the order of 20 mm Hg within the duct. 
In all the experiments, air velocities were sufficiently low so 
that forced-convection effects could be neglected (5). 

The influence of oxygen concentration on combustion at 
atmospheric pressure was investigated in the same apparatu: 
by enriching the air with oxygen before admittance to the 
glass cylinder. Burning rates in air at atmospheric pressure 
were obtained with the spheres removed from the cylinder and 
exposed to room air. 


6 The Alundum spheres (catalyst supports) were supplied by 
the Norton Company, Worcester 6, Mass. 


JET PROPULSION 


| 
| 
| 
| 
‘ 


| 


)x 10° 


etero- 
ed (4, 
tudies 
tation 
enter 
sur- 
th an 
phere 
rate. 
vithin 
ntake 
em?/ 
ipply 
OW, a 
ment 
duct. 
yW SO 


on at 
ratus 
o the 
ssure 
rand 


by 


SION 


EXHAUST 
CaCi2 TRAP 
AND VACUUM PUMP 
UPPER 
Hg MANOMETER 
FUEL 
RESERVOIR 
6 INCH DIAMETER | 


GLASS CYLINDER 
POROUS ALUNDUM —}——>(_}——— 
SPHERE 


LOWER FLANGE 


SPARK GAP 


SCREEN 


VOLTAGE 
= TRANSFORMER 
IAG — THREE-WAY 
GLASS WOOL STOPCOCK 
PRECISION 
BAFFLE BORE SYRINGE 
FLOWMETERS 
SERVO- DRIVE 
MECHANISM To ac 
CALIBRATED 
02 ‘iti FLUID DELIVERY 


RATE CONTROL 


Fig. 1 Steady-state burning apparatus 


Results 


For several different sphere sizes, measurements were made 
of the burning rates of 2,2’-oxydiethanol (diethylene glycol), 
hexadecane (cetane), and xylene in air at subatmospheric 
pressures. The data obtained are presented in Table 1. 
Initially the burning rates diminish progressively as the pres- 
sure is decreased, but there is some indication of an increase in 
the burning rates with further reduction in pressure. Such a 
trend is to be expected because of the variation of latent heat 
of vaporization with pressure (see next section). 

Lowering the pressure causes changes in the flame size and 
geometry (Figs. 2, 3). The distance of the combustion zone 
from the drop surface increases with decreasing pressure and 
the flame envelope approaches a spherical shape at very low 
pressures. Also, the thickness of the luminous region ap- 
pears to increase. The distance from the surface of the 
sphere to the outer surface of the visible flame at the forward 
stagnation point (i.e., in this case the bottom-most point) was 
determined from photographs of a number of flames, and the 
results obtained are presented in Table 2. For comparison 
purposes flame distance measurements were also obtained 
from motion pictures of burning fuel drops suspended from 
silica fibers in air at atmospheric pressure, Table 2. 

The mass burning rates of a number of fuels in air at at- 
mospheric pressure were obtained for a series of sphere sizes, 
Table 3. These data are useful in interpreting the influence 
of natural convection (see next section). The typical trend 
of the data is illustrated in Fig. 4, where the burning rate of 
ethanol (95.5 per cent) is displayed as a function of sphere 
diameter. Table 3 also includes burning rate data at 1 atm 
obtained for ethyl alcohol burning in air enriched with oxygen. 


Discussion 


Similarity considerations in heat transfer problems involv- 
ing natural and forced convection have led to dimensionless 
groups which enable the correlation of data obtained under 
varied conditions. This procedure has been applied to evapo- 
rative processes under forced and natural convection condi- 
tions by Fréssling (6), Ranz and Marshall (7), and others. 
The combustion of liquid spheres under forced convection con- 
ditions has similarly been treated analytically as a heat-trans- 
fer controlled process (3). Experimental work has yielded the 
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Table 2 Flame distance at forward stagnation point 


T T 
| Diameter, Pressure, Flame Distance, 
Fuel D, P. (D,.- D)/2, Gr 
| cm mm Hg cm ‘ae a 
| Porous Spheres | 
| 2, 2*-Oxydiethanol 0.940 15 1.44 | 3.83 10.5 0.46 | 
| 15 1.38 3.70 10.5 0.40 
| 27 i 0.992 2.93 10. 25 1.24 
| 50 0.063 2.2? 10.0 4.14 
| 81 0.495 } 1.93 9.8 10.15 
| lez 0. 386 } ult 9.5 35.0 
| 302 0.330 | 1.60 9.21 124.0 
726 0.274 | 1.49 8.75 688.0 
| 1.208 15 1 54 4.29 13.3 0.858 
| 15 1.48 4.16 13.3 0.858 
| 27 1.045 3.28 13.1 2.64 
0.715 | °2.64 12.65 13.6 
o4 0.658 | 2.93 12.65 13.6 
155 0.495 2.72 12.15 75.7 
423 0.440 | 25005 11.60 529.0 
423 0.330 | 1.87 11.60 529.0 
596 0. 384 1.97 11.43 995.0 
| 722 0.220 1.64 11.30 | 1438.0 
| Hexadecane 0.940 18 1.66 | 4.26 24. 0.51 
| | 28 1.10 3.14 24.15 1.19 
| 28 0.99 |} 2.92 24.1 as 
58 0.772 | 2.48 | 23.3 4.86 | 
Diameter, Pressure, Flame Distance 
Fuel P, (D,-D)/2, D, D Gr 
cm mm Hg cm eo a 
Hexzadecane 0.940 101 0.551 2.05 22.55 14.2 
187 0.440 1.615 21.6 46.3 
273 0.330 1.60 21.05 98.6 
273 0.274 1.485 21.05 98.6 
538 0.330 1.60 20.0 356.0 
1.208 17 1.73 4.65 31.9 0.956 
30 1.29 3.78 30.9 2. 
33 1.20 3.61 30.8 3.53 
33 1.21 3.63 30.8 3.53 
33 1.27 3.75 30.8 3.53 
33 1.13 3.47 30.8 3.53 
33 1.16 3.53 30.8 3.53 
62 0.848 2.90 29.9 11.6 
62 0.824 2.85 29.9 11.6 
103 0.658 2.53 28.9 31.6 
220 0.467 2.14 27.8 136.0 
411 0.302 1.815 26.4 455.0 
604 0.274 1.755 25.5 958.0 
604 0.274 1.755 25.5 958.0 
Xylene 1.208 122 0.490 2.19 32.8 63.4 
18) 0.431 2.07 32.2 135.0 
431 0.346 1.71 31.2 707.0 
760 0. 287 1.78 30.4 2060.0 
Heptane 1.208 760 0.195 1.60 25.8 2380.0 
Diameter Pressure, Flame Distance, 
Fuel D, P, (D,-D)/2, D ea Gr 
cm mm Hg cm 
Suspended Drops 
Ethanol (95. 5%) 0.12 760 0.0707 0.26 2.50 2.55 
0.092 0.0719 0,25 1.92 1.15 
0.092 0.0779 0.25 1.92 1.15 
2, 2'-Oxydiethanol 0. 106 760 0.0707 0.25 0.991 1.055 
0.0636 0.0637 0.19 0.595 0.228 
Pentane 12 760 0.0920 0.30 4.04 2.80 
0.0708 0.0779 0.23 2.39 0.575 
Hexane 0.113 760 0.0920 0.30 3.53 1.92 
0.0708 0.0920 0.26 2.21 0.473 
Heptane 0.12 760 0.0779 0.28 3.56 2.32 
0.092 0.0920 0.28 2.73 1.05 
Hexadecane 0. 106 760 0.0779 0. 26 2.18 0.956 
0.092 0.0849 0.26 1.90 625 


SPHERE DIAMETER: 0.94 cm 
FUEL: 2,2°- OXYDIETHANOL 
OXIDIZER: AIR 

PRESSURE OF AIR: 26 mmHg 
AIR FLOW RATE = 6100 cm¥min 
m = 5.82 x 10°°g/sec 


SPHERE DIAMETER: 0.94 cm 
FUEL: 2,2’-OXYDIE THANOL 
OXIDIZER: AIR 

PRESSURE OF AIR: 726 mmHg 
AIR FLOW RATE = 6100 cm¥Ymin 
m?9.4x10 g/sec 


Fig. 3 Flame surrounding porous sphere at approximately 


atmospheric pressure 
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relation between the dimensionless burning rate m/mo and the 
Reynolds number (5, 8). 

In the systems under study in this investigation, heat trans- 
fer by natural convection is considered to be the mechanism 
controlling combustion. For such systems Spalding has sug- 
gested that the burning rate is some function of the transfer 


Table 3 Steady-state mass burning rate as a function of diam- 
eter of porous sphere and concentration of oxygen in ambient 
gas; pressure = 1 atm, ambient temperature 20 C 


Sphere Buraia 
g Rete, 
Fuel —— Oxidizer m x 10° 
g/sec 
Ethanol (95.5%) 0.440 Air 3.22 
0.638 Air 5.49 
0.940 Air 8.86 
1.208 Air 13.1 
1.59 Air 19.3 
1.906 Air 27.1 
1.208 33 wgt % O, 14.6 
1.208 46 wgt % O, 16.4 
0.637 80 wgt % O, 8.9 
0.940 80 wgt % O, 13.8 
1.208 80 wgt % O, 19.5 
Butanol-1 0.440 Air 4.1 
0.638 5.83 
0.638 5.66 
0.940 9.8 
1.208 13.9 
1.59 21.6 
1.906 28.2 
Pentanol-1 0.440 4.64 
0.440 4.28 
0.638 6.01 
0.638 5.85 
Sphere 
Burning Rate, 
Fuel Diameter, Oxidizer m x 103, 
g/sec 
Pentanol-1! 0.940 Air 9.73 
1.208 14.1 
1.906 28.9 
1.590 22.1 
Octanol-1 0.440 3.83 
0.638 5.39 
0.940 9.20 
1.208 13.6 
1.590 20.4 
1.906 
Pentane 0.638 7 
0.940 14.2 
1.208 18.2 
Hexane 0.638 8.0 
0.940 12.8 
1.208 16.7 
Heptane 0.440 5.60 
0.440 5.75 
0.638 7.10 
0.638 7.25 
0.940 12.2 
1.208 16.5 
1.590 24.9 
Kerosene 0.44 4.16 
(Cy2.8 Hx) 0.638 7.28 
p = 0.81 0.940 9.6 
13.8 
20.7 
27.4 


number B, the Grashof number and the Prandtl number (3) 


m 
In the absence of mass transfer Ranz and Marshall (7) have 
indicated the applicability of an equation of the form 


Nu 
Nuo 
for heat transfer by natural convection (spheres). On the 
basis of these equations a relation of the following form is pro- 
posed and employed in the correlation of the data of this in- 
vestigation 
[3] 
Mo 
where a, M, N and P are empirical constants. The burning 
rate in the absence of convection mp is given by 
mo moCe 
— Pr = — = 2z In(1 4] 
The ratio m/mp is equivalent to Nu/Nup because the Nusselt 
number is proportional to m/D, (5). Since the Prandtl num- 
ber for gases is essentially constant (usually close to unity), 
Pr may be omitted in Equation [3]. 

It is apparent from Equation [3] that when free convection 
is reduced to zero (Gr = 0), the mass burning rate m becomes 
identical to the burning rate mo for a spherically symmetric 
flame. 

Burning rate determinations over a range of either pressures 
or sphere diameters are adequate to evaluate the exponent 1. 
In the present analysis this exponent was determined from the 
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Fig.4 Mass burning rate of ethanol (95.5 per cent) vs. diameters 
of porous spheres 


Sphere 
Burning Rate, 
Fuel — Oxidizer m x 103, 
g/sec 
Hexadecane 0.440 Air 3.41 
0.940 9.0 
1.208 7 
1.590 18.4 
1.906 23.6 
2, 2'-Oxydiethanol 0.312 2.46 
0.440 3.58 
1.590 20.4 
1.906 26.6 
Xylene 1.208 23.3 
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Fig. 5 Mass burning rate data. Porous spheres at atmos- 
pheric pressure. Oxidizer: air (unless otherwise noted) 
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Fig. 6 Correlation of mass burning rate data. Porous spheres 
at atmospheric pressure. Oxidizer: air (unless otherwise noted) 
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Fig. 7 Correlation of mass burning rate data. Porous spheres 
at subatmospheric pressures. Oxidizer: air 
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burning rate data of a number of different fuels at atmospheric 
pressures, for a series of porous sphere diameters. For the 
individual fuels, graphical presentation of the parameter 
(m/mo — 1) against Gr on log paper (mp is calculated by means 
of Equation [4], Table 4) gave parallel straight lines with a 
slope corresponding to N = 0.30. In Fig. 5 the linear relation- 
ship between m/D and Gr°-® for various fuels is presented. 
The straight lines drawn through the data intersect the or- 
dinate at the calculated value of mo/D based on Equation [4]. 
From the slopes of the lines of Fig. 5, numerical values for the 
function of B can be evaluated. It is thus found, empirically, 
that the value of M is equal to —0.44. Replotting the data 
obtained at atmospheric pressure (Fig. 6) yields a uniform 
correlation for all the fuels which can be represented by 


[(m/mo) — 1] B-44 = (5 


Since In(1 + B) is approximately equal to 0.87B°-*4 for the 
range of values of B in consideration, Equation [5] can be 
given alternately as 


[= In (1 + B) = [6] 
mo 

The burning rate data (precision of experiments +4 per 
cent) obtained with hexadecane, 2,2’-oxydiethanol, and xylene 
at reduced pressures (Table 1) are presented in Fig. 7 with the 
straight line based on Equation [5]. The experimental data 
summarized in Table 1 are plotted in Figs. 8, 9 also. These 
curves approach the theoretical curve for burning without 
convection (Equation [4], Table 5) as pressure is reduced. 
In the case of hexadecane (Fig. 9) some of the measured values 
of m fall below the theoretical curves. This overlap suggests 
that the method for evaluating mp is not completely satis- 
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Fig. 8 Mass burning rate as a function of pressure 
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Fig. 9 Mass burning rate as a function of pressure 
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Table 4 Numerical values of parameters and physicai properties employed in the calculation of m /D and Grashof number at 1 atm 


4 


FUEL AND OXIDIZER 
cal/g C 


g/sec cn 


, 
Gut: | 


cal/gC 


e 
— x 104, |— «10° 
c : D 


cal/g cal /g g/sec cm 


Ethanol (95. 5%)* 4.15 0.625 
Ethanol (95.5%)? 4.3 0.625 
Ethanol (95.5%)° 4.48 7 0.625 
Ethanol (95.5%)4 0.625 
Butanol-1* 0.555 
Pentanol-1° (0.682 
Octanol-1* 0.521 
Ethylene glycol" 0.575 
2, 2' -Oxydiethanol* 0.550 
Pentane* 2300 0.560 
Hexane* 2300 0.542 
Heptane* 2305 0.535 
506 
513 


2,2-4 Trimethyl pentane 0 

Hexadecane® 2310 0 

Xylene (m,p)* 2330 0.405 
0 


483 


, a 
Kerosene (Cro | 2320 


240 6,630 1. 4.82 
6,630 | 1. : 6.00 
6,630 | 1. ; 145 
6,630 | 1. ‘ 9.30 
8,051 | 2. 3:29 
8,420 | 2. : 5.36 
9,060 | 2. : 5.03 
4,290 | 1. : 3.85 
5,220 2. 4.80 

10,720 | 3. 7.75 

10,720 | 3. ; 1325 

10,720 | 3. 6.88 

10,672 | 3. 7.26 

10,600 | 3. : 4.88 
9,850 | 3. 6.45 

10,400 | 3. 5.62 


Adiabatic Flame Temperature (Reference 9) 


OXIDIZERS: a = air, Af = 0.232: b = 33 weight % 0,. x. = 


factory. Nevertheless, it is believed that the trend displaying 
a minimum at low pressures is valid. This point is reached 
when the conductive effects expressed by mo begin to control 
the process of heat transfer. From Equation [4] it is evident 
that mp increases with falling pressure because L, and hence B, 
decreases. Therefore the burning rate attains a minimum 
value and then increases with further decrease in pressure. 
The deviations of data measured with the porous spheres at 
low Grashof numbers, corresponding to small sphere diameters 
(Figs. 5, 6), point to a limitation in the experimental tech- 
nique. These deviations are accentuated progressively 
with increasing volatility or increasing B as observed in a 
similar series of tests on octane, heptane, hexane, pentane, and 
diethy] ether, all burning in air at atmospheric pressure. When 
the diameter is diminished, the heat flux (or burning rate per 
unit area) at the surface increases (see Equations [4,5]). It is 
believed that high volatility and high heat flux, combined 
with a factor such as incomplete wetting of the porous sur- 
face, account for this unexpected behavior. Because of this 
apparent limitation to the technique, the applicability of 
Equation [5] to the lighter hydrocarbon fuels could not be 
tested. 
In the calculation of mo, B, and Gr it was necessary to assign 
values to temperature levels and physical properties. For 
reasons of expediency, the flame temperature 7, was taken as 
the equilibrium adiabatic value (9) even though this practice 
is basically inconsistent with Equation [4] derived for con- 
ditions of nondissociation. It has also been shown that as 
a result of finite thickness of the reaction zone, the actual 
flame temperature is somewhat lower than the adiabatic 
value (3, 11). The terms \ and C were evaluated for air at a 
log mean temperature between 7, and T,. The use of the log 
mean instead of the arithmetic mean yielded lower (by 15 per 
cent) values of mp which were more in accord with the ob- 
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= 0.461; = 80 weight 7, = 0.736. 


served burning rates at low pressures (Figs. 8, 9). In the 
calculation of Gr, the properties « and p were similarly taken 
for air at the log mean temperature between 7’, and the tem- 
perature of the surroundings. However, since the heated, 
rising gases are bounded by the liquid surface assumed at the 
boiling point 7’, and the ambient gas at 7'., the temperature of 
the surrounding 7’, was taken simply as the arithmetic mean 
of 7, and 

Drop combustion studies above atmospheric pressure (1, 2) 
indicate a proportionality between burning rate and pressure 
raised to a fractional power 7, where 7 has values between 0.2 
and 0.5. Such dependencies, when considered to apply to 
narrow pressure ranges, are suggested by Equation [6] re- 
arranged in the following form 

m = mo + 0.17 mK P?-/In (1 + B) 
where P is the pressure and K replaces the pressure-indepen- 
dent terms of the Grashof number. Since mo/In(1 + B) is 
practically pressure independent, Equation [4], m becomes 
proportional to P®-” at high Grashof numbers, yielding an 
upper limiting value of 7 = 0.60. 

At constant pressure the distance of the flame at the for- 
ward stagnation point from the surface of a sphere or sus- 
pended drop (D, — D)/2 (Table 2) decreases with decreasing 
diameter, but the ratio of flame diameter to sphere or drop di- 
ameter increases. Similarly, for a given sphere or drop 
diameter the flame distance increases with decreasing pressure 
(Figs. 2,3). The data presented in Table 2 are useful in relat- 
ing the natural convection parameter Gr to flame distance. 
If one assumes that to a first approximation the mass burning 
rate of a drop is proportional to D.D/(D, — D) (Reference 12), 
Equation [6] may be rewritten 

in (1 + B) = 0.17Gr°-» 


| 
0.85 
ee 4 0.85 | | 
0.82 
avg.) | 
2 0.890 
0.785 | 
0.750 | 1 | 
0.580 
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0.504 
0.665 
Jet PROPULSION 


1 atm 


| 


Table 5 Numerical values of parameters and physical properties employed in calculation of m /) and Grashof number at sub- 
atmospheric pressures 


Ty ave L, a,/D, My * 104, Q, 
B 
mm Hg cal/g°C cal/g cal/g g/sec cm g/cm sec cm 3 
Hexadecane 15 432 0.58 67.0 148.0 4.55 6.22 4.34 0.429 
20 438 0.58 66.4 150.4 4.47 6.16 4.35 0.756 
(A/e, = 5.78x10°4 30 448 0.59 65.4 157.0 4.26 6.02 4.35 1.67 
g/sec cm) 50 462 0.59 64.4 164.4 4.05 5.87 4.36 4.46 
; 70 472 0.60 63.5 170.5 3.90 5.76 4.37 8.51 
100 483 0.60 62.6 176.6 ao 4e 5.64 4.38 16.9 
150 497 0.61 61.2 185.2 3.54 5.50 4.38 36.8 
200 506 0.62 60.2 192.2 3.39 5.36 4.39 64.0 
300 521 0.62 58.4 199.4 3.26 5.29 4.40 138.5 
500 542 0. 63 55.6 212.6 3.02 5.05 4.41 367.0 
760 561 0.641 52.8 224.6 2.83 4.88 4.42 803.0 
ee 
2,2'-oxydiethanol 15 390 0.60 134.0 197.0 3.94 Sto 4.31 0.486 
20 398 0.60 133.6 202.0 3.83 5.66 4.32 0.843 
A/c, = 5.75x10°4 30 408.5 0.60 132.4 207.0 33te 5.59 4.32 1.86 
g/sec cm) 50 422.6 0.61 130.4 215.0 Ke 5.46 4.33 4.97 
70 432.4 0.61 129:1 219.0 3.48 5.40 4.33 9.55 
100 443.3 0.62 127.6 225.0 3.38 ee S| 4.34 18.9 
150 456.3 0.62 125.7 232.0 3:27 5:23 4.35 41.0 
| 200 466.1 0.63 124.4 237.0 3.19 5.15 4.36 70.9 
300 480.5 0.63 122.4 245.0 3.07 5.05 4.37 154.0 
500 500.1 0.64 120.0 255.0 2.92 4.92 4.38 406.0 
| 760 Sit.3 0.65 118.0 264.0 2.81 4.80 4.39 886.0 
| ee 
| Xylene (m,p) 100 352.0 0.43 92.8 119.6 5.88 6.99 4.20 25.0 
| 122 357.0 0.43 92.2 121.3 5.78 6.95 4.20 36.6 
| (Ae, = 5.79% 10°* 181 368.0 0.43 90.8 124.9 5.60 6.84 4.22 78.0 
! g/sec cm) 400 392.0 0.44 87.0 132.0 5.24 6.64 4.25 354.0 
| 431 394.0 0.44 86.0 132.6 A 6.63 4.25 409.0 
760 413.0 0.46 83.6 139.1 4.92 6.45 4.28 1192.0 


OX(DIZER: Air T, = 293°K Y, = 0.232 (air) 


Since (Reference 4) 


Deo In (1 + B) 


In{ 1 + — 

i 

Equation [8] becomes 
Dea D. 

] = 10 
n(1 4 22) = 2 [10] 


Fig. 10 was prepared utilizing the parameters of Equation 
[8]. The resulting correlation obtained for different fuels and 
with suspended drops and porous spheres is rather encourag- 
ing although the numerical constants applying to the plotted 
straight line differ from those predicted by Equation [10]. 

It is evident from the data presented in Figs. 8, 9 and Table 
2 that although the burning rates observed at low pressures 
(15-100 mm Hg) approach rather closely the values calcu- 
lated for combustion without convection, the measured and 
calculated flame diameters, Equation [9], differ appreciably. 
It may be shown (13) that a possible cause for these large 
differences is the assumption, in the derivation of Equation 
[9], of the absence of chemical dissociation of combustion 
products and the consequent abnormally high flame tempera- 
tures. Since this assumption (also used in the derivation of 

quation [4]) has little influence on the heat flux at the sur- 
face, the mass burning rate calculated does not differ appre- 
ciably from that applicable to nonconvective systems with 
dissociation. However, the introduction of the adiabatic 
flame temperature markedly reduces the calculated flame dis- 
tance (12), and for some of the fuels examined the values ob- 
tained compare favorably with the flame distances measured. 
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‘ L,, heat of vaporization, was calculated using method of Reference 10 
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Fig. 10 Correlation of flame distance measurements 


Concluding Remarks 


A correlation represented by Equation [5] was obtained for 
a series of liquid organic hydroxy compounds and for several 
heavier hydrocarbons burning (in air and oxygen-nitrogen 
mixtures) on the surface of porous Alundum spheres of various 
sizes at atmospheric and subatmospheric pressures. This 
correlation demonstrates the contribution of natural convec- 
tion to the mass burning rate. In the limiting case at high 
pressures, the burning rate varies linearly with pressure 
raised to the 0.6 power (cf. Equation [7]), substantiating 
trends observed by other investigators in drop burning 
studies at higher pressure levels. As the pressure is reduced 
heat transfer by conduction becomes the controlling mecha- 
nism. 

Diameter measurements were made of the flames surround- 
ing both suspended drops at atmospheric pressure and 
wetted porous spheres at various pressure levels (atmospheric 
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and below). These data (Fig. 10) also demonstrate the in- 
fluence of free convection. The flame distances obtained 
with porous spheres and with suspended drops show con- 
formity in behavior. 

It was observed that the steady-state mass burning rate of 
a liquid sphere goes through a minimum as pressure is re- 
duced. At subatmospheric pressures the predicted pressure 
dependence for small drops (cf. Equation [5]) is so modest 
that for practical purposes the burning rate of hexadecane 
drops (for example) less than 1000 in diameter, may be con- 
sidered independent of pressure under free convection 
conditions. 
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Technical Notes 


Solar Sailing—aA Practical Method of 
Propulsion Within the Solar System 


RICHARD L. GARWIN! 


IBM Watson Scientific Laboratory, Columbia University, 
New York, N. Y. 


It is shown that commercially available metallized plas- 
tic film can be used as a solar radiation pressure sail for 
propulsion of space vehicles within the solar system. 
The method of propulsion is of negligible cost and is per- 
haps more powerful than many competing schemes. 


T IS difficult to exaggerate the importance of solar radiation 
pressure for the propulsion of satellites or space ships 
within the solar system; but since I have never seen any allu- 
sion to this powerful method, while less practical and more 
difficult schemes are frequently cited, I feel it desirable to 
publish this paper. 

The principle involved is simply to make use of the pressure 
of the sun’s light on a sail to propel a space ship as desired 
through the solar system. What is important is the area of 
sail per unit mass, and if, just as a practical example, we use a 
commercially available 0.1-mil-thick plastic sail equal in 
mass to the rest of the “space ship,”’ the mass per unit area 
will be 5 X 10-4 g/cm?. The sail may be aluminized without 
significant additional mass, in which case the sun’s radiation 
may be reflected, thereby doubling the maximum pressure 
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which may be exerted on the sail. The magnitude of this 
pressure on a normally oriented sail is P = 2W/c, where W is 
the energy incident per cm? per sec, and c is the velocity of 
light. At the earth’s orbit W = 2 cal/em?/min or 1.3 X 106 
erg/sec/cm?, so that P = 0.8 X 10-4 dyne/em?. The ac- 
celeration of the space ship is, therefore, a = p/o, where a is 
the mass/cm? and in the example is 5 X 1074 g/em?. There- 
fore a = 0.16 em/sec?. Thus on a space ship so constituted 
one may obtain an acceleration of as much as 1.6 X 10~‘ times 
that of gravity even with a gross 0.1-mil sail. 

The thrust may be varied from the maximum calculated, 
smoothly almost to zero by manipulation of the “shroud 
lines’; and the thrust may be directed almost without loss in 
magnitude to any angle within 45 deg of the radius vector 
from the sun, while smaller thrust (as cos 6) is available up to 
90 deg from the radius vector. 

Although the acceleration is numerically small, the ve- 
locity changes in reasonable times by significant amounts. 
Thus Av = 1.4 X 104 cm/see/day, so that the escape velocity 
from the earth (8 X 10° cm/sec) may be acquired in two 
months, and a velocity large enough to escape from the solar 
system entirely may be reached in a time of the order of a 
year. 

Before comparing the merits of solar sailing with those of 
ion propulsion, hydrogen expulsion by solar batteries, or 
nuclear or thermonuclear reactors, let us discuss two situa- 
tions in which solar sailing might well be used. To have a 
definite example, consider an earth satellite with reflective 
sails, each weighing 20 lb (10 kg). The sail will have an area 


2 X 104 
_ 2X = 4X 107 em? 
5 X 107-4 g/cm? 
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and will be more or less a circle of plastic film 70 meters 
in diameter. There are, of course, no requirements as to ab- 
sence of holes, tears, flaws, etc. The reflective coat plays a 
perhaps important role in eliminating electrostatic forces 
which might easily prevent opening of the sail. The sail may 
be attached to the ship by “shroud ribbons” ~200 meters 
long. Since the forces on the shroud ribbon will not exceed 2 
grams weight, the ribbons may be narrow strips of the same 
material as the sail. Let us for the moment grant that the 
sail may be expelled from its packing container after the 
“space ship” is in a satellite orbit at reasonable altitude. It 
is of interest to compute the time required for the sail to “fill 
with sunlight.” This is the time for the sail to go from a flat 
circle to a dished configuration of the order of 10 meters deep 


and is 
2 103 1/2 
ty = Cx ") = 80 sec 


Thus the sail may be furled and unfurled quickly compared 
with the 90-min period of its orbit around the earth. 

The simplest program for increasing the altitude of the 
satellite, at the same time increasing its energy and reducing 
its velocity about the earth, is to furl the sail by slacking 
half the shroud lines while the satellite approaches the sun 
and to unfurl the sail while receding from the sun. Clearly 
then, the mean projected acceleration along the velocity vec- 
tor isa! X 0.16 cm/sec? or 5 X 107? em/sec?, and the rate of 
decrease of orbital velocity is then 8 = 5 X 107? cm/sec’. 
Thus, an initial velocity of 8 X 10° cm/sec will decrease to 
zero in a time on the order of (8 X 10°5)/(5 XK 107?) = 1.6 X 
10’ sec, or 6 months. The orbit eccentricity grows during this 
time, and, depending on details, the altitude at perigee may 
remain constant or increase. 

It should be noted parenthetically that the torque exerted 
on a satellite by radiation pressure is ordinarily larger than 
that due to differential gravity effects. Indeed, if half of a 
20-in.-diam Vanguard sphere is painted black, while the other 
half is reflecting, and if the moment of inertia is 2 * 10° gm 
em?, the torque is of magnitude ~1 dyne-cm, and angular 
acceleration will be ~10~* radians per sec. The sphere will, 
therefore, revolve one radian in ~30 min due to this cause. 
For this same sphere with maximum gravitational quadrupole 
(two 5-kg masses at opposite ends of a diameter), the differen- 
tial gravity torque is 10 dyne-cm at 45 deg inclination. On 
this small-area device, the torque due to radiation pressure is 
not negligible. For less symmetrical devices (cones, ete.) a 
stable equilibrium may be reached pointing into the sun, 
while damping of long-period oscillations about this equi- 
librium may be achieved by slightly viscous fluids in some 
portion of the satellite. 

The energy required for operating the shroud lines is at 
most (2 < 10% dynes) X (104 em) or ~2 joules each 90 min, 
even if a nonconservative actuator is used. This is 20 
milliwatts, while the kinetic energy delivered to the satellite 
by the sail averages 60 watts. The maximum power falling 
on the sail is ~4 X 106 watts (4 megawatts) so that it is 
evident that the propulsion efficiency is not very high (~10~5 
or, in fact, v/c). A solar battery ~1cm? area supplies enough 
power to manipulate the shroud lines. 

As a second example let us consider a space ship of the 
same fairly high mass per unit area (0.1-mil sail of mass equal 
to that of the rest of the ship) ~5 X 10~4g/cm?. Now we 
want to go from the earth’s orbit to that of Venus and back 
again. 

The optimum sail tilt to go slowly with a small sail from 
one circular orbit to another is 6 = sin—! V/: 3/3, since this 
angle gives the maximum component along the velocity vector 
and thus allows the ship to lose or gain energy as rapidly as 
possible. In order to move in to Venus we must brake the 
ship to increase the orbital velocity (because the potential 
energy is twice the kinetic energy and negative). In the field 
of our sun with the space ship given above we have 


Marcu 1958 


SATELLITE 
WITH SAIL 
FURLED 


SATELLITE 


WITH SAIL 

UNFURLED 

Fig. 1 Use of the solar sail to increase the altitude of a satellite 
and to escape the earth 


SATELLITE WITH 
SAIL SET FOR 
APPROACHING 
THE SUN. 

Fig. 2 Cruising within the solar system by use of the sail 


(0.8) X (0.8 X 10~4) (: 2 2 
(5 X 10~) r 


where r/7 is the distance of the ship from the sun measured 
in earth’s orbit radii (~1.5 cm). But 


v = 2.4 X 106 (") 
r 


Thus we find for the rate of change of r with time 


1077 
r 


showing that the orbit is an equiangular spiral, opening or 
closing by a large factor each turn (indeed, by a factor ~20 
per year if the initial conditions were right for the motion to 
be an equiangular spiral in this noncentral-force field). 

The orbit equations are not difficult to solve in detail, but 
we may here be content with observing that the radiation 
pressure in our example is one-half the gravitational force be- 
tween sun and space ship, so that the space ship is entirely 
free to go “tacking’’ about in the solar system, arriving at 
any specified radius in times considerably less than a 
“vear’’ for the planet at that radius. For instance, one can go 
to Venus and back in less than an ordinary year with this 
first gross sail. 

Compared with escape from the earth, the efficiency of the 
sail rises for operation near the earth’s orbit, from ~10~° for 
earth satellite propulsion to ~10~‘ for solar navigating, since 
the pertinent velocity is higher and the duty cycle more favor- 
able. We may compare alternative schemes for propulsion, 
the one currently most popular being, perhaps, acceleration 
of the interplanetary ions (~10*/cem?) by an accelerator in the 
vehicle, powered by solar energy or nuclear or thermonuclear 
reactor. To achieve the 4 X 10* dynes thrust contemplated 
here with gentle use of a sail, one would have to accelerate to, 
say, 100 kev energy one ampere of interplanetary protons. 
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This would require only 100 kilowatts of power rather than 
the 4 megawatts of sunshine reflected (lower, of course, by the 
ratio of proton exhaust velocity to light velocity); but the 
space ship swallows up only (2.4 X 10°) X (10%) X (A) pro- 
tons per sec, where A is the frontal area of the accelerator. 
Thus A > 10° cm?, and the space ship must have a maw 
larger than the sail contemplated here. The use of a thermo- 
nuclear reactor to eject ~1 gram of hydrogen per day at 100 
kev (in our 20-kg ship) would be an eminently satisfactory 
solution if such reactors existed and were of zero mass. Even 
so, one must eject one’s own mass every 30 years to compete 
with the solar sail for cruising. 

It is obvious that there are considerable difficulties con- 
nected with space travel, but those connected with the sail 
appear relatively small. We note that a reasonable chemical 
fuel with exhaust velocity ~2 km/sec, requiring a ratio of 
take-off mass to final mass of e’/”, were V is the jet velocity 
and v the desired final rocket velocity, or a mass ratio of e4 = 
250 to reach 11 km/sec (escape velocity from the earth) would 
require to reach 24 km/sec a mass ratio of 1.6 X 10°, if rocket 
engines and tanks had zero mass. In actuality the mass 
ratios are very much larger. 


Possible Improvements 


It should be no great trick to make large areas of sail of 
thickness 2 X 10~* cm (as is used for beta-ray source back- 
ing), thus reaching a space ship area per gram twelve times 
smaller than that considered above. With this improved 
sail a satellite could escape the earth within a week, and could 
fall to Venus in less than a month and come back in a week 
(at high velocity, albeit). 


Combustion Instability in Solid 
Propellant Rocket Motors 


E. W. PRICE! and J. W. SOFFERIS? 
U. S. Naval Ordnance Test Station, China Lake, Calif. 


ECENT results of an experimental program have re- 
vealed some new and crucial information regarding the 
mechanism of excitation of combustion instability in solid 
propellant rocket motors, and the effect of oscillatory behavior 
on the quasi-steady burning rate of the solid propellant. The 
nature of this information and the tests which revealed them 
are described in the following. 

In conventional rocket motors the structure and orientation 
of the flow oscillations that correspond to acoustical modes of 
vibration of the combustion chamber cavity are very complex 
and not necessarily stationary with time (1,2). Thus it is 
difficult to know where to place a pressure gage for meaningful 
measurements of the oscillation; difficult to know what aspect 
of the burning and what region of the charge surface are re- 
sponsible for, or associated with, the excitation process; and 
difficult to assess the effects of the different aspects of oscilla- 
tory behavior (e.g., velocity oscillation and pressure oscilla- 
tion) on a partially burned charge. 

To circumvent these difficulties, an experimental burner 
was designed to favor oscillations in a stationary mode, and 
studies were made of the conditions favorable for combustion 
instability. A burner having a combustion chamber with 
circular cylindrical geometry was chosen in order to obtain 
longitudinal oscillations of very simple geometry, with fre- 
quencies comparable to those of the transverse oscillations in 
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Fig.1 No.1 research burner for studies of combustion instability 
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Fig. 2 Envelopes of standing pressure and velocity waves in the 
longitudinal mode of the burner cavity 


much larger motors. The burner is sketched in Fig. 1. In 
order to minimize the attenuation of oscillations in the longi- 
tudinal mode, the nozzle vent was placed midway between 
the ends of the burner, at a point of zero pressure oscillation 
of the fundamental mode of longitudinal vibration of the 
cavity (Fig. 2). The ends of the motor were made as flat as 
practicable and the propellant charge was made a very close 
fit in the burner to avoid dissipation of oscillation energy in the 
annular space between burner and inhibited outer surface of 
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the charge. In this design (Fig. 1) the pressure gage is always 
located at a point of maximum pressure oscillation (antinode) 
on the standing longitudinal pressure wave. The blowout 
disk at the opposite end of the burner provides protection for 
the equipment and leads to partial burns of the charges when 
the (normally progressive) pressure reaches blowout pressure. 

The experimental burner described has been used in a num- 
ber of tests designed to verify or disprove (for one propellant 
so far) the following contentions of various writers regarding 
combustion instability: 

1 Combustion oscillations are excited by periodic fluctua- 
tions in mass burning rate induced by pressure oscillations (as 
opposed to velocity oscillations or periodic fluctuations in 
turbulence in the region of velocity oscillations) (1,2,3,4). 

2 The most readily excited mode of oscillation of a circular 
cylindrical cavity is a traveling tangential mode (2). 

3 The flow at the downstream end of the charge (where 
channel discontinuities occur) generates noise which in turn 
excites combustion oscillations (5). 

4 Quasi-steady burning rate may or may not be increased 
by the action of the oscillatory flow and combustion (6). 

5 Changes in burning rate during unstable performance are 
favored by conditions of high quasi-steady stream velocity in 
the conduit of the burner (7). 

The conclusions regarding these issues and the supporting 
evidence were as follows: 

1 Flow oscillations were excited by action of pressure os- 
cillations on combustion but not by velocity oscillations. To 
determine this, tests were run on charges made up partly of 
dummy (cellulose acetate) propellant as shown in Fig. 3. The 
charges with dummy propellant in the middle would presum- 
ably oscillate much like full charges in the first mode if pres- 
sure (or density or temperature) oscillations were of primary 
importance to the excitation process, since the burning con- 
tinues to be present in the principal region of oscillation of 
pressure. 

On the other hand, the charges with propellant in the middle 
section would be expected to oscillate like full charges if 
velocity oscillation were primarily responsible for the excita- 
tion of oscillations. The pressure records of two tests on these 
charges are shown in Fig. 4, with a pressure record from a full 
charge for comparison. The records show that only relatively 
minor oscillations occurred in the test with propellant in the 
vicinity of the velocity oscillations, while the oscillations with 
the charge having propellant at the ends of the burner where 
pressure oscillation predominates were almost as severe as 
those obtained with the full propellant charge. 

Probably equally significant is the fact that the charge with 
propellant in the central section ‘‘chose”’ to oscillate in a mode 
that placed pressure antinodes in the region of the propellant 
material. From the observed frequency, the mode was found 
to be the third mode of the longitudinal mode (about 8700 
cps), the structure of which is sketched relative to the charge 
configuration in Fig. 2c. It is presumed that the second mode 
was not ‘‘chosen’’ because that mode has a pressure antinode 
adjoining the nozzle, and acoustic pressure tends to “leak 
out”’ through the nozzle in that mode. 

2 Contrary to predictions of Grad (1) and Cheng (3), the 
mode of oscillation most readily excited in the tests reported 
here was the longitudinal one. However, the time lag concept 
used by both these authors could permit this result in any 
particular design because the geometry change during burn- 
ing might not lead to any tangential frequencies compatible 
with the time lag involved. It was found that use of a longi- 
tudinal resonance rod (which would tend to suppress any 
transverse oscillations) did not affect the oscillations ma- 
terially. The changes in amplitude of the oscillations during 
burning, although to some extent reproducible, were not 
found to be associated with any multiple relationship between 
calculated frequency of the longitudinal modes and frequency 
of transverse modes (thus there apparently wasn’t much 
coupling between modes). The instrumentation used was 
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Fig. 3 Partial charges to localize the region of excitation 


PROPELLANT IN MIDDLE 


Fig. 4 Longitudinal pressure oscillations in the motor shown in 
Fig. 1 using the charge configuration shown in Fig. 3; double 
base propellant, 80 F 


linear only to 15 ke, so that no transverse oscillations would 
have been recorded if they had been present. However, the 
above comments support the belief that no transverse oscilla- 
tions were present. 

3 Test results indicated that the intensity of fluid dy- 
namically excited noise was not a conspicuous factor in ex- 
citing combustion instability. Flow velocities in the combus- 
tion chamber were very low in the tests reported as compared 
to conventional motors, being in the vicinity of Mach 0.06 at 
the midpoint of the burning period when the majority of in- 
stabilities developed. The configuration did not have channel 
discontinuities comparable to conventional motors where 
vortex shedding processes would readily excite fluid dynamic 
noise. 

4 For the double-base propellant used in these tests, the 
burning rate during periods of oscillatory behavior was ob- 
served to decrease by as much as 50 per cent at lower pressures 
(500 to 1700 lb/in.*), and to increase more than 50 per cent at 
higher pressures. The first effects are illustrated in Fig. 5, 
where an observed pressure-time curve is compared with a 
theoretical one. The pressure level is severely depressed dur- 
ing the period of oscillation, rises some about two-thirds of the 
way through burning when the oscillations temporarily de- 
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Fig. 5 Comparison of theoretical pressure-time curve with one 
obtained on test with combustion instability; double base pro- 
pellant, 80 F 


crease in amplitude, and then becomes further depressed be- 
low normal as the oscillations grow again. 

5 Results of tests with full charges indicate that high 
mean stream velocity is not essential to the modification of 
quasi-steady burning rate accompanying flow-combustion 
oscillations. The motor used in these experiments had an 
effective ratio of throat area to port area of 0.21 at the start 
of burning, 0.05 at the end of burning. This ratio should re- 
sult in Mach numbers of 0.12 and 0.03 at the start and end of 
the burn, respectively. These values are very low compared 
to most rockets (which have initial ratios of throat area to 
port area of 0.4 to 0.8). The modification of burning rate 
does not seem to follow the change in mean gas velocity, but 
rather primarily the amplitude of the pressure oscillations. 
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A Preliminary Study of 
Vortex-Boiling Burnout Heat Fluxes 


W. R. GAMBILL! and N. D. GREENE? 
Oak Ridge National Laboratory,? Oak Ridge, Tenn. 


Nomenclature 
Ai = internal surface area of test section, ft? 
Di = internal diameter of test-section tube, in. 
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L_ = test section heated length, in. 
Py = power required for fluid flow across test section, hp 
Pq = power corresponding to heat absorption rate g, hp 
AP = test section pressure drop, psi 
q = rate of heat transfer, Btu/hr 
V = axial fluid velocity, fps 
W = water weight flow rate, lb/hr 
x = steam quality at test section exit, per cent by weight 
pi = liquid density evaluated at exit temperature, lb/ft* 
Subscripts 

burnout condition 


il 


internal 


ECENT experiences with source-vortex flow of fluids 
(air and water) in forced-convection heating without 
boiling (to be published) led the authors to believe that this 
type of flow field might be especially advantageous in a boiling 
system. By “source-vortex flow’ is meant fluid flow with 
both axial and tangential velocity components; more specifi- 
cally, it refers in the present instance to the flow of a swirling 
or whirling fluid from one end of a tube to the other. With 
source-vortex flow, steam bubbles formed at the inner surface 
of the heated tube wall are surrounded by a rapidly rotating 
body of liquid of much higher density, and the gravitational 
acceleration arising from the vortex motion would be ex- 
pected to cause a much stronger inward radial displacement 
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Table 1 Results of Investigation 
Inlet Test 
Axial Degree of static Heated exit 
Voltage velocity, subcooling pressure, length, Mate- = quality, Py o (q/ Ai), 
Run variation fps atexit,°F psig ID, in. OD, in. in. L/D; rial % Pq] Btu/hr-ft? 
1 Step 12.1 36 93 0.275 0.300 3.063 11.14 Copper 0.0 0.26 5,940,000 
taps 

2 Variacs 11.8 0 94 0.275 0.300 3.187 11.60 Copper 1.24 0.19 8,240,000 

3 Variacs 12.35 0 99 0.275 inlet—0.2875 2.875 10.44 Copper 1.49 0.20 9,690,000 
exit—0.300 

4 Variacs 13.61 72 100 0.250 0.300 1.20 4.80 Inconel 0.0 0.43 10,980,000 

5 Variacs 13.18 135 248 0.300 0.375 3.05 10.16 Copper 0.0 0.37 8,630,000 

6 Variacs 16.3 150 248 0.250 inlet—0.270 2.25 9.00 Inconel 0.0 0.39 10,810,000 
exit—0. 281 

7 Variacs 24.0 49 285 0.250 inlet—0.270 1.625 6.50 Inconel 0.0 1.01 17,250,000 
exit—0. 278 


o: bubbles away from the inner wall toward the tube center 
than in the case of linear liquid motion. Were this effect 
sufficiently strong, the inner surface would be freer of bubbles 
at a given heat flux than in standard straight-through flow, 
thereby permitting the attainment of a greater peak flux 
lefore burnout. These considerations led the authors to 
carry out the study described herein. The large measured 
heat fluxes empirically demonstrate the qualitative validity of 
these early concepts. 

A schematic representation of the experimental system is 
given as Fig. 1. The test-section inlet, shown as Fig. 2, 
consisted of a spiral-ramp fluid introduction zone followed 
by a convergent cone section which accelerates the rotation 
rate of the water through reduction of radius. This same 
vortex generator, with orifice diameter of 0.104 in., was used 
with all of the test section tubes. The available flow rates, in 
conjunction with the design of the vortex generator, permitted 
the attainment of several thousand gravities (calculated) at 
the test section inlet. The resistance-heated, thin-walled 
tubes (of copper or Inconel) were cooled by the vortex water 
flow within them. 

The tubes for runs 3, 6 and 7 were tapered, possessing 
thinner walls at the inlet than at the exit. The use of this 
geometry was motivated by an attempt to more efficiently 
utilize the greater circulation in the freshly formed vortex at 
the test section inlet by generating more power in this zone, 
and reducing the power generation with length to at least 
cualitatively match the decreasing vorticity of the water 
flow. This sort of variation of heat-generation rate with 
length is similar to that obtainable by selective placement of 
fuel in a nuclear reactor. It might be noted here that in all 
tests about four tube diameters separated the inlet orifice 
and the initial heated portion of the test section, so that some 
vortex decay had taken place before heating began. 

The degree of subcooling listed is the minimum for the sys- 
tem (at exit conditions) and is equal to the saturation tem- 
perature of water at the exit static pressure (always just 
above atmospheric, except in run 6) minus the actual exit 
water temperature. The inlet static pressure was measured 
at the small orifice of the vortex generator at the tube inlet. 

The results of this investigation are summarized in Table 14 
for the seven experiments believed to be of highest accuracy 
(estimated maximum over-all error in burnout heat flux of 
+10 per cent, probable error of 3 to 4 per cent). 

The quantity designated (P;/P,) per cent represents the 
ratio of calculated power required to force the water through 
the test section to the rate of heat absorption by the water 
at the burnout point, expressed in the same power units. 
This was calculated from 


_WaP_ 
18,750 pz 
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(total 

The flow in all tests exhibited high stability right to the 
burnout point, and the system could apparently be operated 
indefinitely a few per cent below the burnout condition. 

The heat fluxes reported in this paper are extremely large 
for the conditions of operation; in fact, at least five times 
greater than any values previously reported for comparable 
system conditions. The largest experimenta! values of 
burnout heat flux that have been reported in the literature, 
to the knowledge of the authors, for flow through straight, 
round tubes, are those of Rohsenow and Clark (1)5 (3.1 X 
10° Btu/hr-ft? for water at 2000 psia, 30 fps, 240 F average 
subcooling, and L/D = 52), and of Buchberg et al. (2) (3.8 X 
10° Btu/hr-ft? for water at 507 F exit temperature, 44 fps, 
236 F exit subcooling, and L/D = 110). This latter value 
corresponds to about 5.3 X 10° Btu/hr-ft? for L/D = 7, 
according to the approximate rule that halving L/D in- 
creases (q/A;), by ~10 per cent (3). The highest value for 
any geometry or experimental conditions is apparently that 
of Gunther (4), who obtained 11.4 < 10® Btu/hr-ft? for water 
at 114 psia exit pressure, 40 fps, and 256 F subcooling at exit. 
Gunther obtained his values with a thin (0.004 in.) elec- 
trically heated metal strip in a rectangular channel of L/D, = 
6; the strip was cooled from both sides. Since Gunther’s 
values are consistently among the highest reported in the 
literature, a comparison of our run 4 and his test 5, for which 
conditions are almost identical, is revealing (see Table 2). 

Were all of Gunther’s experimental conditions identical 
with those of our run 4, it is estimated that his (¢/A)» would 
be 2.3 X 10° Btu/hr-ft?, which gives a 377 per cent higher 
peak heat flux for the vortex case. When comparisons are 
made with other literature data, the advantage becomes more 
nearly 500 per cent for the vortex case. 

It might be noted that the burnout heat flux for run 7 
corresponds to an average volumetric heat release rate within 


* Added in proof: Burnout heat fluxes for runs 2 to 7 of 
Table 1 (Variac control) have been empirically correlated, with a 
4.8% average error, by 


(q/Ai)s = (15,750 V? + 8,070,000) (1.18 — 0.03 L/Di) 


For the bulk exit subcooling range of 0 to 150 F, there was only 
a very small dependence of (¢/A;), on degree of subcooling. 

Recent preliminary comparisons on an over-all power expendi- 
ture basis indicate that linear and vortex flow patterns will give 
approximately equal burnout fluxes at equal flow power in the 
low flow power range, but that vortex boiling exhibits an increas- 
ing superiority with increase of flow input power. 

5 Numbers in parentheses indicate References at end of paper. 
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Table 2 
Gambill-Greene Gunther (4) 
Run 4 Test 5 

Axial water velocity, fps 13.6 12.3 
Pexit, psia 15.0 14.4 
Subcooling at exit, °F 72 71 

L/De 4.8 6° 

10.98 10® 2.0 x 10° 


(q/A)o, Btu/hr-ft? 


® Based on hydraulic diameter of } in. given by Gunther (4). 


the tube wal! of about 173 kw/cc, based on water heat ab- 
sorption alone. 

The ascending magnitude of burnout heat flux with run 
number (Table 1) is primarily attributable to improvements in 
experimental technique as more was learned of the nature of 
the process of burnout in this new flow system. Run 2, 
e.g., gives a (q/A)» greater than that in run 1 because a 
transition was made from the use of step transformer taps for 
voltage control to Variacs, by which power level could be 
increased gradually and smoothly rather than by finite incre- 
ments only. The result is still higher in run 3 because of 
slower Variac adjustment as a result of having learned the ap- 
proximate burnout point from run 2. Another factor tend- 
ing to further increase burnout heat flux in run 3 was the use 
of the tapered-wall test section, which provided an area for 
electrical current flow at inlet about half that at the exit. 
The calculated heat flux at the inlet for the burnout condition 
is 14.8 X 10° Btu/hr-ft? for this run. The further increase 
in (g/A)» of run 4 is primarily associated with the reduction 
of L/D from ~10 to ~5. 

In run 5, the tube (thick-walled copper) probably failed 
because of a combination of undesirable mechanical stress 
conditions that existed. Runs 6 and 7 serve as the basis of 
one important conclusion that has been reached concerning 
the relative importance of the pertinent variables. In run 6, 
the exit pressure was held at a constant 100 psig but the burn- 
out flux showed little change, except as expected from dif- 
ferences in velocity and L/D, from that of run 4. In run 7, 
exit pressure was reduced to atmospheric, a maximum inlet 
pressure was attained, and the effect of the greater velocity 
showed up in a much greater burnout flux. This comparison 
indicates that the velocity effect is more important than any 
pressure level effects. This is to be expected, since, in the 
case of a vortex, increased circulation or vorticity also brings 
about greater static pressure of the fluid at the wall, greater 
local degree of subcooling, and possibly greater radial bubble 
migration. 

Heat flux at the inlet (zero length) of the tapered tube of run 
7 is computed to be 20.9 X 10° Btu/hr-ft? at the burnout con- 
dition. 

Burnout in run 7 was prematurely triggered by a small tem- 
porary decrease in water flow. The water circuit used in all 
tests was open and subject to small fluctuations in supply 
pressure. 

The large heat fluxes obtained for the case of net steam 
generation, runs 2 and 3 of Table 1, suggest another possible 
advantage of vortex boiling—that any gas released or vapor 
generated will be transported to the center of the tube where 
it will do the least harm. This action occurs to some extent 
in linear flow because of Bernoulli forces, but many investiga- 
tors have nevertheless concluded that greater amounts of dis- 
solved gas will decrease peak flux in subcooled boiling with 
linear flow. 

Some idea of the excessive velocities required in standard 
axial flow systems for attainment of these high heat fluxes 
can be obtained from the following comparison. Both 
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Buchberg (2) and Gunther (4) indicate that for constant sub- 
cooling and pressure in a given tube 


By a correlation proposed by Gunther (4) the (¢/A)» to be 
expected for axial flow under the conditions of our run 7 is 
1.68 X 10° Btu/hr-ft?. The vortex value of 17.25 xX 108 
is 10.27 times greater and, by Equation [3], a linear velocity 
of (24) (10.27)? = 2540 fps would be required to attain 
17.25 X 10° for an exit subcooling of 49 F. 

In practice, of course, this velocity would be neither avail- 
able nor necessary, since the high pressure required would at 
the same time increase the subcooling and thus decrease the 
required velocity. Nevertheless, the comparison does illus- 
trate the variation reyuired in the velocity variable alone. 

A new heat-transfer mechanism apparently is acting in 
vortex boiling, and the high fluxes obtained are by no means 
due solely to higher velocities and degrees of subcooling at the 
wall. The burnout fluxes obtained are much larger than 
values calculated for linear flow even when the orifice velocit, 
(maximum tangential in the system) is used in available 
burnout equations. All comparisons and analyses point to 
the existence of some mechanism such as postulated previously. 

The burnout values presented here can quite likely be in- 
creased to considerably higher levels by the use of greater 
velocities, pressures and degrees of subcooling, all of which 
decrease the maximum attainable bubble size and increase 
the burnout heat flux. None of the conditions of the present 
study (except Z/D) approach the optimum for maximum 
(g/A)», which further emphasizes the immense heat fluxes 
possibly attainable. 

An effort will be made to extend this study of heat-transfer 
effects within a vortex flow field to other ranges of experi- 
mental variables. It is felt that a (¢/A)» value of 40,000,000 
Btu/hr-ft? is not unattainable, and that further refinements 
of experimental technique may be expected to reduce the 
maximum error for (q/A)» to below +5 per cent. Further 
studies will involve determination of more than (q/A), and 
the additional data will permit an attempt to formulate a 
generalized correlation. 

Further results will be published as they are produced. 
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Condensed phase temperature distributions, experi- 
mentally determined for several burning liquid sys- 
tems, have been compared with distributions predicted for 
yarious modifications of the thermal combustion model. 
Reasonably good agreement exists between the theory and 
experiment, and the results have been analyzed with re- 
spect to the relative importance of various types of heat 
production. In general, the temperature distribution can 
be largely accounted for by thermal conduction alone. 
The effects of radiation absorption and chemical reaction 
are discussed. 


Introduction 


REVIOUS papers (1,2)* in this series have dealt with the 
measurements of the combustion wave temperature dis- 
tribution in burning liquid systems and the application of the 
method to a fundamental study of the mechanism of combus- 
tion. This paper describes a portion of these studies in which 
the experimental condensed phase temperature profiles are 
compared with those expected for a thermal model of the 
combustion process, i.e., one in which the flow of heat is the 
process of primary importance. 

According to the thermal theory, the temperature at any 
point in the condensed phase of a burning system is deter- 
mined by three possible processes: The conduction of heat, 
the absorption of radiation from hot gases above the surface, 
the generation (or absorption) of heat by chemical reaction. 

The desired procedure would be to describe the situation in 
mathematical terms and then to inspect the resulting expres- 
sions for the distribution of temperature. Unfortunately, the 
presence of the chemical reaction term gives rise to a nonlinear 
differential equation, which must be solved to obtain the tem- 
perature distribution. Because of this nonlinearity, one can- 
not obtain an analytical solution giving the temperature as an 
explicit function of distance. Numerical solutions of this 
equation have been obtained by Olds and Shook (3), and an 
analysis of their results indicates that the chemical reaction 
term influences the over-all temperature distribution only in 
the region extremely close to the burning surface. For ex- 
ample, a numerical solution containing a chemical term of 
13.6 X exp (—35400/RT’) cal em~ sec, a surface tem- 
perature of 450 C, an initial temperature of 20 C, a thermal 
diffusivity of 6.9 X 10-4 cm? sec~! and a burning rate of 1.0 
cm sec! showed that the effect of the chemical term was felt 
only at temperature above 420 C. 

The major portion of the temperature-distance curve (cov- 
ering the region from 20 to 420 C) could, however, be fit quite 
accurately by solutions to the above heat flow problem ob- 
tained by considering only heat conduction and radiation ab- 
sorption. The nature of this “conduction-radiation”’ solution 
has been described by Reid (4). It thus appears that useful 
information about the radiation and conduction parameters 
may be obtained by comparing the experimental temperature- 
distance curves with the ‘“conduction-radiation” solutions. 
The method used in making this comparison follows. 


The Curve-Fitting Method 


The experimental profile data were obtained by having a 
liquid burn past a thermocouple, as described in (1). From 
the observed liquid consumption rate and the thermocouple 
readings, one has values for the temperature at different dis- 
tances from the burning surface. Attempts were made to fit 
the data so obtained by means of functions of the sort 


[2] 
From (4) it is seen that 
I I 
A=(7,-T. +—.—— B = — ———_- 
( pe(hta pe(h?a — v) 
v 
B = h? C= Y=a 


* Numbers in parentheses indicate References at end of paper. 


Marcu 1958 


©00 

100 

50} 

° 
5 10 20 25 50 
DISTANCE, ABITRARY. UNITS 
Fig. 1 ‘‘Best fit’? to the ethyl nitrate data at 7.7 atm 


where 7) is the maximum temperature, 7’, is the initial tem- 
perature, J is the radiation intensity at the burning surface, a 
is the radiation absorption coefficient, h? is the thermal dif- 
fusivity, v is the liquid consumption rate, and p and c are 
density and specific heat, respectively. 

The temperature distribution is given by Equation [1] when 
heat conduction and radiation absorption are the contributing 
factors, and by Equation [2] when there is heat conduction 
alone. 

The procedure used in attempting to fit the data by the 
above functions was to consider first a series of points well 
below the surface and obtain a best fit (in the least square 
sense) by means of Equation [2]. The curve fitting was then 
extended to progressively higher temperatures so as to in- 
clude points closer and closer to the surface. For each set of 
points a curve was obtained, and the sum of the squares of 
the deviations of the points from the curve was calculated. 
Upon dividing this value by N — 1 (where N is the number of 
points) and taking the square root, one obtains the rms (root 
mean square) error. This value was used to measure how 
well the equation fit the data. An excellent fit could usually 
be obtained when just points well away from the burning 
surface were used. In addition, good fits could sometimes be 
obtained almost up to the burning surface. The results ob- 
tained with ethyl] nitrate at lower pressures, for instance, were 
of this type. In some cases the fit to the data by means of 
Equation [2] became unsatisfactory when too many points 
were included. In these cases it was found that the use of 
Equation [1] made it possible to obtain a better fit over a 
wider range of points. 

From the mathematical standpoint, it is not easy to obtain 
a best fit of a set of points by means of a sum of two ex- 
ponentials. Several procedures were tried before a satisfac- 
tory one was found. The method used was as follows: First 
trial values for 8 and y were taken, and with them “‘best’’ 
values for A, B and C were determined by the method of 
least squares. Values for 7 were then computed at each 
point by means of Equation [1] and compared with the 
experimental data. Again, the rms error was used as a meas- 
ure of how well the curve fit the data. The trial values of 6 
and y were then systematically varied until a minimum value 
of the rms error was obtained. That is, the 8 and y were such 
that if either was increased or decreased, the rms error would 
increase. The curve obtained with these values for 6 and y 
was considered to be the ‘“‘best”’ fit in the least square sense. 
Fortunately a high speed electronic computer (IBM 701) was 
available for the calculations. A sample fit is shown in Fig. 1. 
The circles are experimental points and the heavy line repre- 
sents the calculated curve. 
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Table'1 Curve-fitting the ethyl nitrate data 


Cons. Equa- rms 
Pressure, rate, tion h? X 104, error, 
atmosphere cm sec~! fitted 8, cm~! °C 


4.3 0.028 [2] 34.4 8.1 2.3 
“if 0.052 [2] 73.1 7.1 2.2 
Be 0.083 [2] 140 5.9 2.2 
14.5 0.122 [1] 164 7.4 1.0 


Results and Discussion 


Ethyl nitrate data. The extent to which the ethy] nitrate 
condensed phase profiles obtained at various combustion pres- 
sures can be fitted by Equation [1 or 2] is shown in Table 1. 
Shown from left to right are the combustion pressure, the 
liquid consumption rate, the equation to which the data were 
fitted, the value of the parameter 8 obtained from the curve- 
fitting, the thermal diffusivity calculated from 6 and the rms 
error. In deciding whether to fit the data with Equation 
[1 or 2], the curve-fitting was considered adequate whenever 
the rms error was less than 3 deg. This estimate was based 
on a consideration of the reproducibility of a given set of ex- 
perimental results. 

It can be seen that up to a combustion pressure of 11.1 atm, 
the measured temperature distributions can be satisfactorily 
fitted by the single exponential or “conduction only” solution. 
At all four pressures, essentially the entire liquid phase profile 
can be fitted. Although it appears to be necessary to bring 
in the radiation term at the highest pressure, an analysis of the 
data shows the rms error to be extremely sensitive to small 
changes in the conduction parameter, but relatively insensi- 
tive to changes in the radiation parameter. Thus the tem- 
perature distributions, even at 14.5 atm, can be accounted for 
primarily by thermal conduction, with radiation absorption 
being a second order effect, at best. For this reason, no sig- 
nificance is attached to the calculated values of the radiation 
parameters, although they appear to be of the correct order of 
magnitude. 

A further check on the validity of the curve-fitting method 
can be made by comparing the predicted thermal diffusivity 
values, as in Table 1, with values obtained by independent 
measurement. An analysis of twelve sets of profile data ob- 
tained at various combustion pressures gave values for h? that 
agreed to within 15 per cent and for which an average value 
of 7.1 X 10-4 em? sec! was obtained. Combining this value 
with an average specific heat (5) of 0.45 cal g-! deg and an 
average density (5) of 1.05 g cm~, one obtains a predicted 
value of 3.4 X 10~4 cal em sec! deg for the thermal con- 
ductivity of liquid ethy] nitrate. The average value of the 
thermal conductivity over the temperature range involved, 
as obtained by direct measurement (6) is 4.7 X 10~4 cal em! 


Table 2 Curve-fitting the metriol trinitrate data 
rms error, 
To, °C Equation [1], °C 
187 1.8 
208 0 
226 2.2 
242 2.8 
256 3.7 
264 4.6 
271 5.8 


sec! deg~!, Considering the error involved in the tempera- 
ture profile work, the agreement between the predicted and 
observed thermal conductivities is felt to be satisfactory and 
would indicate that the experimental results are compatible 
with the thermal theory. 

Metriol trinitrate data. In treating the temperature data 
obtained on the metriol trinitrate‘ (82 per cent)-triacetin® (18 
per cent) system at a combustion pressure of 34.9 atm, only 
the portion of the profile below about 240 deg can be fitted 
satisfactorily. The surface temperature is about 290 dev. 
Table 2 shows the way in which the rms error begins to diverge 
when the data are fitted to Equation [1] and the curve-fitting 
is extended toward the surface. Given in the table are tlie 
maximum temperature attained in the curve-fitting and the 
rms error. When the data were fitted to the ‘conduction 
only” solution, rms errors of 3.0 and 4.6 deg resulted at 7) 
values of 187 and 208 deg. Even so, an analysis of the vari:- 
tion in the rms error with small changes in 6 and y indicates, 
as before, that the shape of the profile is determined largely 
by thermal conduction and that radiation absorption is, «t 
most, a second order effect. 

The data of Table 2 are in qualitative agreement with the 
results of a series of thermal stability measurements made on 
the metriol trinitrate mixture. For these observations, the 
liquid sample was pressurized to 35 atm in the combustion 
chamber and temperature-time curves were recorded as the 
liquid was heated at varying rates with stirring. In all cases, 
rapid self-heating due to exothermic chemical reaction oc- 
curred at temperatures of 190-210 deg. Thus the failure o! 
Equation [1] to fit the profile data above 240 deg can be ac- 
counted for on the basis of interference from chemical reac- 
tion. 

Here also a comparison can be made between the liquid 
thermal conductivity value predicted from the curve-fitting 
and that obtained by direct measurement. An analysis of four 
separate sets of profile data obtained at 34.9 atm gave values 
for the thermal diffusivity agreeing to within 8 per cent, the 
average being 6.7 X 1074 em? sec~!. Combining this value 
with an average specific heat (5) of 0.43 cal g~! deg and an 
average density (5) of 1.20 g em~*, one obtains a predicted 
thermal conductivity of 3.5 107-4 cal 
The average value of the liquid thermal conductivity over the 
temperature range of interest, as obtained by direct measure- 
ment (5), is 3.2 X cal deg Again, the 
comparison indicates that the experimental data are in reason- 
able agreement with the thermal theory. 
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Technical Comments 


Terminology in Rocket Combustion 
Instability 


E. W. PRICE! 
U. S. Naval Ordnance Test Station, China Lake, Calif. 


'EXHE problem of oscillatory performance of rocket motors 

has been studied and discussed to some extent during the 
past 15 years by every agency in the rocket business. The 
isolated nature of the various studies has lead to a diversity of 
terminology which often results in misinterpretation of the 
reported information, e.g. (1,2).2. This is further complicated 
by the fact that much of the terminology was established in 
the early days when knowledge of the phenomenon was little 
more than speculation. Confusion continues to be injected 
into the issue by the fact that the majority of measurements 
are made by instruments which are unresponsive to the rapid 
oscillations which occur during combustion instability. In 
view of this situation, some clarification of terminology seems 
inorder. The remarks below on this subject pertain primarily 
to solid propellant rockets, although the situation is much the 
same with liquid propellant rockets. 

The phenomenon of combustion instability with solid pro- 
pellants is advantageously viewed as consisting of two major 
aspects. The first is the self-excited oscillation of the gas 
flow-combustion system. The second aspect is the ‘cumula- 
tive effect”’ of this oscillatory phenomenon on the mean burn- 
ing rate of the solid propellant (averaged over several cycles 
of oscillation). The first aspect is referred to as ‘‘sonance”’ by 
Green (3), and is referred to as “dynamic instability” by 
Cheng (4). 

The terminology ‘“‘resonance’’ was originally applied to the 
combustion instability phenomenon at a time when actual 
oscillations had not been successfully measured. Thus the 
term ‘resonance’ was introduced because of a speculation 
that oscillatory behavior was present, but was in practice 
associated primarily with very different, low-frequency 
measurements properly involving the second aspect of com- 
bustion instability noted above. It may end up that this 
trend toward application of the term “resonance” to the 
mean burning rate aspect of the problem was fortunate, since 
it is not accurate when applied to the self-excited oscillations,* 
and since recent work by Green (3) suggests that the mean 
burning rate effect may actually be a true resonance phe- 
nomenon between the oscillating gas phase system (forcing 
function) and the “natural frequency”’ of the solid phase com- 
bustion process determined by the reciprocal of the time lag 
of the solid phase process. At present the terms “resonance” 
and “resonant burning’ are widely and indiscriminately 
applied to all aspects of combustion instability. 

The terminology “combustion instability” or “unstable 
burning” has become increasingly evident. Its meaning is 
sufficiently vague and general to avoid actual incorrect inter- 
pretation. However, it implies that it is the combustion that 
is unstable, whereas, it is the interaction between the com- 
bustion and the flow oscillations that is actually unstable. By 
selective reference to the combustion, this terminology is mis- 
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3 Resonance is a concept which applies to systems with forced 
oscillations rather than self-excited ones (5). 
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leading, and blocks the way to subsequent terminology dis- 
tinguishing the high speed and the low speed aspects of the 
phenomenon mentioned above. 

In the case of solid propellant rocket motors, the oscillatory 
phenomena occur at very high frequencies (pressure oscilla- 
tions higher than 50 ke have been observed), and the pressure 
sensing and recording equipment used in routine testing is 
usually unresponsive to these high frequencies. As a result, 
only the effect (if there is any) due to the change in mean 
burning rate is detected. This often leads to failure to detect 
combustion instability, or erroneous identification of other 
performance anomalies with combustion instability. The re- 
sulting looseness in meaning and interpretation further con- 
fuses the meaning of combustion instability. In view of this 
situation it is undoubtedly wise to use increased discretion in 
interpretations of test results and wording in the description 
of results and interpretations. This writer does not consider 
any of the terminologies currently in use to be accurately 
descriptive of the phenomenon and at the same time unam- 
biguous with respect to the aspect of the phenomenon. It is 
suggested that terminology be sought which covers the three 
concepts indicated in the diagram below. 


] 
General title for the ‘combustion instability’’ or | 
‘Tesonance burning’’ phenomenon—all aspects | 


| 
| | 


Stability or instability of the 
combustion-flow oscillation 
interaction 


Modification of burning rate 
(averaged over several cycles 
of oscillation) attendant to in- 
stability 
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A Remark on the Free Vibration of a 
Variable Mass 


GEORGE LEITMANN'! 


University of California, Berkeley, Calif. 


Nomenclature 


mass of body (tank) 

initial mass of body 

position vector of center of mass of body relative to an 
inertial reference system 

position vector of center of mass of body relative to a 
reference system fixed in the invariant part of the 
body (body system) 


— 
o 
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R, = position vector of “effective center’’ of exit surface rela- 
tive to body system 

1 = R, — 

® = angular velocity of body system 

F = resultant of external forces 

Ve = exit velocity (of out-flowing mass) relative to body 
system 

t = time 

d/dt = differentiation with respect to time in inertial system 

0/dt = differentiation with respect to time in body system 

x = linear displacement of mass center of body 


N A recent note in this Journal, Miesse (1)? discussed the 
solution of a differential equation which was to describe 
the free vibration of a variable mass, a tank in which the mass 
of liquid decreases linearly with time. However, the basic 
equation, Equation [1] of (1), is incorrect for two reasons: 
1 The linear momentum of a body of variable mass, i.e., of 
a body which loses mass across a portion of its bounding sur- 
face is not merely M(dr,/dt), (2 and 3), but is given by 


2 Newton’s Second Law, which was originally stated for 
particles, cannot be applied directly to a variable mass such as 
the one under study. A correct application of Newton’s 
Second Law, as used in (2 and 8), yields the equation of 
motion 
d*r, | oR, dM 


= 14+ —-2— — 
dt at | dt 


The second of these reasons is probably the more important 
one. 

If it is assumed that the “effective center” of the exit sur- 
face and the center of mass of the body are stationary in the 
body system 

aR, _ _ 
If, in addition, the body system does not rotate relative to the 
inertial system 


= 0 
Then for linear mass flow 
dr, dM 
M— =F + V, 2 
dt? dt 2] 


Equation [2] also results under the assumption that the linear 
momentum of the body is approximated by Mdr,/di, (4). If 
now 

dM 


<F 
at 
then for the one-dimensional case 
dx 
I(t) de F (xz + (3 


which'differs from Equation [1] of (1) by the term (dM/dt)- 
(dx/dt). Equations [8 and 6] of (1) then become, respec- 
tively 


vf” — uf’ + (v2? — p? — 2p)f = 0........ [5] 
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On the Free Vibration of a Variable 
Mass 


J. J. MARCINEK! 


Convair Division of General Dynamics Corp., 
Pomona, Calif. 


N OCTOBER 1957 Jer Propunsron, C. C. Miesse (1)? pre- 
sented work on the free vibration of a variable mass. In 
establishing the force balance for the system, Newton's 
Second Law was used incorrectly. The Second Law, as com- 
monly stated, applies only to mass conservative systems. 

In a variable mass system momentum changes are the resuit 
of (a) velocity changes caused by the impressed force and (}) 
mass changes that add or remove momentum of the system. 
Thus: 


time rate of (system momentum change + momentum 
entering or leaving system) = force 


A formal explanation and proof can be obtained in (2). 
Using Miesse’s nomenclature throughout, the force equa- 
tion for a vibrating system of decreasing mass is 


dx d t\ dx dz 
—=—— 4h —'h@—...... 1 
where 
Mu _ Mo 
so that 
t\ dx 
Substituting 
K t 


into Equation [2] results in 
for which the solution is 


which differs from Miesse’s solution by one order. 
The frequency approximation, Equation [20] of (1), is then 


written as 
1 4(p + 1)? 1 


More important, the amplitude variation for large values of 
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v is (Equation [22] of (1)) 


2 


Therefore the amplitude will decrease with time for all 
positive values of p. The amplitude can only increase if p 
is a velocity-dependent forcing function with a value less than 


—}. 


then Equation [9] becomes 
+ of’ + [v2 — (p — = 0.......... {12} 


The approximate frequency is then given as 


1 41 — p?-1 


The amplitude variation is obtained as 


For the case of the system to which mass is added, Equation 


us 

#+e ae: + eran eee eee 9] Since v increases with time in this case, the amplitude of the 
oscillations will decrease for values of p > 3, and increase for 
If now p< }. 
K t 

v = 2G = (: (10) References 

’ 1 Miesse, C. C., “The Free Vibration of a Variable Mass,” 
and JET PRopuLsioN, vol. 27, Oct. 1957, p. 1103. 

2 Rosser, J. R., Newton, R. R., Gross, G. L., “Mathematical 
[11] Theory of Rocket Flight,’’ McGraw-Hill, New York, 1947, p. 6. 
George F. McLaughlin, Contributor 
Tracking system (2,817,081). ; W. Van B. a When the rudder is out of its neutral posi- 
toberts, Princeton, N. J., assignor to the sf | —- . tion, the pressure actuates linkage con- 
U.S.A. ; ; — nected to a tab at the rudder trailing edge. 
System for measuring the speed and F Self-supporting shell for flying body 
course of a rocket. A continuous wave (2,817,484). F. M. B. Stenzel, Alamo- 
transmitter is tuned to resonate with the ht. gordo, N. Mex., assignor to Machine Tool 

length of the rocket. Recorders located , Works Oerlikon. 
below the expected trajectory are syn- st 4 Hollow metal wire helically coiled to 


chronized with the firing of the rocket to 
be observed. 


Variable-area inlet constrictor for ramjet 
engine (2,817,209). C. W. Besserer, Sil- 
ver Spring, Md., assignor to the U. S. 
Navy. 

A bellows connected to a pitot-tube 
presented forwardly adjacent to the in- 
take. Expansion of the bellows, in re- 
sponse to velocity and atmospheric den- 
sity, shifts a central body forward within 
the intake. 


Airplane design (181,858). A. Kartvelli, 
Huntington, N. Y., assignor to Republic 
Aviation. 

Fighter bomber known as the F-105 
Thunderchief, driven by a P&W J-75 
turbojet engine of 15,000-lb-thrust class. 
It features swept forward air intake ducts 
to shatter shock waves. 


Non-tip off launcher (2,818,779). C. J. 
Koeper, Hunstville, Ala., assignor to the 
U.S. Army. 

Steel section with parallel upper and 
lower plates carrying a flange forming a 
double rail receiving notches in the rocket 
fin section to guide the missile during ini- 
tial flight. After the rocket has traveled 
beyond the end of the double rail to a sin- 
gle rail extension, a collar encircles the 
forward end of the rocket. 

Firing pin for rocket fuze (2,817,245). 
H. B. Wylie, Barrington, Vt., assignor to 
the U.S. Navy. 

An intermediate thin-walled section 
crushable upon application of a sudden 
axial force. A spring-loaded shutter in a 
compartment containing an explosive ele- 
ment is prevented from movement to an 
armed position until the inertial force has 
subsided to a predetermined level. 

Aircraft control force modifier (2,817,483). 
T. C. Hill, Harford County, Md., assignor 
to The Martin Co. 

Ram air led from an intake port at the 
leading edge of a vertical fin surface to 
pressure responsive means within the fin. 


form a body shell. Longitudinal stringers 
contact the shell; an agent bonds adja- 
cent turns of the wire to one another and 
to the stringers. 

Jet powered adjustable pitch helicopter 
rotor (2,818,122). J. G. Saco, Palm- 
dale, Calif. 

Fluid lines extending radially outward 
from a center post to a pair of rotor blades. 
A valve in the post selectively and simul- 
taneously opens the lines to fluid passage 
in either or: both blades for controlling 
their thrust and pitch. 


Directing apparatus (2,818,782). J. R. 
Jasse, Paris, France, assignor to Societe 
Nouvelle des Etablissements Brandt. 

System of launching self-propelled pro- 
jectiles. An anemometic device driven 
by electric motors receives corrections in 
elevation and direction, one of which is 
dependent on the wind component in the 
plane of fire and the other on the transverse 
wind component. 


Jet propulsion of helicopters (2,818,223). 


Eprtors Note: Patents listed above were selected from the Official Gazette of the U.S. Patent Office. 


Printed copies of Patents may 


be obtained from the Commissioner of Patents, Washington 25, D. C., at a cost of 25 cents each; design patents, 10 cents. 
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Aerodynamics & Propulsion 


Information 
manual about 
APL and its 
programs 
now available 


The Applied Physics Labora- 
tory (APL) of The Johns 
Hopkins University is unique 
in that we are neither an in- 
dustrial nor an academic 
organization, but rather a 
composite, having drawn 
freely from the methodolo- 
gies of each. 

For thirteen years APL has 
pioneered in guided missiles. 
Today we are engaged in a 
broad pro: of R & D for 
the Navy; in addition, we are 
responsible for technical di- 
rection of industrial and 
academic contractors in de- 
veloping the Terrier, Talos 
and other major weapons and 
weapons systems. Our staff 
members enjoy not only the 
stimulus of association with 
their immediate colleagues at 
APL, but also with those in 
other organizations of con- 
siderable stature. 


NEW 30-PAGE PUBLICATION 


A few positions for senior engi- 
neers and scientists are now 
open. Information on our ac- 
complishments and goals is 
available in a new 30-page pub- 
lication, just off the press. 

In it staff leaders representing 
each of the various disciplines 
and fields outline the nature of 
their programs. Information on 
our new laboratory in Howard 
County, Md. (equidistant be- 
tween Baltimore and Washing- 
ton) is also included, together 
with facts on the outstanding 
communities in which our staff 
members live. 


Quantity 1s somewbat limited. 
May we suggest you send now to: 
Professional Staff Appointments, 


The Johns Hopkins University 
Applied Physics Laboratory 


8617 Georgia Avenue, Silver Spring, Md. 


F. L. Doblhoff, Zell am See, Austria, as- 
signor to The Fairey Aviation Co., Ltd. 
(England). 

Turbine gases directed through at least 

one rotor blade propulsion nozzle to cause 
blade rotation. Gases are also directed 
against the rudder to achieve orientation 
about a vertical axis during hovering. 
Jet engine fuel, pressure ratio and nozzle 
area control (2,818,703). I. W. Victor, 
Cincinnati, Ohio, assignor to General Elec- 
tric Co. 

Means responsive to the ratio of com- 
pressor discharge pressure to turbine ex- 
haust pressure to control the variable area 
discharge section. A ratio is maintained 
at a first predetermined value when the 
shutoff valve is open, and at a second pre- 
determined value when the valve is closed. 
Turbine engine starting system (2,- 
818,704). H. C. Karcher, Mansfield, 
Ohio, assignor to General Motors Corp. 

Control valve with a radial array of 
pivotal overlapping shutter flaps at the in- 
let of the combustion chamber. The 
valve closes the flaps upon ignition of the 
combustible charge. 

Rotary jet propelled engine (2,818,705). 
L. L. Ledet, Charleston Heights, S. C. 

Rotor including a jet exhaust within a 

frame, and mounted in a housing. Louver 


openings in the housing permit air under 
ressure to enter the combustion chamber. 
eated air escaping to the atmosphere 
creates thrust, imparting motion to the 
rotor. 


Firing Chamber 


2,818,808 


W. S. 


Jet perforating gun (2,818,808). 
Dill, Midland, Tex. 

Gun for perforating casing in oil wells. 

Series of shaped charges are stacked in a 
carriage so each directs a jet outwardly 
through the casing at an angle of 120 deg 
from the jets of adjacent shaped charges. 
A firing Primacord has parallel burning 
paths of travel, firing charges in regular 
unidirectional sequence. 
Valve arrangement for afterburner igniter 
control (2,819,587). R. J. Coar, Hart- 
ford, Conn., assignor to United Aircraft 
Corp. 

Two pistons displace fuel into the com- 
bustion chamber. The smaller piston is 
moved by the larger one to supply addi- 
tional fuel through nozzles into the after- 
burner. 

Ventilated suit refrigeration unit (2,819,- 
590). F. H. Green, Palos Verdes Estates, 
Calif., assignor to The Garrett Corp. 

Heat exchanger for cooling hot com- 
pressed air supplied to an airtight en- 
closure. An expansion turbine further 


cools and reduces the pressure of a portion 
of the cooled air. Temperature respon- 
sive means controls a valve in response to 
the temperature of air flowing into the 
enclosure. 


= ZA 


2,818,914 


WS 


System of containers for at least two 
liquids and a pressure gas (2,818,914). 
H. Thomann and A. Steck, Vaduz, Bartle- 
grosch, Liechtenstein, assignors to Ma- 
chine Tool Works Oerlikon. 

Central tube with two lids slidably in- 

serted from each end to form a pressure 
gas chamber. Two liquid containers sur- 
round the gas chamber and one another. 
Pistons in the containers put gas into con- 
trollable communication with the liquid 
containers. 
Variable area turbine entrance nozzle 
(2,819,732). R. A. Faetz, University 
Heights, Ohio, assignor to Thompson 
Products, Ine. 

Movable fluid directing vanes between 
a pair of radially spaced coaxial rings. 
A third ring operates to controllably move 
pivotal vanes, regulating flow between the 
movable and fixed vanes. ane 
Turbojet engine fuel and nozzle control 
system (2,820,340). J. Dolza, O. P. 
Prahar, and J. B. Wheatley, Indianapolis, 
Ind., assignors to General Motors Corp. 

Means for varying the fuel rate coor- 
dinately with the effect of the ambient 
air condition on fuel requirements. 
Braking and reverse turbine for gas tur- 
bine engines (2,820,341). C. A. Amann, 
Birmingham, Mich., assignor to General 
Motors Corp. 

Exhaust duct receiving motive fluid 
from a forward turbine and a radial inflow 
reverse turbine. A valve controls the 
flow of fluid from the turbines to the duct. 
Air-jacketed annular combustion chamber 
for a jet-propulsion engine (2,821,066). 
J. &. Clarke, S. R. Jackson and G. J. 
Hudson, Burnley, England, assignors to 
Joseph Lucas Industries Ltd. 

Walls with corrugations shaped and 
arranged to form inner and outer bound- 
aries of combustion zones in the combus- 
tion chamber. Walls are positioned so air 
from the jackets can enter the combustion 
zones and set up a circular motion around 
the axis of each burner nozzle. 


2,820,397 


77770) 


Airborne missile dispenser (2,820,397). 
T. H. Durkin, Roslyn, Pa. 

Rotatable container from which missiles 
may be released from an aircraft by cen- 
trifugal force in a sequential order, or 
simultaneously over a wide area beneath 
the line of flight. 
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For records you can read at once... 
new Kodak Linagraph Direct Print Paper 


This paper is designed for use when you must know 
what’s happening while it’s happening . . . in such ap- 
plications as rocketry, control, computing, product 
design, component testing and nuclear research. It re- 
quires no processing, comes out of the recording in- 
strument ready to read. 

To get such instant monitoring of high-speed varia- 
bles you eliminate the powerful energy amplification 
factor of photographic development. Traces must be 
put down by a high-intensity u-v light source. 

Kodak Linagraph Direct Print Paper gives you a 


EASTMAN KODAK COMPANY 
Graphic Reproduction Division 
Rochester 4, N. Y. 
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wide range of recording speeds—from 0 to 3,000 cps. 
Records have adequate permanence for extended stor- 
age. If used under sunlight or fluorescent lamps they 
should be stabilized. The surface of the paper takes 
notations in pencil or ink. Records are suitable for 
use in recording system readers. 

Several major manufacturers now offer recording 
instruments based on the characteristics of Kodak 
Linagraph Direct Print Paper. We'll gladly send you 
their names on request. If you want more detailed in- 
formation on the paper, we'll send that too. 
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Book Reviews 


Heat Transfer—Volume II, by Max 
Jakob, Wiley, New York, 1957, 652 
+ xxxii pp. $15. 

Reviewed by J. Ep>warp SUNDERLAND 

Purdue University 

The second volume of Jakob’s treatise 
on heat transfer completes the most com- 
prehensive survey of heat transfer in 
existence. Dr. Jakob was particularly 
well qualified for writing a treatise on heat 
transfer as he was intimately familiar 
not only with American research in heat 
transfer, but also with German research. 
He has made use of many German refer- 
ences and has presented the most signifi- 
cant German contributions to the field of 
heat transfer. He has done a remarkable 
job in giving credit to those who have 
made significant contributions in heat 
transfer. Although the book was essen- 
tially complete before Jakob died in Jan- 
uary 1955, considerable work was neces- 
sary in order to prepare it for press. This 
work was carried out by one of Jakob’s 
colleagues, S. P. Kezios. 

As in Volume I, the author uses a partly 
historical approach in order to show the 
development of the ideas leading to our 
present state of knowledge. The book 
does not lend itself to the easy solution of 
problems by merely substituting into 
convenient formulas, but its particular 


Ali Bulent Cambel, Northwestern University, Associate Editor 


asset comes from a clear and concise 
presentation of theoretical, mathematical 
and experimental analyses. 

The book is divided into three parts: 
Heat Radiation in Spaces of Simple Con- 
figuration, Selected Fields of Application, 
and Supplements to Volume I. The 
first part, dealing with radiation, continues 
the material on radiation which is covered 
in Volume I. In the first volume some 
basic laws and properties of radiation are 
given, but the treatment is limited to 
cases where surfaces are separated by 
nonabsorbing mediums. In Volume II 
the treatment is extended to cover various 
boundary configurations and the effects 
of absorbing mediums. 

In the second division of Volume II 
Jakob selected certain specialized topics 
of heat transfer applications. The topics 
are: Applications of thermometry, heat 
exchangers, regenerators, cooling towers, 
falling liquid films, cooling of surfaces ex- 
posed to hot gases, heat transfer through 
laminar and turbulent boundary layers 
at high fluid velocities, heat transfer in 
liquid metals, and steady state heat trans- 
fer in packed columns. In the chapter on 
thermometry the author discusses meth- 
ods for the accurate determination of tem- 
perature and heat flow measurements. 
Following the chapter on ordinary heat 
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exchangers (recuperators), in which a 
general theoretical analysis is made of 
unidirectional and counterflow heat ex- 
changers, there is a most. complete tre:t- 
ment of regenerator theory. Due to the 
extreme complexities encountered in re- 
generator theory, the author discusses 
only some of the relatively simple analyti- 
cal solutions, but treats with more detuil 
the numerical and graphical solutions. 
An analysis is given for the working con- 
ditions of cooling towers. The theoretical 
and experimental results for both laminar 
and turbulent film cooling are presented in 
the chapter on cooling by falling liquid 
films. Asa result of the increasing impor- 
tance of cooling surfaces which are exposed 
to hot gases, the author includes a chapter 
covering the following topics: The effect 
of radiation from furnace gases to boiler 
tubes, the cooling of hollow gas turbine 
blades by a fluid in forced or free convec- 
tion inside the blade, and transpiration 
cooling. Boundary layer theory is given 
for the case of high velocity flow with 
laminar boundary layers; for instance, a 
method is given for determining the tem- 
perature distribution in a fluid flowing 
parallel to a plane plate with a uniform or 
an arbitrary surface temperature. The 
theoretical methods and experimental 
results are presented for recovery factors 
for both laminar and turbulent boundary 
layers. An excellent survey is given for 
heat transfer in liquid metals which in- 
cludes a detailed presentation of Lyon’s 
modification of Martinelli’s analogy, Jen- 
kins’ and Deissler’s modified mixing length 
theories, a discussion of contact resistance 
between solid surfaces and liquid metals, 
and a comparison of theoretical and ex- 
perimental data. In the chapter on packed 
columns, heat transfer is studied between a 
fluid and a bed of packed solids and the 
container wall. The fluid flow and heat 
transfer are assumed to be under steady 
state conditions. 

In order to incorporate in Volume I 
many of the recent papers associated with 
heat transfer, a short supplement to 
Volume I is included as the last major 
division of Volume II. It must be noted, 
however, that the most recent references 
quoted in the Supplement as well as the 
rest of the book have been written in 
1953 and 1954. Although there has been 
much work done since the book was writ- 
ten, it will continue for many years to be a 
most outstanding reference for those en- 
gaged in the field of heat transfer. 


Spectroscopy at Radio and Microwave 
Frequencies, by D. J. E. Ingram, Philo- 
sophical Library, New York, 1956, 
332 pp. $15. 

Reviewed by S. F. SIncER 
University of Maryland 


Since the war the development of micro- 
wave techniques for radars has stimulated 
research in a new region of spectroscopy. 
It turns out that the microwave region is 
as useful as the optical region for many 


JET PROPULSION 


: 
ce | 
350°F 
AIR B j 
MENTS, 
SON 


DOUGEAS ENGINEERS AND 
DOUGLAS 
te tre:t- 2 
e to the 4 eee your 
re detiil 
sora can be as big as 
oretical 
laminar your talents 
ented in 4 
Now in its 17th year, the Douglas missiles 
| is program is projected far into the future 
chapter by such exciting new projects as THOR 
e effect 
> boiler Out of such veteran projects as Nike and 
turbine . 3 y Honest John are coming fantastic new mis- 
sile systems to challenge the finest engineer- 
s given ing talents in the land. 
_ Since early in World War II, Douglas has 
epi : ae been engaged in missile projects of prime 
flowing importance. New engineering teams are con- 
orm or stantly being formed for research, design, 
The BS Cee development and production. Engineers 
advance rapidly as Douglas expands its 
indary — leadership in this challenging field. 
a You are stimulated to accelerate your 
career by the importance of each assignment 
-, Jen- ... by the help of your associates who are 
length — recognized experts in missile work... by the 
™ vastness of opportunity for engineers in this 
company that is run by engineers. 
acked oe There is no more promising future than 
the that which awaits you in the Douglas 
| heat Missiles Divisions. 
teady THOR —7n intermediate range ballistics mis- 
me I sile now going into mass production—has top 
with : priority in our country’s program for national 
For complete information, write: 
ences E. C. KALIHER, 
s the MISSILES ENGINEERING 
been eo DOUGLAS AIRCRAFT COMPANY, 
wave 
hilo- 
956, 
icro- 
ated 
is 
any 
Marcu 1958 


, 
4 
‘ 
. 
4 
ys 
: 
: 
203 
F 


Spheres of titanium alloy represent 
the optimum configuration:and 
material of maximum strength/weight 
for the storage of gases and liquids 
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Rheem Aircraft, presently supplying 
pressure units to current ICBM pro- 
grams, is directing its years of experi- 
ence to the development and pro- 
duction of all types of vessels and 
propellant chambers. 


RHEEM MANUFACTURING CO. 
AIRCRAFT DIVISION 


11711 woodruff avenue, downey, california 


Engineers: Join Rheem in challenging technical work. 


204 


applications. Research workers, although 
nonspecialists in this field, will want to 
know about the various applications to 
their own problems. Microwave spectros- 
copy can be used for gases as well as for 
solid state work. One of the most im- 
portant applications, for example, is the 
detection and measurement of free radicals. 
Although not specifically discussed in the 
book, one of the important applications 
might be to combustion problems. 

The first portion of the book discusses 
techniques, the production and detection 
of microwaves and various wave guide 
techniques. For the person who buys a 
ready-made microwave spectroscopy setup 
these chapters are of lesser intportance; 
but in many cases one will wish to build his 
own setup. In any case the information 
presented is unusually detailed; many 
clear circuit diagrams are given and refer- 
ences to the literature are quite up-to-date. 

Next, microwave spectroscopy systems 
are discussed, together with information 
on resolution and sensitivity. Then there 
follow discussions on the theory and results 
of gaseous spectroscopy, paramagnetic 
resonance and ferromagnetic resonance. 
Next follows a discussion on radiofre- 
quency spectroscopy covering the region 
up to about 300 mc (10 electronvolts per 
quantum). Radiofrequency spectroscopy 
has been uséful for solid state problems. 

To give a view of the wide field applica- 
tions of these techniques, one can list de- 
termination of nuclear spins, magnetic 
moments and electrical quadrupole mo- 
ments of nuclei, methods for aligning nuclei 
by different techniques, and finally the 
wide field of measurement of various 
chemical parameters. 

Applications of a more practical nature 
include the construction of frequency 
standards by means of atomic clocks, the 
measurement of magnetic fields using, for 
example, nuclear precession methods, 
stabilization of magnetic fields by reso- 
nance methods, chemical analysis by vari- 
ous nondestructive methods, detection of 
free radicals, the study of irradiation 
damage and applications in solid state re- 
search, gas discharges and radio as- 
tronomy. 


Analysis of Feedback Control Systems, by 


Robert A. Bruns and Robert M. 
Saunders, McGraw-Hill, New York, 
1955, 383 pp. 


Scott H. CAMERON 
Armour Research Foundation 


This book is concerned more with 
the practical or applied aspects of 
control system design than with the 
theoretical. Also no prior knowledge of 
servomechanism theory is needed in order 
to deal easily with the material presented. 
In fact, the only background required is a 
knowledge of the steady state theory of 
a-c networks, the fundamental laws of 
dynamics and a familiarity with ele- 
mentary differential equations. 

The first ten chapters of this book 
treat the characteristics of the components 
of which control systems are comprised, 
since it is the author’s view that a real 
understanding of feedback control theory 
can be most readily obtained after one is 


already familiar with the operational 
characteristics of the components. This 
portion of the text is particularly good in 
that it treats not only the usual electrical 
and electromechanical components, but 
the hydraulic and pneumatic components 
of which so many of our modern systems 
are composed. The operation and con- 
struction of each of the components is 
explained and, where appropriate, the 
transfer function derived. The transfer 
function is treated throughout as a simple 
generalization of steady state a-c network 
theory rather than a ratio of Laplace 
transforms. 

The following chapters deal with various 
aspects of feedback theory; while thes» 
chapters are adequate, they certainly find 
their counterpart in any number of recent 
texts. The material which is presented 
is adequate in detail and each chapter i: 
concluded with a set of excellent problems. 
The two concluding chapters presen( 
rather sketchy discussions of the de- 
scribing function and phase plane methods 
of nonlinear analysis. This reviewer 
thought that the chapter on experimenta! 
techniques for the determination of trans- 
fer functions was particularly good. 

This book is apparently intended and, 
in the opinion of the reviewer, would serve 


very well as a text for an introductory 


course in servomechanism theory and 
practice. 


Books Received 


Electro-Technology, by M. G. Say, 
Philosophical Library, New York, 1955, 
167 pp. $6. 

Principles of Metal Casting, by Richard 
W. Heine, McGraw-Hill, New York, 
1955, 639 pp. $7.50. 

Magnetic Materials in the Electrical 
Industry, by P. R. Bardell, Philosophical 
Library, New York, 1955, 288 pp. $10. 

The Sun, by Giorgio Abetti, MacMillan, 
New York, 1957, 336 pp. $12. 

The Inexplicable Sky, by Arthur Con- 
stance, Citadel Press, New York, 1957, 288 
pp. $3.50. 

Introduction to Electrical Applied 
Physics, by N. F. Astbury, Philosophical 
Library, New York, 1957, 241 pp. $10. 

Introducing Astronomy, by J. B. 
Sidgwick, MacMillan, New York, 1957, 
259 pp. $3.50. 

Automation: Its Purpose & Future, by 
Magnus Pyke, Philosophical Library, 
New York, 1957, 191 pp. $10. 

Zirconium, by G. L. Miller, Butter- 
worths Scientific Publications, London, 
1954, 382 pp. 

Electro-Magnetic Machines, by R. 
Langlois-Berthelot, Philosophical 
brary, New York, 1955, 535 pp. $15. 

Atomic Energy Research at Harwell, by 
K. E. B. Jay, Philosophical Library, 
New York, 1955, 144 pp. $4.75. 

Electricity and Magnetism, by J. New- 
ton, Philosophical Library, New York, 
1956, 613 pp. $10. 

Radio Telemetry, by M. H.-Nichols and 
L. L. Rauch, Wiley, New York, 1956, 
461 pp. $12. 

Basic Physics, by Alexander Efron, 
John F. Rider Publisher, Inc., New York, 
1957, 724 pp. $8.95. 
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Technical Literature Digest 


Materials of Construction 


Studies of Structural Failure Due to 
Acoustic Loading, by Robert W. Hess, 
Robert W. Fralich and Harvey H. Hub- 
bard, NACA TWN 4050, July 1957, 11 pp. 

Lithium Emerges as Major Aircraft 
Material, by William Beller, Amer. 
Aviation, vol. 21, July 15, 1957, pp. 47- 
55. 

Effect of Stress on Creep at High 
Temperatures, by H. Laks, C. W. Wise- 
man, O. D. Sherby and J. E. Dorn, J. 
. ppl. Mech., vol. 24, June 1957, pp. 207- 

13. 

Cermets and Metallic Refractories, 
by Robert Steinitz, pp. F-6-7, Metal 
Strengths at High Temperatures (Table), 
pp. F-8-9, Ceramic Materials (Table), 
pp. F-8-9, Thermal Limits of Non-metals 
and Chemicals (Graphs), pp. F-10-11, 
Beryllium Properties (Table), p. F-12, 
Aviation Age, Res. and Dev. Tech. Hand- 
book, 1957-1958, Section F. Materials: 
Research and Development Trends, pp. 
(Metals-non-metals-chemicals). 

Corrosion and Ignition of Titanium 
Alloys in Fuming Nitric Acid, by John B. 
Rittenhouse, Nicholas D. Stolica, Stephen 
P. Vango, Julia S. Whittick and David 
M. Mason, Wright Air. Dev. Center, TR 
56-404, (ASTIA AD 118028), Feb. 1957, 
of pp. 

Tensile Properties of Inconel X Sheet 
Under Rapid-heating and Constant-tem- 
temperature Conditions, by Ivo M. Kurg, 
NACA TN 4065, Aug. 1957, 20 pp. 

Compressive Stress-Strain Properties 
of 17-7 PH and AM 350 Stainless-steel 
Sheet at Elevated Temperatures, by Bland 
A. Stein, NACA TN 4074, Aug. 1957, 
21 pp. 

Torquing of Bolts for Nuclear Reactors, 
by F. J. Mehringer, Atomic Energy 
Comm., KAPL-1645, Aug. 22, 1956, 124 
pp. 

Uranium Alloys for High-temperature 
Application, by Henry A. Saller, Ronald 
F. Dickerson and William E. Murr, 
Atomic Energy Comm., BMI-1098, June 
25, 1956, 45 pp. 

Optimum Stresses of Structural Ele- 
ments at Elevated Temperatures, by 
Arthur Schnitt, M. A. Brull and H. S. 
Wolko, Trans. ASME, vol. 79, July 1957, 
pp. 959-966. 


Instrumentation, 
Telemetering, Data 
Recording 


Electroanalogic Methods, 3, by T. J. 
Higgens, Appl. Mech. Revs., vol. 10, no. 
2, Feb. 1957, pp. 49-54. 

An Evaporated Gold Bolometer, by E. 
Archbold, J. Sci. Instruments, vol. 34, 
no. 6, June 1957, pp. 240-241. 

A Thermistor Probe, by W. Belfield and 
R. W. Johnson, J. Sci. Instruments, vol. 
34, no. 5, May 1957, pp. 209-210. 


Eprror’s Note: Contributions from E. R. 
G. Eckert, J. P. Hartnett, T. F. Irvine Jr. and 
P J. Schneider of the Heat Transfer Labora- 
tory, University of Minnesota, are gratefully 
acknowledged. 
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The James Forrestal Research Center, Princeton University 


Simple Equipment for a Small Low Tem- 
perature Laboratory, by R. F. Broom and 
A. C. Rose-Innes, J. Sci. Instruments, 
vol. 33, Nov. 1956, pp. 420-421. 

Measuring Instruments for Resistance 
Thermometers, by Jaroslav Edl, Elektro- 
techniky, Obzor, vol. 45, Dec. 1956, pp. 
609-616 (in Czech). 

Temperature and Pressure Measure- 
ments in an Overcompression Cloud 
Chamber, by D. W. Hadley, J. C. Bower 
and J. Walker, J. Sci. Instruments, vol. 
33, Dec. 1956, pp. 507-509. 

Design of Total Temperature Probes, 
by D. H. Henshaw and D. F. Daw, 
Canada, N AE LR-184, Jan. 1957. 

Edge Effects in Electromagnetic Flow- 
meters, by J. A. Shercliff, J. Nuclear 
Energy, vol. 3, no. 4, Nov. 1956, pp. 
305-311. 

Measurements of Heat Transfer from 
Fine Wires in Supersonic Flows, by John 
Laufer and Robert McClellan, Calif. 
Inst. Tech., Jet Propulsion Lab., Rep. 
20-101, June 29, 1957, 26 pp., 10 ref. 

Recent Developments in Electronic 
Instrumentation, by P. G. Hubbard, 
reprinted from Hydraulics Conference, 
6th, 1956 (Iowa State Univ., Studies in 
Engnq., Bull. 36), pp. 101-100, 16 ref. 

Potential-flow Analogs and Computers, 
by S. C. Ling, reprinted from Hydraulics 
Conference, 6th, 1956 (Iowa State Univ., 
Studies in Engng., Bull. 36), pp. 167-186, 
10 ref. 

All Metal Apparatus for the Determina- 
tion of the Density of Liquefied Gases, 
by Roger L. Jarry, Rev. Sci. Instruments, 
vol. 28, Aug. 1957, pp. 641-642. 

Effects of Thermal Relaxation and 
Specific-heat Changes on Measurements 
with a Pneumatic-Probe Pyrometer, by 
P. W. Kuhns, NACA TN 4026, July 
1957, 67 pp. 

Temperature Transducers for High 
Speed Gases, by 8S. G. Darevskii, Gt. 
Brit., Royal Aircraft Estab., Translation 
598 (ASTIA AD 121650), (Translation of 
USSR Akademia Nauk Doklady, vol. 107, 
no. 3, 1956, pp. 373-376), Oct. 1956, 5 
pp., 4 fig. 

The Development of a Mass Flow Probe, 
by E. R. Wolff, Univ. Toronto, Inst. Aero- 
phys., TN 12, Jan. 1957, 18 pp. 

Automatic Data Plotter for F-M/F-M 
Telemetering, by H. B. Riblet, Electronics, 
vol. 30, Aug. 1, 1957, pp. 182-187. 

Note on Strain Gauge Recording Equip- 
ment for the RAE Intermittent Supersonic 
Wind Tunnels, by K. G. Winter and 
E. J. Petherick, N. Atlantic Treaty Organ- 
ization, Advisory Group Aeron. Res. 
Dev., Rep. 64, Aug. 1956, 11 pp. 

Accuracy of CEP Estimates, by M. L. 
Eaton, Naval Air Missile Test Center, 
Tech. Mem. Rep. 105, Jan. 1957, 49 pp. 

A Generalization of the Nyquist Stability 
Criterion with Particular Reference to 
Phasing Error, by R. H. Merson, Gt. 
Brit., Aeron. Res. Council, Current Paper 
330 (formerly ARC Tech. Rep. 18747; 
Royal Aircraft Estab. TN GW316), 
1957, 8 pp., 1 fig. 

Determination of Equilibrium Gas 
Emissivities from Spectroscopic Data, 
by S. S. Penner and A. Thomson, Calif. 
Inst. Tech., Daniel and Florence Guggen- 


heim Jet Propulsion Center, Tech. Rep. 
25, July 1957, 63 pp. 

Microflowmeters at a Loss of D.F.G. 
Charge, by F. Bebeauvais, R. Feidt 
and C. Gebel, France, Ministere de I’ Air, 
Pub. Sci. Tech. Notes Techniques 64, 1957, 
23 pp. (in French). 

Response of Mechanical Systems to 
Random Excitation, by W. T. Thomson 
and M. V. Barton, J. Appl. Mech., vol. 
24, June 1957, pp. 248-251. 

Research and Development Trends, 
pp. H-3-4 (Research Testing-Develop- 
ment Testing-Flight Testing-Service Test- 
ing-the Future), Airborne FM-FM Tele- 
metering Units (Table), pp. H-5-9, Air- 
borne Data Recorders (Table), pp. H-10- 
17, Environmental Testing Facilities (List), 
pp. H-18, 22-36, Aviation Age, Research 
and Development Technical Handbook, 
1957-1958: Section H: Testing. 

Research and Development Trends, 
p. I-3, Computers for Aviation Research 
and Development (Table), pp. I-4-11, 
Data Review, pp. I-12-16, Aviation Age, 
Research and Development Technical Hand- 
book, 1957-1958, Section I: Computers 
and Data Processing. 


Guidance Systems and 
Components 


Theoretical Investigation Based on 
Experimental Frequency-Response Meas- 
urements of an Automatic Altitude Con- 
trol in Combination with a Supersonic 
Missile Configuration, by Ernest C. Sea- 
berg, Edward S. Geller and William W. 
Willoughby, NACA RM L54F04, Aug. 
1954, 28 pp. (Declassified from Confiden- 
tial by authority of NACA Res. Abs. 
116, p. 22, 7/5/57.) 

Principles of Self-contained Naviga- 
tion, by H. H. Bailey, Aeron. Engng. 
Rev., vol. 16, Aug. 1957, pp. 68-73. 

Cytac—A Multi-purpose Precision Navi- 
gation System, by Wilbert P. Frantz, 
Sperry Engng. Rev., vol. 10, May-June 
1957, pp. 2-6. 

Supersensitive Giant Antennas Will 
Track Long Range Missiles, by A. J. 
Garon, J. K. Van Hock and C. H. New, 
Avtation Age, vol. 28, Aug. 1957, pp. 
116-119. 

Hughes Latest Falcon Has Infra-red 
Seeker, Missiles and Rockets, vol. 2, Aug. 
1957, p. 52. 

Project Vanguard Report No. 18, 
Minitrack Report No. 1—Phase Measure- 
ment, by C. A. Schroeder, C. H. Looney, 
Jr., and H. E. Carpenter, Jr., Naval Res. 
Lab., NRL Rep. 4995, July 26, 1957, 28 
pp. 


Flight Vehicle Design and 
Testing 


Missile Progress by the Corps of Engi- 
neers, by Maj. Gen. David H. Tulley, 
USA, Missiles and Rockets, vol. 2, July 
1957, pp. 89-90. 

Progress in Army Missile Research and 
Development, by Maj. Gen. August Schom- 
burg, USA, Missiles and Rockets, vol. 2, 
July 1957, pp. 91-94. 

Navy Rockets and the Future of Man, 
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@ Motor Cases 
Solid and Liquid 
propellants 


@ Jato Cases 

@ Nozzles 

@ Plenum Chambers 
@ Blast Tubes 


@ Fuel Injectors 


ee e we have 
the “Know-How:” 


We have developed new tech- 
niques, new methods, new 
processes that effect production 
economy so necessary to a success- 
ful missile program. 


Here at Newbrook you will find 
men with experience gained 
from doing ...a modern plant 
with up-to-date equipment... 
precision inspection to meet your 
most exacting quality control 
requirements. 


And most important, Newbrook 
specialization results in strict 
reliability ! Let us help you with 
your Missile Hardware problems. 


Phone or Write 


NEWBROOK MACHINE CORPORATION 
20 Mechanic Street 
SILVER CREEK, NEW YORK 


Phone 45 


by Rear Adm. Rawson Bennett, USN, 
Missiles and Rockets, vol. 2, July 1957, 
pp. 97-98. 

Navy Missile Progress, by Rear Adm. 
John E. Clark, USN, Missiles and Rockets, 
vol. 2, July 1957, pp. 101-102. 

Progress in Air Force Scientific Research, 
by Brig. Gen. Hollingsworth F. Gregory, 
AFOSR, Missiles and Rockets, vol. 2, 
July 1957, pp. 105-109. 

Missile Encyclopedia, Missiles and 
Rockets, vol. 2, July 1957, pp. 123-164. 
US Army: Nike Ajax, Hawk, Nike 
Hercules, Dart, LaCrosse, Little John, 
Honest John, Corporal, Sergeant, Red- 
stone, Jupiter, pp. 124-136. US Navy: 
Bullpup, Corvus, Diamondback, Side- 
winder, Zuni, Sparrow, Torpedoes, Terrier- 
Tartar, Regulus I, Talos, Triton, Regulus 
II, Polaris, pp. 137-150. US Air Force: 
Mighty Mouse, Falcon, Bomarec, Rascal, 
Matador, Snark, Navaho, Thor, Atlas, 
Titan, Miscellaneous Missiles, pp. 151- 
164. 


Rockets and Satellites and Cranfield, 
Flight, vol. 72, Aug. 2, 1957, pp. 153-156. 

Trident Blents Simplicity, Performance, 
by David A. Anderton, Aviation Weei-, 
vol. 67, Aug. 19, 1957, pp. 50-59. 

RTV-N-12a Viking Rockets 8-10, Design 
Summary, Martin Co., Rep. ER-6534, 
Aug. 1955, 199 pp., 78 ref. 

Development of a High-velocity Free- 
flight Launcher—the Ames Light-gas Gun, 
by A. C. Charters, B. Pat Denardo and 
Vernon J. Rossow, NACA RM A55G11, 
Dec. 1955, 97 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abs. 118, p. 6, 8/16/57.) 

SBAC Display Suggests British Missile 
Policy Is Premature, by David A. Ander- 
ton, Aviation Week, vol. 67, Sept. 9, 
1957, pp. 26-30. 

Some Observations on Guided Missiles, 
by Th. von Karman, Interavia, vol. 12, 

Missiles at the Paris Show, by M. 
Decker, Interavia, vol. 12, Aug. 1957, pp. 
792-794. 

The West’s Guided Missiles, [nteravia, 
vol. 12, Aug. 1957, pp. 795-798. 

Air-to-Air Rocket Weapons, by F. I. 
Ordway, Interavia, vol. 12, Aug. 1957, pp. 
806-807. 

Britain’s Guided Weapons, by D. Wood, 
Interavia, vol. 12, Aug. 1957, pp. 808- 
809. 


High-Velocity, Light-Gas Gun, by W. D. 
Crozier and William Hume, J. Appl. 
Phys., vol. 28, Aug. 1957, pp. 892-894. 

Research and Development Trends, pp. 
A-3-4 (Aircraft Systems-Ballistic Missile 
Systems-The Future), Weapon System 
Programs Need Better Organization, by 
Milton W. Mohr, pp. A-5-7, Missiles 
(Table), pp. A-8-13, Aviation Age, Res. 
and Dev. Tech. Handbook, 1957-1958: 
Section A: System Engineering. 

Guided Missile Progress, by John H. 
Sides, Ordnance, vol. 42, July-Aug. 1957, 
pp. 44-46. 


Combustion, Fuels and 
Propellants 


Bibliography on Aerosols, 1950-1955, 
Supplement to Report SO-1003, by W. 
J. Scheffy, Atomic Energy Comm., COO- 
1016, June 1, 1956, 109 pp. 

Burning of a Liquid Droplet, III. 
Conductive Heat Transfer within the 
Condensed Phase During Combustion, 
by Henry Wise and C. M. Ablow, J 
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USN, Chem. Phys., vol. 27, Aug. 1957, pp. 389- 
1957, 393. 
Heats of Combustion of Some Peroxides | Ow 4 Re 
Adm. and the Heats of Formation of Acetate, 
ockels, Propionate, and Butyrate Radicals, by 


L. Jaffe, IX. J. Prosen and M. Szware, 
416-420. | 


gory, 
ol. 2 Hydrogen Bonding in Radical Reactions, | 


Some Remarks on CF; Radicals, by M. 
and Szware and H. Smid, J. Chem. Phys., | ces the TC-106 AT C. A. E. 


164. vol. 27, Aug. 1957, pp. 421-422. 

Nike Free Energy Functions for 54 Gaseous 
John, Elements, by Harwood G. Kolsky, Ruth PORTABLE a 
Red- M. Gilmer and Pau! W. Gilles, J. Chem. | 


avy: Phys., vol. 27, Aug. 1957, pp. 494— | STARTING UNIT 


On the Existence of Higher Than 


rrier- Normal Detonation Pressures, by Donald | 
gulus R. White, J. Fluid Mech., vol. 2, July | for 
orce: 1957, pp. 513-514. 


emmy Effect of Plastic Viscosity and Yield | LARGE JET 


Value on Spray Characteristics of Mag- 


151- 4 

nesium-slurry Fuel, by George M. Prok, | 

NACA RM E56J19a, Jan. 1957, 23 pp. | | AIRCRAFT 
field, (Declassified from Confidential by author- 


156. ity of NACA Res. Abs. 116, p. 20, 7/5/57.) 


Bee Experimental Performance of Chlorine 
Trifluoride-hydrazine Propellant Combi- 
nation in 100 pound Thrust Rocket Engine, 
sign by Paul K. Ordin and Riley O. Miller, 
NACA RM E9FO1, Aug. 1949, 22 pp. | 
(Declassified from Confidential by author- | 
7 ity of NACA Res. Abs. 116, p. 14, 7/5/57.) | 
Gun, Status of Combustion Research on 
and High-energy Fuels for Ram Jets, by Walter 
G11, T. Olson and Louis C. Gibbons, NACA | 
from RM ¥:51D23, Oct. 1951, 75 pp. (De-! 
Res. classified from Confidential by authority 
of NACA Res. Abs. 116, p. 14, 7/5/57.) 
ssile Preparation and Physical Properties of MODEL 141 
der- Metal Slurry Fuel, by James B. Gibbs and | TURBO-COMPRESSOR 
Preston N. Cook, Jr., NACA RM| ENGINE 
March 1952, 36 pp. (Declassi- 
iles, fied from Confidential by authority of 
12, NACA Res. Abs. 16, p. 14, 7/5/57.) 
Study of the Physical Properties of 
M. Petrolatum-stabilized Magnesium-hydro- 
pp. carbon Slurry Fuels, by Murray L. Pinns 
and Irving A. Goodman, NACA RM, 
253316, Jan. 1954, 58 pp. (Declassified | 
= from Confidential by authority of NACA | Typical of the fine results of Continental development is 
the TC-106 portable starting unit for large jet aircraft. This 
A Preliminary Study of the Preparation 
pp. of Slurry Fuels from Vaporized Magne- ae advanced new model, with a high performance turbine 
sium, by Walter R. Witzke, George M. thi 
yod, Prok Thomas J. Walsh, NACA RM as its heart, weighs third less than 
08- E53K23, Feb. 1954, 21 pp. (Declassified predecessor, yet has 17 per cent higher output, and in 
from Confidential by authority of NACA | 
. Res. Abs. 116, p. 16, 7/5/37) | addition, other important qualities: greater mobility, less 
ppl. Study of Some Dielectric Properties of neise, and a completely automatic control system. ... It 
in is now in volume production at the Continental Aviation 
pp. inera il, by Aubrey P. Altshuller, | P P 
sile NACA RM E53L09, Feb. 1954, 42 pp. | and Engineering Toledo plant. 
em (Declassified from Confidential by author- | e 
by ity of NACA Res. Abs. 116, p. 17, 7/5/57.) | 
The Magnesium of | 
Various Equivalent Diameters on Some 
58: Physical Properties of Petrolatum-Stabil- | C.A.E. gas the the J69- 
ized Magnesium-Hydrocarbon Slurries, T-19A are being built for Cessna’s T-37A twin jet trainer, Temco's 
H. by Joseph M. Lamberti, NACA RM | TT-1 Navy jet trainer, the Beech jet Mentor trainer, and the Ryan 
57, April 1954, 48 pp. _ (Declassified | 


from Confidential by authority of NACA 
Res. Abs. 116, p. 17, 7/5/57.) 


Space Flight 
55, Satellite, by Carl Gazley, Jr. ‘. god 
‘. J. Masson, Aviation Age, vol. 28, 1 
CONTINENTAL AVIATION & ENGINEERING CORPORATION 
The Atmosphere of Mars, by R. M.| im a 12700 KERCHEVAL AVENUE, DETROIT 15, MICHIGAN = 
IT. Goody, J. Brit. Interplanet. Soc., vol. 16, | 
he April-June 1957, pp. 69-83. 
‘‘Electro-gravitics’’: What It Is—or, 
J. Might Be, by A. V. Cleaver, J. Brit. | 
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Interplanet. Soc., vol. 16, April-June 1957, 
pp. 84-94. 

Soviet IGY Rocket Program Revealed, 
Missiles and Rockets, vol. 2, Aug. 1957, 
pp. 40-41. 

Vanguard Engines Pass Qualification 
Test, Missiles and Rockets, vol. 2, Aug. 
1957, p. 50. 

Space Travel and Our Technological 
Revolution, by Wernher von Braun, 
Missiles and Rockets, vol. 2, July 1957, p. 
75. 


Meteoric Dust Erosion Problem and 
Its Effect on the Earth Satellite, by Stuart 
A. Hoenig, Aeron. Engng. Rev., vol. 16, 
July 1957, pp. 34-38. ~ 

Spiral and Elliptical Orbits, by Brian 
Egerton, J. Royal Aeron. Soc., vol. 61, 
June 1957, pp. 422-423. 

Algae Hazard Found Facing Space 
Travel, Aviation Week, vol. 67, Sept. 9, 
1957, p. 65. 


Astrophysics, Aerophysics 


A Consideration of Radio Star Scintilla- 
tions as Caused by Interstellar Particles 
Entering the Ionosphere, Part III: The 
Kind, Number, and Apparent Radiant of 
the Incoming Particles, by G. A. Harrower, 
Canadian J. Phys., vol. 35, July 1957, pp. 
792-798. 

The Study of the Upper Atmosphere by 
Means of Rockets, at the Academy of 
Sciences, USSR, by S. M. Poloskov and 
B. A. Mirtov, J. British Interplanet. Soc., 
vol. 16, April-June, 1957, pp. 95-100. 


Nuclear Propulsion 


. Promising Nuclear Powerplants, by 
Leo Seren, Aviation Age, vol. 28, July 
1957, pp. 42-49. 

Dose Rates Produced from Gamma Ray 
Sources, by J. G. Lewis, Ind. Engng. 
Chem., vol. 49, July 1957, pp. 1197-1202. 

Oscillations of a Cylindrical Plasma, 
by Thomas H. Stix, Physical Rev., vol. 
106, June 15, 1957, pp. 1146-1150. 

The Lithium Hydride, Deuteride, and 
Tritide Systems, by F. K. Heumann and 
O. N. Salmon, Atomic Energy Comm., 
KAPL-1667, Dec. 1, 1956, 54 pp., 25 ref. 

Decomposition and Recombination in 
Irradiated Static Water Systems at 
High Temperatures, by J. R. Humphreys, 
E. K. Abers and Y. Solomon, Atomic 
Energy Comm., ANL-5004, March 1953 
(declassified Feb. 23, 1957), 55 pp. 

High-energy Particle Data, vol. 1, by 
Beverly Hill Willis and Charles V. Stable- 
ford, Atomic Energy Comm., UCRL 2426, 
rev. Nov. 1956, 87 pp. 

Criticality Calculations for Hydrogenous 
Systems, by J. R. Knight, Atomic Energy 
Comm., K-1260, Nov. 25, 1955 (declassi- 
fied April 22, 1957), 20 pp., 14 ref. 

Neutron Shielding Materials, by T. 
M. Snyder, Atomic Energy Comm., KAPL- 
704, Feb. 12, 1952 (declassified April 30, 
1957), 6 pp. 

Investigation of Burnout Heat Flux in 
Rectangular Channels at 2000 Psia, by 
H. S. Jacket, J. D. Roarty and J. E. 
Zerbe, Atomic Energy Comm., WAPD- 
AIW(IM)-3, Dee. 5, 1955, 52 pp. 

Reactor Heat Transfer Progress, No. 
13, ed. by John E. Viscardi (Nuclear Dev. 
Corp. of America), Atomic Energy Comm., 
NDA-33 (2008-18), Feb. 14, 1957, 27 pp. 

Conference on the Use of Organic 
Materials in Nuclear Reactors, Aug. 5-6, 
1953, Downey, Calif. (North American 
Aviation, Atomic Energy Res. Dept.) 


Atomic Energy Comm., TID-10041, Dec. 
1, 1953 (declassified April 23, 1957), 
146 pp. 

Nuclear Ramjet Study, Aviation Week, 
vol. 67, Aug. 19, 1957, p. 27. 

Nuclear Heat Won’t Be Used Directly 
in Near Future, by A. A. Kucher, T. F. 
Nagey and H. J. Ogorzaly, SAE J., 
vol. 65, Aug. 1957, pp. 24-25. 

Heat Generation by Inelastic Scatter- 
ing Gamma Rays, by J. P. Davis, Comb. 
Engng., Inc., Reactor Dev. Div. (U.S. 
Atomic Energy Comm.), GERDSIC-108, 
April 1956, 8 pp. 

Liquid Metals and Nuclear Power, 
by C. D. Boadle, Atomics and Nuclear 
Energy, vol. 8, Feb. 1957, pp. 41-45, 
March 1957, p. 83. 


Aerodynamics of Jet 
Propelled Missiles 


Performance Comparison of Three 
Canard-Type Ram-jet Missile Configura- 
tions at Mach Numbers from 1.5 to 2.0, 
by Evan A. Fradenburgh and Emil J. 

emzier, NACA RM E53F11, Aug. 1953, 
31 pp. (Declassified from Confidential 
by authority of NACA Res. Abs. 116, 
p. 16, 7/5/57.) 

A Method for Simulating the Atmos- 
heric Entry of Long-Range Ballistic 

issiles, by A. J. Eggers, Jr., NACA RM 
A55115, Dec. 1955, 19 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abs. 116, p. 14, 7/5/57.) 

Flight Measurements of Boundary- 
Layer Temperature Profiles on a Body of 
Revolution (NACA RM-10) at Mach 
Numbers from 1.2 to 3.5, by Andrew G. 
Swanson, James J. Buglia, and Leo T. 
Chauvin, NACA TN 4061, July 1957, 
40 pp. 

An Experimental Investigation of the 
Flow over Simple Two Dimensional and 
Axial Symmetric Bodies at Hypersonic 
Speeds, by I. E. Vas, S. M. Bogdonoff 
and A. G. Hammitt, Princeton Univ., 
Dept. Aeron. Engng., Rep. 382 (Wright 
Air Dev. Center, TN 57-246; ASTIA 
AD 130826), June 1957, 11 pp., 9 fig. 

Two-dimensional, Steady, Cavity Flow 
about Slender Bodies, by Hirsh Cohen 
and Robert Gilbert, J. Applied Mech., 
vol. 24, June 1957, pp. 170-176. 


Research and Development Trends 
Aviation Age, Research and Development 
Technical Handbook, 1957-1958, Section 
B. Aerodynamics and Structures. pp. 
B-3-5. (Space Dynamics Aerodynamics- 
Structures. ) 

Calculating Thermal Stresses in Sand- 
wich Panels, by W. C. Schoeller, Avia- 
tion Age, Research and Development 
Technical Handbook, 1957-1958, Section 
B. Aerodynamics and Structures, pp. 
B-6-8. 

Data Review, Aviation Age, Research 
and Development Technical Handbook, 
1957-1958, Section B. Aerodynamics and 
Structures, pp. B-10, 12. 

Analysis of Inequalities of Flow Field, 
by L. Frager, Fusées et Recherche Aeron., 
vol. 2, May 1957, pp. 133-134 (in French). 

The Prediction of Nonlinear Pitching 
and Yawing Motion of Symmetric Missiles, 
by Charles H. Murphy, Jr., J. Aeron. 
Sci., vol. 24, July 1957, pp. 473-479. 

Optimum Nose Shapes for Missiles in 
the Superaerodynamic Region, by W. J. 
Carter, J. Aeron. Sci., vol. 24, July 1957, 
pp. 527-532. 

Two-Dimensional Bow Shock Wave 
Detachment Distances, by R. A. A. 
Bryant and J. N. G. Grant, J. Roy. 


Aeron. Soc., vol. 61, June 1957, pp. 424— 
426. 

The Effect of a Spike Protruding in 
Front of a Bluff Body at Supersonic 
Speeds, by D. Beastall and J. Turner, 
Gt. Brit., Aeron. Res. Council, Rep. and 
Mem. 3007 (formerly ARC Tech. Rep. 
14892; Royal Ai:craft Estab., TN Aero 
2137) 1957, 26 pp. 

Measurement of Projectile Motion 
during Engraving Process, by J. A. Peter- 
sam, H. G. McGuire, and T. E. Turner, 
Aberdeen Proving Ground, Ballistic kes. 
Labs., Mem. Rep. 1070, May 1957, 53 pp. 


An Introduction to the Equations of 
Magnetogasdynamics, by J. A. Steketve, 
University of Toronto, Inst. Aerophys., 1'N 
13, May 1957, 13 pp. 

Investigation of Some Wake Vortex 
Characteristics of an Inclined Ogive-Cylin- 
der Body at Mach Number 1.98, by Leland 
H. Jorgensen and Edward W. Perkins, 
NACA RM A55E31, Aug. 1955, 47 pp. 
(Declassified from Confidential by authwr- 
ity of NACA Res. Abs. 118, p. 6, 8/16/57.) 


A Flight Investigation of Mach Num- 
bers from 0.67 to 1.81 of the Longitudinal 
Stability and Control Characteristics of 
a 60° Delta-Wing Missile Configuration 
Having an All-Movable Tail, by Martin 
T. Moul and Hal T. Baber, Jr., NACA 
RM _ 153G29, October 1953, 43 pp. 
(Declassified from Confidential by author- 
ity of NACA Res. Abs. 118, p. 12, 8/16/57.) 


A Pressure-Distribution Investigation 
of a Fineness-Ratio-12.2 Parabolic Body 
of Revolution (NACA RM-10) at M = 
1.59 and Angles of Attack up to 36°, 
by Morton Cooper, John P. Gapcynski 
and Lowell E. Hasel, NACA RM 
L52G14a, Oct. 1952, 89 pp. (Declassi- 
fied from Confidential by authority of 
NACA Res. Abs..118, p. 11, 8/16/57.) 


Investigation of Flow Around Simple 
Bodies in Hypersonic Flow, by Toshi 
Kubota, Calif. Inst. Tech., Hypersonic 
Res. Project, Mem. 40, June 25, 1957, 91 


Bodies of Revolution Having Minimum 
Drag at High Supersonic Speeds, by A. J. 
Eggers, Jr., Meyer M. Resnikoff, and 
David H. Dennis, NACA Rep. 1306, 
1957, 12 pp. (supersedes NACA TN 3666). 


Performance Characteristics of Canard- 
Type Missile with Vertically Mounted 
Nacelle Engines at Mach Numbers 1.5 to 
2.0, by Leonard J. Okery and Howard §. 
Krasnow, NACA RM _ E52H08, Sept. 
1952, 25 pp. (Declassified from Confi- 
dential by authority of NACA Res. Abs. 
119, p. 18, 8/16/57.) 

Free-Flight Tests to Determine the 
Power-on and Power-off Pressure Distri- 
bution and Drag of the NACA RM-10 
Research Vehicle at Large Reynolds 
Numbers between Mach Numbers 0.8 
and 3.0, by Sherwood Hoffman, NACA 
RM 1L55H02, Sept. 1955, 55 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abs. 119, p. 23, 8/16/57.) 


A Theoretical Analysis of a Simple 
Aerodynamic Device to Improve the 
Longitudinal Damping of a Cruciform 
Missile Configuration at Supersonic 
Speeds, by James E. Clements, NACA RM 
L55H31, Oct. 1955, 36 pp. (Declassi- 
fied from Confidential by authority of 
NACA Res. Abs. 119, p. 23, 8/16/57.) 


Experimental Investigation of Thick, 
Axially Symmetric Boundary Layers on 
Cylinders at Subsonic and Hypersonic 
Speeds, by Ronald L. Richmond, Calif. 
Inst. Tech., Hypersonic Res. Project, Mem. 
39, June 20, 1957, 66 pp. 

Measurement of the Upwash Factor 
on a Cone-Cylinder Model, by Oscar C. 
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Holder, Army Ballistic Missile Agency, 
Rep. DA-R-11, June 20, 1957, 23 pp. 

A Correlation of Free-flight Transition 
Measurements on Various Blunt Nose 
Shapes by Use of the Momentum-thick- 
ness Reynolds Number, by W. R. Witt, 
Jr., and J. Persh, (Aeroballistic Range 
Symposium, Proc., Jan., 1957), Aberdeen 
Proving Ground, Ballistics Res. Lab., Rep. 
1005, Pt. 1, March 1957, pp. 23-43. 


Jet and Rocket 
Propulsion Engines 


Experimental Investigation of an Axial- 
Flow Supersonic Compressor Having 
Rounded Leading-Edge Blades with an 
8-per cent Mean Thickness-Chord Ratio, 
by Theodore J. Goldberg, Emanuel 
Boxer and Peter T. Bernot, NACA RM 
L53G16, Dec. 1953, 43 pp. (Declassified 
from Confidential by aahesity of NACA 
Res. Abs. 118, p. 8, 8/16/57.) 

Effects on the Weight-Flow Range and 
Efficiency of a Typical Axial-Flow Com- 

ressor Inlet Stage That Result from the 

se of a Decreased Blade Camber or 
Decreased Guide-Vane Turning, by Robert 
J. Jackson, NACA RM E52G02, Sept. 
1952, 41 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
18, p. 8, 8/16/57.) 

Experimental Investigation of a 16- 
inch Impulse-Type Supersonic-Compres- 
sor Rotor, by Guy N. Ullman, Melvin J. 
Hartmann and Edward R. Tysl, NACA 
RM E51G19, Oct. 1951, 29 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abs. 118, p. 7, 8/16/57.) 

The 20,000-lb. s.t. de Havilland Gyron 
D.Gy.2, Aeroplane, vol. 93, Aug. 2, 
1957, pp. 154-158. 

A Theoretical Investigation of the 
Effects of Changes of Shape on the 
Operation of the R.P.I. Wave Engine, 
by R. S. Rainhill, Rensselaer Polytech. 
Inst., Res. Inst., TR-AE-5703, June 1957, 
33 pp., 24 fig., 3 tab. 

A Method for Eliminating the Effects 
of Residual Nitric Acid in Rocket Engines 
and Related Equipment, by Edward L. 
Harris, Wright Air Dev. Center, Tech. 
Rep. 57-31 (ASTIA AD 1107327), Dec. 
1956, 6 pp. 

Screen-Type Noise Reduction Devices 
for Ground Running of Turbojet Engines, 
by Willard D. Coles and Warren J. North, 
NACA TN 4033, July 1957, 23 pp. 

Foreign-Object Retention and Flow 
Characteristics of Retractable Engine- 
Inlet Screens, by Fred W. Steffen and 
Lewis A. Rodert, NACA RM E57A15, 
July 1957, 30 pp. 

Modern Developments in Solid-Pro- 
pellant Rocket Engineering, by Richard 
D. Geckler and Robert E. Davis, Aeron. 
Engng. Rev., vol. 16, Aug. 1957, pp. 
42-46, 51. 

Progress in Air Force Rocketry, by 
Frederick I. Ordway III and Ronald C. 
Wakeford, Missiles and Rockets, vol. 2, 
July 1957, pp. 83-87. 

Propulsion Progress, Missiles and Rock- 
ets, vol. 2, July 1957, pp. 81. 

Investigation of Supersonic-Compres- 
sor Rotors Designed with External Com- 
pression, by Lawrence J. Jahnsen and 
Melvin J. Hartmann, NACA RM 
£54G27a, Sept. 1954, 41 pp. (Declassi- 
fied from Confidential by authority of 
NACA Res. Abs. 118, p. 9, 8/16/57.) 

Small Solid Rockets for Commercial 
Use, by G. E. Rice, Missiles and Rockets, 
vol. 2, Aug. 1957, pp. 80-81. 

Solid Propellant Rocket Testing, by 
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Men who know 


TURBINES, PUMPS, COMPRESSORS 


Get in on the development of the most 
advanced high-speed rotary equipment ever built 


Here’s your opportunity to step up to a new, higher level of 
turbo-machinery technology—the Large Rocket Engine. 

Whatever type of rotating machinery you know best, your 
experience could be extremely valuable in the important develop- 
ments now going on at Rocketdyne. The seasoned and ambitious 
man who has cut his teeth on jet engines, steam or gas turbines, 
or other elements of rotating machinery, is urgently needed to 
apply mechanical principles to meet the increasing demands of 
power plant performance. 

The combination of high speed, light weight, heavy loadings 
and exceptional pressures required in rocket engine work is lead- 
ing to an entirely new breed of high-performance rotating 
machinery ...and a new breed of engineer. You can be one of 
this advance guard of the turbo-machinery field—if you have the 
desire to build your professional status by accepting new 
challenges. 

Rocketdyne is building high-thrust rocket engines for the 
nation’s major missiles. You'll work with the leading producer in 
the nation’s fastest growing industry. You and your work will 
be recognized as a vital part of the overall achievements. Testing 
facilities are among the world’s finest. The power produced is 
beyond anything ever before thought possible. If you would like 
to tackle new assignments working alongside some of the finest 
minds in turbo-engineering, write and tell us about your back- 
ground: A. W. Jamieson, Rocketdyne Engineering Personnel 
Dept. R-3, 6633 Canoga Avenue, Canoga Park, California. 


ROCKETDYNE 


A DIVISION OF NORTH AMERICAN AVIATION, INC. 


BUILDERS OF POWER FOR OUTER SPACE 
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H. E. Westgate, Missiles and Rockels, 
vol. 2, Aug. 1957, pp. 87-90. 

Determination of Surge and Stall 
Limits of an Axial-Flow Turbojet Engine 
for Control Applications, by Ross D. 
Schmidt, George Vasu, and Edward W. 
McGraw, NACA TN 3585, Sept. 1957, 
29 pp. (supersedes RM E53B10). 


Visualization of Rotor Tip Secondary 
Flows with Blade Tip. Air Discharge and 
Suction in a Low-Speed Turbine, by Mil- 
ton G. Kofskey and Hubert W. Allen, 
NACA RM E56E16, August 1956, 28 
pp. (Declassified from Confidential by 
authority of NACA Res. Abs. 119, p. 20, 
8/16/57.) 

Aerodynamic Forces on a Vibrating 
Unstaggered Cascade (Luftkrifte an einem 
schwingenden Gitter), by H. Sdhngen, 
NACA TM 1412, Aug. 1957, 16 pp. 
(translation from ZAMM vol. 35, March 
1955, pp. 81-88). 

Preliminary Investigation of Hollow- 
Bladed Turbines Having Closed and Open 
Blade Tips, by Gordon T. Smith and 
Robert O. Hickel, NACA RM E55F27a, 
Aug. 1955. 18 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abs. 119, p. 20, 8/16/57.) 

Performance of a Supersonic Mixed- 
Flow Rotor with a Swept Leading Edge 
and 0.52 Inlet Radius Ratio, by Arthur W. 
Goldstein and Ralph L. Schacht, NACA 
RM E53H27, Nov. 1953, 34 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abs. 119, p. 19, 8/16/57.) 


Performance of a Swept Leading Edge 
Rotor of the Supersonic Type with Mixed 
Flow, by Arthur W. Goldstein and Ralph 
L. Schacht, NACA RM _ E52K03, Jan. 
1953, 34 pp. (Declassified from Confi- 


dential by authority of NACA Res. 
Abs. 119, p. 19, 8/16/57.) 

Investigation of a 24-inch Shock-in- 
Rotor Type Supersonic Compressor De- 
signed for Simple Radial Equilibrium be- 
hind Normal Shock, by Haroid Lown and 
Melvin J. Hartmann, NACA RM 
E51H08, Dec. 1951, 25 pp. (Declassi- 
fied from Confidential by authority of 
NACA Res. Abs. 119, p. 17, 8/16/57.) 

Spectre—de Havilland’s Latest Rocket 
Motor, Flight, vol. 72, Aug. 23, 1957, p. 
245. 

The Rocket Comes of Age, by W. von 
Braun, Interavia, vol. 12, Aug. 1957, pp. 
789-790. 

Rocket Propulsion Units, by Vernet- 
Lozet, Interavia, vol. 12, Aug. 1957, pp. 
799-801. 

Solid Propellant Rocketry, by J. K. 
Elder, Interavia, vol. 12, Aug. 1957, pp. 
801-802. 

Principles of Control of Aircraft Gas 
Turbines, by L. Beck, Z. Flugwissen- 
schaften, vol. 5, July 1957, pp. 199-204. 


Heat Transfer and 
Fluid Flow 


One-dimensional and Two-dimensional 
Gas Dynamics Analogies, by R. A. A. 
Bryant, Aust. J. Appl. Sci., vol. 7, no. 4, 
Dec. 1956, pp. 296-313. 

Study ef Stagnation Temperature and 
Pressure in a Free-Molecular Flow Re- 
gime by Means of the Revolving Arm 
Method, by F. M. Deviennce. Lab. 
Mediterranean de Rech., Thermodynamiques 
TR (AFOSR TR 57-32), (AD 126495), 
March 1957, 102 pp. 


Velocity, Temperature, and Heat-trans- 
fer Measurements in a Turbulent Bound- 
ary Layer Downstream of a Stepwise 
Discontinuity in Wall Temperature, by 
D. 8. Johnson, J. Appl. Mech., vol. 24, 
March 1957, pp. 2-8. 

Effect of Adverse Pressure Gradients 
on Turbulent Boundary Layers, by J. M. 
Robertson and J. W. Holl, J. Appl. 
Mech., vol. 24, June 1957, pp. 191-196. 


Film Condensation of Vapor in Hori- 
zontal Tube, by J. G. Chaddock, Refriger. 
Engng., vol. 65, April 1957, pp. 36-41, 
90, 92-95. 

Heat Transfer in Condensation of Steam 
Flowing in a Horizontal Group of Tubes, 
by S. N. Fuks, Jeploenergetika, vol. 4, 
Jan. 1957, pp. 35-38 (in Russian). 

Study of Boiling in Short Narrow Chan- 
nels and Its Application to Design of Mag- 
nets Cooled by Liquid H:2 and Nz, by S. ‘:. 
Sydoriak and T. R. Roberts, J. Appl. 
Phus.. vol. 28, Feb. 1957, pp. 143-148. 


Heat Transfer Media, by P. M. Krishna 
and D. Venkateswarlu, Chem. and Process 
Engng., vol. 37, Nov. 1956, pp. 382-388. 

Production and Use of Very High Tem- 
perature, by W. Lochte-Holtgreven, (V DJ 
Zeitschrift, No. 23, 1955, pp. 785-788), 
Gr. Brit., RAE Lib. Transl. 630, Jan. 
1957. 

Some Structural Penalties Associated 
with Thermal Flight, by J. W. Mar and 
L. A. Schmit, Trans. ASME, vol. 79, 
July 1957, pp. 990-1004. 


Generation of High Temperatures (up 
to 55,000°) in the Laboratory, by . 
Preining (Ost. Chemiker-Zeit., no. 5/6, 
1957, pp. 67-72), Gt. Brit. RAE Lib. 
Transl., 629, Feb. 1957, 11 pp., 20 refs. 
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CHARACTERISTICS 


ANALYSIS 


Stainless Steel 
Ball and Race 


Chrome Alloy { For types operating under a radial 
ultimate loads (3000-893,000 Ibs.). 


For types operating under normal loads 
with minimum friction requirements. 


Steel Ball and Race 


Bronze Race and 
Chrome Steel Bail 


Thousands in use. Backed by years of service life. Wide variety 
of Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in 
similar size range with externally or internally threaded shanks. 
Our Engineers welcome an opportunity of studying individual 
requirements and prescribing a type or types which will serve 
under your demanding conditions. Southwest can design special 
types to fit individual specifications. As a result of thorough 
study of different operating conditions, various steel alloys 
have been used to meet specific needs. Write for revised Engi- 
neering Manual describing complete line. Dept. JP-58. 


SOUTHWEST PRODUCTS CO. 


- 1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 
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RECOMMENDED USE 


For types operating under high temper- 
ature (800-1200 degrees F.). 


BOVANCED DESIGN TRANSDUCER 


NEW Self-Powered 
Speed Sensing Pickup | 
operates relays directly’ 


| Ferrous metals passing pickup operate 
| relay without amplifier 
No tubes or transistors needed with = 
the Powr-Mite #3040 self powered “yo 
| magnetic pickup. 
Simply feed the Powr-Mite #3040 » 
| output into any relay having 100 milli- 
watts sensitivity. : 

You get extra dependability through ~~ 
circuit simplification, plus important * 
savings in weight and space. 

Ideal for use with over and under 


speed controls. 


*100 milliwatts of power at surface speeds of 250 in./sec. 


Electro 


TRANSDUCERS FOR COUNTING 
CONTROL AND MEASUREMENT 


Powr-Mite #3040 $25.00 list 
@h 


ELECTRO LABORATORIES 


4501-J North Ravenswood, Chicago 40 
Canada: Atlas Radio Ltd., Toronto 
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Heat Transfer Rates for Gases, by 
F. L. Rubin, Petroleum Refiner, vol. 36, 
March 1957, pp. 226-228. 

Mean-Temperature Difference in One, 
Two, and Three-Pass Crossflow Heat 
Exchangers, by R. A. Stevens, J. Fernan- 
dez and J. R. Woolf, Trans. ASME, vol. 
79, Feb. 1957, pp. 287-297. 

On the Heat Transfer Coefficients of 
Cylindrical Bodies, II, by T. Tsubouchi, 
S. Sato, and H. Usami, Tohoku University, 
Rep. Inst. High Speed Mech., vol. 7, 

1956, pp. 31-44. 

A New Method to Measure Prandtl 
Number and Thermal Conductivity of 
Fluids, by E. R. G. Eckert and T. F. 
Irvine, Jr., J. Appl. Mech., vol. 24, March 
1957, pp. 25-38. 

Predict Thermal Conductivity, II. How 
to Estimate Engineering Properties, by 
Wallace R. Gambill, J. AIChE, vol. 64, 
March 1957, pp. 271-276. 

Heat Transfer to Turbine Blades, by 
S. J. Andrews and P. C. Bradley, Gt. 
Brit., Aeron. Res. Council, Current Paper 
294, (formerly ARC Tech. Rep. 12078; 
Natl. Gas Turbine Estab., Mem. M37), 
1957, 12 pp., 11 fig. 

Investigation of Vanes Immersed in 
the Jet of a Solid-Fuel Rocket Motor, 
by Leo V. Giladett and Andrew R. Wine- 
man, NACA RM _ L52F12, Sept. 1952, 
30 pp. (Declassified from Confidential 
by authority off NACA Res. Abs. 116, 
p. 21, 7/5/57.) 

A Theoretical Investigation of Shock 
Oscillations in One Dimensional Flows 
with and without Pressure Gradients, 
by Hans B. Schechter, Wright Air Dev. 
Center, T'N 57-37, Feb. 1957, 30 pp. 

On the Mixing Theory of Crocco and 
Lees and Its Application to the Interaction 
of Shock Waves and Laminar Boundary 
Layer, Part I: Generals and Formula- 
tions, by Sin I. Cheng and K. N. C. 
Bray, Princeton Univ., Dept. Aeron. 
Engng., Rep. 376 (AFOSR-TN-57-283; 
ASTIA AD 132354), May 1957, 46 pp., 
19 fig. 

Reflection and Refraction of Acoustic 
Waves by a Shock Wave, by J. Brillouin, 
NACA TM 1409, July 1957, 42 pp. 
(Translation from Acustica, vol. 5, no. 3, 
1955, pp. 149-163.) 

Study of Certain Problems of Liquid 
Gravity Waves, by Georges Brillouet, 
France, Ministere de l’Air, Pub. Sci. 
Tech. 329, 1957, 145 pp. (in French). 

One Dimensional Treatment of Non- 
uniform Flow, by R. D. Tyler, Gt. Brit. 
Aeron. Res. Council, Rep. and Mem. 
2991 (Formerly ARC Tech. Rep. 17186), 
1957, 20 pp. 

The Effect of Chemical Nature of Heat- 
ing Surfaces on the Heat Transfer Coeffi- 
cients of Pool Boiling Liquids, by George W. 
Preckshot, Atomic Energy Comm., AECU- 
3183, Dec. 1955, 61 pp. 

Heat Transfer to Lead Bismuth in Tur- 
bulent Flow in an Annulus, by R. A. 
Seban and D. F. Casey, Atomic Energy 
Comm., A ECU-3164, June 1956, 16 pp. 

Heat Transfer to Liquid Metals in Con- 
centrated Annuli, by Raymond V. Bailey, 
Atomic Energy Comm., ORNL-521, June 
1950, 45 pp. 

Heat Transfer; a Bibliography of 
Unclassified Report Literature, by James 
M. Jacobs and Gifford A. Young, Atomic 
Energy Comm., TID 3305, March 1957, 
169 pp. 

Analysis of Heat-driven Oscillations of 
Gas Flows, II: On the Mechanism of the 
Rijke-tube Phenomenon, by H. J. Merk, 
Applied Sci. Res., vol., 6, Sect. A, no. 5-6, 
1957, pp. 402-420. 
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cast 


mandrels 


or cores ? 


Aluminum mandrels for 

forming solid fuel propellant are 
now being cast in production 

by the unusual foundry methods of 
Morris Bean & Company. While 
we assume there is no present 
need for a mandrel as large as the 
one on the left, it can be cast. 
Currently we are working 

on solid and hollow mandrels 

up to 8 feet long. Their smooth 
surfaces and accurate contours 
eliminate much difficult machining; 
cost is drastically reduced. 
In addition to large size, 
we would be happy to 
explore with you ways to 
produce intricate star-lobes. 
Telephone or write. 


Morris Bean & Company, 
Yellow Springs 6, Ohio. 
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Do you” 
think of 
pressure 
transducers? 


ATOMICS 
INTERNATIONAL 


DIBION oF NORTH AMERICAM AWATION, 180. 


does: ee uses Statham 
pressure transducers as-part of the 
contfol instrumentation in the con- 
struction of its nuclear reactors 
which open entirely new fields in 
industrial research and develop- 
ment by providing: an. on-the-spot 
source of rays 


WHEN THE NEED 
IS TO KNOW...FOR SURE 
SPECIFY STATHAM 
Accelerometers 
Pressure Transducers 
Load Cells 
Catalog, complete with prices, 
available upon request. 


INSTRUMENTS, INC. 
LOS ANGELES 64 
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ELectro Propucts LABORATORIES 
Robertson, Buckley & Gotsch, Inc., Chicago, Ill. 


Forp INSTRUMENT ComPaANyY, DIVISION OF SPERRY RAND CORPORATION... 
James Thomas Chirurg Co., New York, N. Y. 


Tue CorPoRATION, AIRESEARCH MANUFACTURING CoMPANY.... 
J. Walter Thompson Co., Los Angeles, Calif. 


LockHEED ArrcraArt ComMpaANy, MIssILe Systems Division 
Hal Stebbins, Inc., Los Angeles, Calif. 


Morris Bean & Company 
Odiorne Industrial Advertising, Inc. 


NEWBROOK MaAcHINE CoRPORATION 
Melvin F. Hall Adv. Agency, Inc., Buffalo, N. Y. 


OLDsMOBILE Division, GENERAL Motors 153 
D. P. Brother & Company, Detroit, Mich. 


THE Ramo-WooLpRIDGE CORPORATION 
The McMarty Company, Los Angeles, Calif. 


Reaction Motors, Inc 
The Aitkin-Kynett Co., Philadelphia, Pa. 


RHEEM MANUFACTURING COMPANY 
Getz & Sandborg, Inc., Beverly Hills, Calif. 


RocketpynE, A Division or NortH AMERICAN AVIATION, INC 
Batten, Barton, Durstine & Osborn, Inc., Los Angeles, Calif. 


SEATON-WILSON MANUFACTURING COMPANY 


SouTHWEsT PRropucts CoMPANY 
O. K. Fagan Advertising Agency, Los Angeles, Calif. 


StatTHaM INsTRUMENTS, INC 
Compton Advertising, Inc., Los Angeles, Calif. 


JET PROPULSION 
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